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ABSTRACT 

The objective of this study was to determine the 
influence of site properties on the response of barley, 
rapeseed, and wheat to phosphate fertilizer in Alberta. 
Yield and site data from 254 field experiments from the 
period of 1969 to 1975 were assembled. Additional 
information, including site classification (agro-climatic 
area, soil zone, and soil order), and laboratory analysis of 
particle size distribution, CaCO, equivalence, and organic 
matter content of the surface depth of the field sites were 
determined. Discriminant analysis was used to determine 
those site properties important for the separation of sites 
into responsive and unresponsive categories. Multiple 
regression procedures were used to determine those site 
variables which could significantly account for the 
variation in yield increase of the responsive sites. 
Principal component analysis was used to identify the 
interrelationships among site properties of the responsive 
sites. 

Analyses of the pooled barley data for 125 site-years 
indicated that the soil test for phosphorus (ASFTL-P) was 
the most important site variable influencing site 
Separation, and for the variation in yield increase of the 
responsive sites. Clay and CaCO; content of the soil and 
growing season precipitation were additional variables which 
appeared to be important for site separation, while soil pH, 


growing season precipitation and organic matter content of 
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soils were additional co-variates that significantly 
accounted for the variation in yield increase of the 
responsive sites. Site classification had a significant 
influence on both site separation and variation in yield 
increase. Principal component analysis indicated an inverse 
relationship between the required phosphate fertilizer rate 
for "optimum" yield and each of ASFTL-P, soil pH, and 
Organic matter content of soils. 

Analyses of the pooled rapeseed data for 91 site-years 
indicated that the soil test for phosphorus (ASFTL-P) was 
the most important site parameter affecting site separation, 
and for the variation in yield increase of the responsive 
Sites. Other site variables which appeared to be important 
for site separation were soil electrical conductivity (E.C.) 
and clay content of soils, while CaCO; was the only 
additional variable to significantly account for the 
variation in yield increase. Site classification appeared to 
be important for site separation, but did not significantly 
account for any of the variation in yield increase. 
Principal component analysis revealed an inverse 
relationship between the required phosphate fertilizer rate 
for "optimum" yield and each of ASFTL-P, soil pH, organic 
matter content of soils and growing season precipitation. 

Results of the analyses of the pooled wheat data for 38 
Site-years indicated that a soil test for phosphorus was the 
most important Site variable to influence site Separation, 


and for the variation in percent yield increase of the 
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responsive sites. However, the specific soil test procedure 
varied among the results of the discriminant analyses. Other 
Site properties influencing site separation included organic 
matter content of soils, and soil E.C., while soil E.C. was 
the only additional variable to significantly account for 
the variation in percent yield increase. Site classification 
did not appear to have a clear influence on either site 
Separation or the variation in percent yield increase. 
Principal component analysis of the responsive sites 
indicated an inverse relationship between the phosphate 
fertilizer rate for "optimum" yield and each of soil pH, 
soil organic matter, and growing season precipitation, but 
the positive relationship with Olsen-P was contrary to the 
results of the barley and rapeseed sites, and the expected 
relationship. 

The results of this study lack verification with data 
external to this study, but suggest that the phosphate 
fertilizer rate for "optimum" yield should decrease as the 
soil test for phosphorus (ASFTL-P), soil pH, organic matter 
content of soils and/or growing season precipitation 


increase. 
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I. INTRODUCTION 

For optimum crop growth and nutrition, the soils of 
western Canada have been generally considered to be low in 
plant available phosphorus. Phosphate fertilizer application 
has proven to be benefical in promoting flowering, seed 
formation, root growth, disease resistance, Straw strength, 
and maturation, in addition to increasing yield. The major 
problem has been the prediction, by means of a soil test, of 
phosphate fertilizer requirements and the yield response to 
phosphate addition. 

Numerous soil test procedures for measuring the 
available phosphorus status of the soil have been developed 
and used, with varying degrees of success. In western 
Canada, studies have been conducted to determine the best 
soil test procedure to meaSure available phosphorus, predict 
crop response to phosphate fertilizer, and to determine 
optimum phosphate fertilizer requirements. In general, 
greenhouse studies have resulted in better correlations 
between soil test phosphorus and yield response than have 
field studies. Field research in Manitoba found that the 
relationship between percent yield and extractable 
phosphorus by a number of methods was not very high 
regardless of how the crops or soils were selected (Soper, 
1967). Poor correlations between percent yield for cereal 
crops and the available phosphorus were also found in field 
studies conducted in Alberta (Robertson, 1967). All three 


western Canadian prairie provinces have recognized 





differences among soils and among climatic areas in regards 
to soil test phosphorus levels and crop response to 
phosphate fertilizer. To improve predictions of phosphate 
fertilizer requirements, attempts have been made at 
developing new soil test procedures. Alternatively, as more 
researchers recognize the influence of environmental factors 
and soil properties other than the fertility status of soils 
on crop response to fertilizer, the use of more elaborate 
statistical and modelling techniques has become increasingly 
common. 

The objectives of this study were, by means of 
discriminant analysis, multiple regression procedures, and 


principal component analysis, 


1. to determine the best soil test procedure 
for predicting crop response to phosphate 
fertilizer and to aid in the prediction of 
"optimum" fertilizer rates. 

2. to determine the influence of various 
chemical and physical soil properties on 
crop response to phosphate fertilizer. 

3. to determine the value of site 
classification systems in predicting crop 


response to phosphate fertilizer. 





II. LITERATURE REVIEW 


A. INTRODUCTION 

It has long been recognized that crop yield, both 
quantity and quality, is a function of the soil on which the 
crop is grown, the climate, management factors, and the crop 
itself (Fitts, 1974). The influence of each factor is 
difficult to discern since each is a broad category 
consisting of several components, each of which may be 
modifying or limiting. The fertility status of a soil is but 
one component of the soil factor, and is composed of several 
individual elements. Thus it is difficult to predict crop 
yield from only one variable such as the available soil 
phosphorus without taking into account the other growth 
factors. These other growth factors can exert strong 
influences on yields and fertilizer effects. 

When compared with carbon, hydrogen, oxygen, and 
nitrogen, the phosphorus content of plants is small, in the 
range of 0.1 to 1.0 %, yet this element is essential for 
plant nutrition. Its most important function within the 
plant is that of energy storage and transfer as adenosine 
triphosphate (ATP) in biochemical processes such as 
photosynthesis, electron transport, active ion transport, 
and sucrose transport. In addition, it iS an important 
structural component of numerous compounds, including 
phospholipids, nucleic acids, phytin, sugar phosphates, and 


coenzymes (Glass et a]., 1980; Wallingford, 1977). All of 
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these compounds and processes are essential for plant 
metabolism which ultimately determines growth, development, 
and crop yield. 

The phosphate content of plant material is controlled 
by two factors, the specific, genetically-fixed nutrient 
uptake potential of plants for phosphorus and the 
availability of phosphorus in the soil (Mengel and Kirkby, 
1978). The ability of the soil to supply phosphorus to 
plants can be separated into several general factors. 
Omanwar (1970) defined these factors as (i) intensity, the 
properties of the soil phosphorus that affects the ease or 
difficulty of phosphorus withdrawal by plants, (ii) 
quantity, the total amount of the nutrient reserve in the 
soil that is available’ to the plant, and (iii) rate, the 
transport of phosphorus to roots. Numerous researchers have 
attempted to relate these factors, either individually or in 
combination with each other, to crop growth. Studies have 
shown that these factors are not independent of each other, 
nor are they independent of the chemical and physical 
properties of the soil. Recognizing this, various techniques 
for evaluating and modelling the influence of soil 
properties on crop response to fertilizer have been used. 


Therefore, the objectives of this literature review are: 


1. to review the concepts of phosphate 
equilibria in soil and the influence of soil 


properties on the reactions of added 
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phosphate fertilizer. 

to examine some of the techniques that have 
been used to evaluate the phosphorus 
fertility status of soils. 

to review the mathematical models used for 


characterizing crop response to fertilizer. 





B. SOIL PHOSPHORUS EQUILIBRIA AND REACTIONS OF ADDED 
PHOSPHATES WITH SOIL 

The immediate source of phosphorus for the plant is the 
soil solution, but the phosphate concentration in this 
solution is very low, in the order of 1 to 0.1 ug ml-'. 
Within the soil solution, the forms of phosphorus are in 
equilibrium governed by protonation reactions and ionic 
complex formations. The ionic species of phosphates that are 
commonly found in the soil solution include H;P0,, H,P0O,°, 
HPO, 7, and PO,~7, with the most abundant being H,PO,~™ and 
HPO,~? (Larsen, 1967). Also, many metallic ions form soluble 
complexes of varying stability with phosphorus but following 
Bremgeneral) Order (Of) Fes s>GAle *soiMne? > Cat*,Mgcte :>wK*,;Na* 
(Sillen and Martell, 1964). In addition to this solution 
equilibrium, solution phosphorus is also in equilibrium with 
the phosphorus in the solid phase, but such that it heavily 
favours the solid phase. Hence, it is this latter overall 
equilibrium which controls the phosphorus concentration in 
solution. As plant roots remove phosphorus from the soil 
solution, phosphorus from the solid phase enters the 
solution phase in an attempt to reestablish the overall 
equilibrium. The rate of phosphorus dissolution from the 
solid phase is controlled by the forms of phosphorus in the 
solid phase which in turn are a function of the physical and 
chemical properties of the soil. 

The influence of the soil environment on the intensity 


factor, and the chemical properties of orthophosphates in 
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the soil, particularly fertilizer phosphates has been well 
studied (Swenson et a]., 1948; Dean, 1949; Wild, 1950; — 
Kurtz, 1953; Olsen, 1953; Hemwall, 1957; Mattingly and 
Talibudeen, 1967; Williams, 1970; Soper and Racz, 1980). 
Phosphates added to soils react strongly with various soil 
components with the most commonly suggested reaction 
mechanisms being physical and chemical sorption, anion 
exchange, surface precipitation, and precipitation as 
separate solid phases. In a review of phosphorus fixation by 
soils, Hemwell (1957) indicated that the recovery of 
fertilizer phosphorus by crops amounts to only 10 to 30 % of 
the quantity added to the soil with the remaining 70 to 90 % 
being primarily chemically precipitated and physiochemically 
sorbed by the soil. The soil properties and components that 
play important roles in these reaction mechanisms include 
pH, aluminum and iron hydrous oxides, alumino-silicate 
minerals, carbonates, non-living organic matter, moisture, 
and the ionic nature of the soil solution. 

Time iS an important aspect of the reactions of added 
phosphates with soil and can be separated into initial and 
long term categories. Soper and Racz (1980) describe the 
dissolution of fertilizer phosphate granules in moist soils 
as being fairly rapid, forming a saturated phosphate 
solution around the granules. As this phosphate rich 
solution moves into the Surrounding soil, alteration of soil 
constituents and solution composition occurs, resulting in 


precipitation and adsorption reactions. The initial reaction 
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products are metastable and are altered to more stable and 
less water soluble products over time, with the rate of 
alteration being controlled by soil properties and 
environmental factors. Thomas and Peaslee (1973) found, from 
fractionation studies, that added phosphates will assume the 
pattern of native phosphates with time and that over a 
number of years, the various fractions merely build up, 
about in proportion to their original content. 

The most important soil property which appears to 
control phosphate behavior in soils, in terms of ionic 
Species, initial chemical reactions and final products, is 
soil reaction (pH). The prevailing soil pH has a definite 
relationship with some predominant reaction mechanisms of 
phosphate retention by soils. As indicated earlier, one of 
the reactions controlling the species of phosphate ions in 
solution is protonation. As a result, H,PO,> tends to be the 
dominant phosphate species under acidic conditions, while 
HPO,~? is dominant under alkaline conditions. It is 
difficult to separate the direct effects of pH on phosphate 
behavior from those of other soil properties such as 
mineralogy and exchangeable cations. Extensive reviews by 
Dean (1949), Wild (1950), Hemwall (1957), Smith (1965), 
Larsen (1967), Ryden et al]. (1973), Parfitt (1978), and 
Soper and Racz (1980) have dealt with the mechanisms of 
phosphate retention and "fixation" by soils. 

In soils with acidic pH, the reactions of added 


phosphates are dominated by Al, Fe, and, to some extent, Mn 





to produce basic reaction products. Al and Fe sesquioxides, 
which can occur as discrete compounds or as coatings on soil 
particles, have been implicated as playing a significant 
role in phosphate retention. Depending on time, temperature, 
pH and phosphate concentration in the soil solution, these 
compounds can retain large quantities of added phosphate 
(Wild, 1950). The suggested mechanism by which Fe and Al 
oxides retain phosphates has been separated into three 
stages of adsorption occurring at different solution 
phosphate concentrations: (i) a high energy chemisorption of 
small amounts of phosphate; (ii) precipitation of a separate 
phosphate phase; and (iii) a low energy sorption of 
‘phosphates onto the precipitate (Bache, 1964). Hingston et 
al., (1967,1968) have suggested and shown a specific 
adsorption mechanism for hydrous Fe and Al oxides by which 
the phosphate is capable of exchanging with edge OH, and OH™ 
groups and becoming coordinated with the Fe or Al ion at the 
surface. Phosphate adsorption has been correlated with 
either extractable Al or Fe (Lopez-Hernandez and Burnhan, 
1974; Evans and Smillie, 1976; Schwertmann and Knittel, 
1973: Myszka and Janowska, 1973) and with exchangeable Al 
and Fe (Udo and Uzu, 1972). Precipitation of phosphorus by 
Al and Fe is also considered to be significant by Ghani and 
Islam (1946) but Hsu (1964) and Fitter and Sutton (1975) 
found this only in soils having pH less than 5, due to the 
low activities of Fe** and Al** in soil solution at pH 


values above 5. 
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In alkaline soils, the reaction of phosphate fertilizer 
and the solubility of phosphates is influenced by Ca*? 
and/or Mg*? and CaCO,, resulting in more stable basic Ca and 
Mg phosphates being formed. Alkaline soils can be separated 
into calcareous and non-calcareous depending on the presence 
of carbonates. The reactions of phosphates in a 
non-calcareous soil would be dominated by exchangeable Ca’*? 
and/or Mg*?. Olsen (1953) demonstrated that if the soil pH 
is raised by additions of NaOH, the solubility of calcium 
phosphate increases, but if Ca(OH), is added to increase the 
pH, the solubility of the calcium phosphate decreases as a 
result of the commom-ion effect. He indicated that the 
mechanisms for this common-ion effect are precipitation 
reactions forming calcium phosphates, and adsorption 
reactions with calcium on clay minerals forming a monolayer. 
Larsen et a]. (1965) found that pH was significantly 
correlated with the half-life of the labile phosphorus 
measured. They suggested the decrease in labile phosphorus 
was due to the formation of a crystalline basic calcium 
phosphate at a rate that increased with pH. The reactions of 
phosphates in a calcareous soil are again dominated by Ca*? 
and/or Mg*? with CaCO, and/or CaMg(CO;), acting as a source 
of calcium and/or magnesium and also as a pH buffer (Soper 
and Racz, 1980). Bell and Black (1970) found the change of 
the initial reaction products to more basic compounds was 
more rapid when CaCO, was present. Added phosphates also 


react with the carbonate particles themselves by forming a 





surface coating on the particles. With time, the layer of 
phosphate on the carbonate particles may be coated by more 
carbonates (Thomas and Peaslee, 1973). Parfitt (1978) 
Suggested three steps involved in the reaction: 

(i) chemisorption of phosphate, accompanied by heterogeneous 
formation of nuclei of amorphous calcium phosphate; (ii) a 
Slow transformation of these nuclei into crystalline calcium 
phosphate; and (iii) crystal growth of calcium phosphate. 
The result is a tendency for the solubility of the 
"adsorbed" phosphate to decline with time and thus decrease 
the phosphate availability. 

The investigation of the retention of phosphates by 
alumino-silicate clay minerals have been extensive. Wild 
(1950) reported that silicate clays could sorb phosphorus by 
several mechanisms. These include an exchange reaction of 
phosphates with OH groups on an edge Al-OH (ligand exchange) 
and/or an anion exchange reaction at a positively charged 
Site developed by the adsorption of protons on -OH groups. 
Dissolution of clay minerals to release Si and subsequent 
precipitation of phosphorus as alumino-phosphate compounds 
has also been proposed, but only at high phosphorus 
concentrations (Low and Black, 1948,1950; Rajan and Fox, 
1975). The rate of phosphate retention by clay minerals 
generally increases with temperature, concentration of 
phosphorus, and decreasing pH, and follows a decreasing 
order of illite, kaolinite, and montmorillonite (Haseman, 


1950). Exchangeable cations also influence the retention 





capacity of clay minerals. As already stated, exchangeable 
Fe and Al have been correlated with phosphate adsorption 
under acidic conditions. Kurtz (1953) pointed out that Ca 
clays retain more phosphates than do Na, NH,, or K clays. It 
is possible that the linkage of phosphates to the clay 
particle may be through exchangeable Ca*? or Mg*? ions 
acting as a bridge. Blanchet (1974) illustrated the 
influence of physio-chemical properties of the soil 
(particularly particle size) on plant nutrition. He compared 
the amount of phosphate absorbed/gram of root with 
increasing phosphate additons for two soils, a sandy loam 
and a clay. The amount absorbed was greater for the sandy 
loam than the clay due to the higher adsorption properties 
of the clay. 

The influence of organic matter on the retention of 
phosphates in soil has been studied by many workers. Soil 
Organic matter, and more specifically, humus, is considered 
to have very little ability to retain phosphates due to its 
normal negative charge. However, because of this negative 
charge, it can hold many cations which can react with the 
phosphate ion. Doughty (1930, 1935) gave evidence that Fe*?, 
Al*?, and Ca*? ions which are associated with the organic 
matter can react with phosphates. Several researchers have 
reported positive relationships between organic matter 
content of soils and phosphate adsorption (Rennie and 
McKercher, 1959; Harter, 1969; Hinga 1973; Lopez-Hernandez 


and Burnhan, 1974; Holford and Mattingly, 1975). By 





contrast, Moreno et al]., (1960) demonstrated that organic 
matter may complex Ca ions and thus increase the phosphate 
concentration in the soil solution from some of the calcium 
phosphates. Replacement of phosphate ions adsorbed by clay 
minerals by the humate ion has been shown by Mattson (1931). 
Nagarajah et a]., (1970) found that organic acids were 
capable of reducing the amount of phosphate adsorbed by 
kaolinite, gibbsite, and goethite by what they believed to 
be a ligand exchange mechanism on the mineral surfaces and 
thus the organic acids compete with phosphates for 
adsorption sites. Phosphate and organic matter competition 
has also been suggested for adsorption on CaCO, surfaces in 
calcareous soils (Holford and Mattingly, 1975). Thus, the 
evidence suggests that organic matter may either decrease or 
increase the ability of soils to adsorb phosphorus. 

Soil moisture influences phosphorus nutrition of plants 
by affecting many soil factors and processes which control 
the supply of phosphorus to the plant. These include 
transport rates, adsorption-desorption rates, mineral and 
precipitate solubility, and mineralization and 
immobilization rates. As the moisture content of soil 
decreases, adsorption-desorption equilibria would favour 
adsorption, the solubility of phosphate minerals and 
precipitates would decrease, and biological activity would 
decrease, reducing mineralization of organic phosphorus 
(Sheppard and Racz, 1980). Olsen et a/]., (1961) concluded 


that reducing the soil moisture content reduced phosphorus 
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uptake because (i) it reduces the movement of phosphorus to 
the root by reducing the thickness of water films which 
increases the diffusion path length, and (ii) it reduces the 
amount of phosphorus absorption by the root by reducing the 
number of root hairs, elongation of roots and turgidity of 
roots. Simpson (1965), Reichman and Grunes (1966), and 
Strong and Barry (1980) found that the availability of 
native phosphorus is more sensitive to soil water content 
than the availability of fertilizer phosphorus. In additon, 
Strong and Barry (1980) suggested that the reduced 
utilization of native soil phosphorus under dry conditions 
was the result of the reduced soil volume exploited by the 
stunted root system. AS a consequence of this and the 
relatively high availability of phosphorus in the fertilizer 
band, there may be a relatively large crop response to 
phosphorus under arid conditions. 

The presence of soluble salts in association with 
phosphate fertilizer materials influences phosphate 
availability. The common-ion effects of Ca salts have 
already been cited as decreasing phosphate availability. An 
increase in phosphate availability may be accounted for by 
an increased stimulation of the plant due to the presence of 
the salts or by chemical effects on the phosphate reaction 
products in the soil. Addition of nitrogen to a phosphate 
fertlizer band has been reported by many workers to increase 
the phosphate absorption by the plant. This has been 


attributed to (i) increased root growth in the vicinity of 
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the band (Duncan and Ohlrogge, 1956; Grunes et a/]., 1958; 
Miller and Ohlrogge, 1958), (ii) increased solubility of the 
phosphate fertilizer (Bouldin and Sample, 1958, 1959; 
Starostka and Hill, 1955), (iii) increased metabolic 
activity of the plant (Cole et al]., 1963; Leonce and Miller, 
1966; Minshall, 1964), and (iv) a reduction in pH at the 
soil-root interface, most likely caused by the exchange of 
Heeiens *from within«the .root for NH;* or -K*-ions in soil 
(Miller et a]., 1970; Riley and Barber, 1971). Bouldin and 
Sample (1958), studying the influence of associated salts on 
plant availability of concentrated superphosphates, found 
the order of effectiveness to be generally KNO, > (NH,).,S0, 
> NH,NO,;, > NH,Cl > KCl. Whatever the mechanism, the 
literature does indicate a definite increase in phosphate 
absorption by plants when nitrogen is in close contact with 
the phosphate fertilizer. Several workers (Mitchell, 1957; 
Olsen et a]., 1954) have demonstrated an appreciable 
increase in plant availability of rock phosphate and other 
Paesphate Carriers from the use of sulphur, while no 
appreciable influence of potassium on phosphorus uptake 
could be demonstrated (Olsen et a] ., 1954; Fine, 1955). 

The interaction of phosphorus with other elements in 
the soil may influence the crop response to phosphate 
fertilizer and the availability or utilization of many other 
elements. Nitrogen effects have already been cited, but in 
addition micronutrient-phosphorus interactions have been 


studied, as reviewed by Adams(1980). Micronutrient 
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deficiencies, induced by phosphate application, have been 
noted. Racz and Haluschuk (1970) reported the effects of 
phosphorus levels on the utilization of Cu, Zn, Fe, and Mn 
by wheat and flax on Manitoba soils. They found that trace 
element content and uptake by these crops were reduced in 
many instances when large amounts of phosphorus were added 
to soils or nutrient solutions. They concluded that the 
reduction in trace element uptake was due to the inability 
of the plant, under high phosphorus levels, to absorb the 
trace elements. For soils having amounts of available 
micronutrients which could be considered as bordering on 
deficiency, addition of phosphate fertilizer could induce 
‘Micronutrient deficiencies. In order to achieve maximum 
plant growth, both macro and micro nutrients must not be 
limiting. Leibig proposed in his Law of the Minimum that the 
amount of plant growth was controlled by the factor present 
in the minimum amount, and implied that if two factors are 
limiting, or nearly limiting growth, adding only one of them 
will have little effect on growth, while adding both 
together could have considerable effect (Russell, 1961). 
Therefore, if a soil is deficient in both phosphorus and a 
micronutrient, addition of phosphorus alone could result in 
a small degree of crop response or have no effect on crop 


growth. 





C. EVALUATION OF PHOSPHORUS FERTILITY STATUS OF SOILS 

The evaluation of the phosphorus fertility status or 
quantity factor of soils has been extensively studied (Olsen 
et al]l., 1954; Miller and Axley, 1956; Robertson, 1962; 
Omanwar, 1970; Alexander, 1973; Gwyer, 1979). Omanwar (1970) 
stated that the use of the term "available" requires that 
some time limit be specified since all soil phosphorus could 
be mobilized and made available to plants over an infinite 
time period. In general, the term has been associated with 
one crop growth period, and implies that prior to crop 
growth, the soil has a particular level of phosphorus 
reserve which could be made available to plants during the 
growing season. 

Various methods for determining the phosphorus 
fertility status of soils have been developed and used. 
These include anion exchange resins, radioisotope 
techniques, and equilibrium isotherm techniques, but the 
most common method is chemical extraction by one of a 
variety of solutions including water, acidic solutions, 
alkaline solutions, and neutral salt solutions. The original 
approach to the problem was to attempt a dissolution of the 
Same amount of phosphate from the soil as would the plant 
roots (Russell, 1961). This concept was soon abandoned and 
the present approach involves selection of a method for 
which there is a high correlation between extractable soil 
phosphorus and phosphate uptake, yield, or yield response to 


phosphate fertilizer. Kamprath and Watson (1980) described 





the objectives of the phosphorus soil tests as being 

(i) grouping of soils into classes for the purpose of making 
phosphate fertilizer recommendations, (ii) prediction of the 
probability of getting a profitable response to application 
of phosphate fertilizer, and (iii) providing an index of the 
amount of phosphorus a soil can supply. These objectives can 
be restated as (i) separation of soils as responsive or 
unresponsive to phosphate fertilizer, (ii) prediction of an 
expected yield response to phosphate fertilizer, and 

(iii) prediction of the phosphate ferilizer rate that needs 
to be applied to attain an optimum yield. 

The two chemical extraction methods used presently in 
western Canada are a modification of the acid fluoride 
solution used by Miller and Axley (1956), and a buffered 
sodium bicarbonate solution developed by Olsen et al., 
(1954). The Miller and Axley procedure uses a 0.03N NH,F in 
0.03N H,SO, solution. The hydrogen of the H,S0O, greatly 
increases the solubility of all calcium phosphates; in 
addition it attacks aluminum and iron phosphates, although, 
the rate of dissolution of the aluminum and iron phosphates 
is somewhat slower than the calcium phosphates. Generally, 
it has been observed that the H* remove phosphates in the 
order Ca > Al > Fe. The SO,°>? forms weak complexes with 
polyvalent metal cations but competes poorly with phosphates 
for iron and aluminum. Sulphate appears to prevent 
readsorption of phosphate released by hydrogen ions. 


Fluoride ions specifically precipitate soluble calcium as 
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CaF, and as a result will extract the more soluble calcium 
phosphates such as CaHPO, from the soil. Fluoride also 
complexes aluminum stongly and frees phosphates bonded to 
aluminum. The fluoride ion is rather harmless to basic 
calcium and iron phosphates unless the fluoride solution is 
acidified (Thomas and Peaslee, 1973). The Miller and Axley 
procedure is considered most Suitable on neutral to slightly 
acidic soils (Olsen and Dean, 1965). Difficulties may arise 
when it is used on calcareous soils because of 
neutralization reaction between carbonates and the acid, 
resulting in low values for extractable phosphorus (Olsen et 
al., 1954). 

The Olsen procedure uses a 0.5M NaHCO, solution 
buffered at pH 8.5. The presence of HCO;” decreases the 
activity of Ca*? by causing precipitation of calcium as 
CaCO;. This results in increased solubility of calcium 
phosphates which are thought to be a major source of plant 
available phosphorus in calcareous soils. In addition, 
bicarbonate ions remove aluminum bound phosphates, probably 
by replacement and by aluminum precipitation because of the 
OH- ion content in the solution (Thomas and Peaslee, 1973). 
Extractable phosphorus by the Olsen method is usually better 
correlated with plant response on calcareous soils than is 
extractable phosphorus by acidic extraction methods (Olsen 
et aj., 1954). This is thought to be a result of the 
buffered nature of the extracting solution making it more 


Suitable for extracting calcium phosphates. 
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The amount of phosphorus extracted by both methods has 
been found to be highly correlated with "A" value 
measurements of plant available phosphorus (Olsen et a]/., 
1954: Omanwar, 1970; Omanwar and Robertson, 1973). As well, 
extractable phosphorus by these methods has been shown to be 
highly correlated with yield response. Robertson (1962) ina 
greenhouse study of 79 Alberta soils, found that the 
response of barley was highly correlated with extractable 
phosphorus as measured by both methods. Correlations ranging 
from R = 0.73** to R = 0.79** for the Miller and Axley 
method and correlations of R = 0.73** to R = 0.82** for the 
Olsen method were found. Numerous other studies in the 
greenhouse have shown high correlations between phosphorus 
extracted by these methods and plant response (Olsen et a/., 
1954; Maclean et a]., 1955; Miller and Axley, 1956; Martar 
and Samman, 1975). Holford (1980) compared several phosphate 
extraction procedures to determine the effects of phosphate 
buffer capacity of a soil under field conditons. The 
phosphate buffer capacity is the resistance of the soil 
solution concentration to change when phosphate is added to 
or removed from the labile pool (Holford and Mattingly, 
1976). Holford's results confirmed that the larger the 
negative effect of buffer capacity on extraction of labile 
phosphate by a soil test, the higher was the correlation 
between the soil test and plant response to phosphate. He 
found that the Bray (ammonium fluoride) method was the most 


sensitive to the buffer capacity of a soil while the sodium 


21 


bicarbonate extraction was less sensitive. Whereas a 
previous study suggested that the ammonium fluoride test was 
over-SsensSitive to buffering, and hence underestimated 
available phosphate in strongly buffered soils, this field 
study showed that the test was correctly sensitive to 
buffering. Consequently critical levels for near-maximum 
wheat yields do not vary for the ammonium fluoride test, but 
increase with the increasing buffer capacity for the sodium 
bicarbonate tests. As a result, an additional measurement of 
buffer capacity is therefore required to give precision in 


the use of the sodium bicarbonate soil test. 
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D. MATHEMATICAL MODELS IN SOIL RESEARCH 

Mathematical models are quantitative techniques for 
expressing the relationship between two or more variables. 
Numerous statistical procedures have been developed to 
evaluate, explain, and model experimental results, ranging 
from purely graphical to multiple regression and 
multivariate procedures. Probably the first method was by 
Simple observation of the data and arbitrary separation. In 
an attempt to separate responsive and unresponsive sites to 
fertilizer application, Cate and Nelson (1965) developed a 
graphical method for partitioning a scatter of percentage 
yield versus soil test level into two groups (i) those for 
which probability of response to added fertilizer is large 
and (ii) those for which probability of response to added 
fertilizer is small. They attempted to find the "critical 
level" soil test value for separating the two groups. In 
1971, these same authors outlined a statistical procedure 
for partitioning soil test correlation data into two classes 
of probable response to fertilizer (low and high), based 
upon maximization of the class sum of squares in a one-way 
analysis of variance. This sum of squares reflects the 
weighted sum of squares of the difference between the 
percentage yield means for the various classes and the grand 
mean. Using this procedure, one finds quantitatively the 
best divisions from the point of view of maximizing mean 
differences among classes. These, in turn, should be the 


best divisions from the point of view of prediction. The use 
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of more elaborate techniques of data analysis have become 
increasingly common due to the recognization of the 
influence of many site properties on the results of field 
experiments or observations. These techniques include yield 
response functions, multiple regression analyses, 
Simultaneous equations, diScriminant function analyses, and 
principle component analyses. Many of these techniques have 
been used in soil fertility studies, while others show great 
promise. 

As indicated earlier, crop yield is a function of many 
growth, or input factors. Dillon (1977) simplified this 
Situation by uSing a theory of response based on the 
important input factors. His theory contained three 
simplifying assumptions, 

1. there 1S a continuous smooth causal relation 

between the X's (inputs) and Y (outputs); 
2. diminishing returns prevail with respect to 
each input factor, X, so that the additional 
output from succeeding units of X becomes 
less and less; 

3. decreasing returns to scale prevail so that 
an equal proportionate increase in all inputs 
results in a less than proportionate increase 
in output. 

Crop response to successive fertilizer nutrient increments, 


a single input variable, follows these assumptions. 





24 


Some researchers have attempted to develop models for 
the effect of nutrient application on crop yield ona 
theoretical basis so that biological and physical meanings 
can be attached to their parameters. However, care is needed 
Since such models could be used to express a particular 
bias. Alternatively, models may be chosen for their 
computational convenience, the statisical estimation of 
functions from data or to permit calculation of optimal 
rates. As yet, there is no fundamental theoretical model for 
the effects of nutrient application on crop yield, but 
rather the model chosen is empirical, based on observations 
and experience (Colwell, 1978). In general, the mathematical 
expressions that have been used to relate crop growth to 
nutrient levels in the soil fall into three categories, 
namely; the straight line, the exponential, or the quadratic 
equation. Response functions for a single nutrient, such as 
phosphorus, have generally been either exponential or 
quadratic expressions. Characteristically, the exponential 
function never reaches a maximum and will never indicate a 
yield decrease. By contrast, the quadratic function does 
reach a maximum yield, followed by a yield decrease which 
could be due to a toxicity level of the factor, induced 
nutrient deficiency or a depletion of soil water by 
excesSive early vegetative growth stimulated by high 
fertilizer applications (Melsted and Peck, 1977; Colwell, 
1978). As a result, polynomial (quadratic) models are more 


popular than exponential models. Polynomial functions are 
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easily fitted to data by standard multiple regression 
procedures and can be made flexible enough to describe most 
trends and rigid enough to smooth out most errors in data 
(Colwell, 1978). Johnson (1953) compared quadratic functions 
with exponential functions for fitting response curves and 
concluded that the quadratic polynomial model generally gave 
the better fit and the best results for purposes of 
interpolation. Anderson and Nelson (1975) however, concluded 
that the use of second order polynomial models may result in 
biases in the estimates of optimal fertilizer rates due to a 
ceiling on the crop yield imposed by environmental or 
management factors and type of crop. 

Multiple regression analysis is an attempt to account 
for the variation in the dependent variable by a linear 
combination of independent variables. As mentioned in the 
previous discussion on yield response functions, multiple 
regression procedures are commonly used to fit equations to 
response data, but a more frequent use of multiple 
regression analysis in soil research has been the combining 
of a series of similar experiments. Studies using multiple 
regression analysis deal primarly with crop yield or 
compositon as influenced by fertility status of soils, 
fertilizer application, soil chemistry, site topography, and 
climate (Laird and Cady, 1969; Cady and Allen, 1972; 
Williams et al]., 1975; Bole and Pittman, 1980a,1980b). 
Agronomic experiments on the same factor or a group of 


factors are usually repeated for a number of years at one or 
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more locations. Because of the variation in the effect of 
many factors due to location and year, the results obtained 
at a single site for a single year are not precise enough as 
a basis for generalization (Leonard, 1962). Yates and 
Cochran (1938) stated that it is impossible to lay down 
rules of procedure for combining several experiments for 
different years which will be applicable in all cases, and 
that the results usually require comprehensive examination 
with special emphasis on certain treatment effects. 

The mathematical analysis of a complex problem can lead 
logically to a system of simultaneous equations (Heapy, 
1971). If the model can be divided into specfic stages such 
that a multi-equation system can be used to describe the 
model and where such models involve jointly determined 
variables, simultaneous equations procedure should be used 
(Dillon, 1977). Heapy et a]., (1976a,1976b) used this system 
of multi-equations to develop a barley yield equation based 
on the effects of soil and fertilizer nitrogen and 
phosphorus. As part of this equation, a moisture stress term 
was included but calculated from a second equation derived 
from data external to the study. 

A special type of statistical analysis that has been 
used in soil research, as well as geology and biology, to 
classify an individual into one of two or more groups is 
discriminant function analysis. The objective of this 
procedure is to find a linear combination of the variables 


that maximally discriminate among groups. The technique was 
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first used by Fisher (1936) as a solution to a taxonomic 
problem and has since found limited application in soil 
Science. Cox and Martin (1937) used the technique to predict 
the presence of Azotobacter on the basis of pH, available 
phosphorus, and total nitrogen content of the soil sample. 
Most of the recent applications of discriminant analysis in 
soil science has been in soil genesis and soil 
classification (Oertel, 1961; Norris and Loveday, 1971; 
Bracewell and Robertson, 1973; Berg, 1980; Henderson and 
Ragg, 1980). 

Thus far, all the statistical techniques discussed have 
been of a single criterion and multiple predictor 
association, with the exception of the discriminant 
analysis, of which only the two group situation fits this 
association, but another type of analysis which has been 
proposed and used in soil fertility studies is an analysis 
of variable interdependence, principal component analysis 
(PCA). Ferrari (1965) illustrated the use of a system of 
Simultaneous equations for modelling the magnesium content 
of herbage and suggested the use of factor anaylsis or PCA 
to obtain these equations. Kyuma and Kawaguchi (1973), and 
Kyuma (1973a,1973b) illustrated the use of PCA as a method 
of fertility evaluation and grading for paddy soils..Using 
the new components formed by PCA, they were able to develop 
a multiple regression equation to account for 57% of the 
yield variation. Principal component analysis has some 


advantages over multiple regression. Interpretation of 
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multiple regression analysis is dependent upon the 
assumption that explanatory variables in the analysis are 
not strongly interrelated (Chatterjee and Price, 1977). 
However, the real world, particularly soils, does not behave 
in this fashion. Even when subjected to various analytical 
chemical procedures, the analytical results will be 
influenced by various soil and environmental properties, due 
to the empirical nature of some procedures. Therefore, use 
of PCA has potential in identifying and evaluating the 


interrelationships among soil properties. 
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III. MATERIALS AND METHODS 


A. BACKGROUND OF FIELD DATA 
1. Cooperators and Site Design 

In 1971, the Risk Adjusted Yield Potential (RAYP) 
project was initiated in Alberta to collect data for the 
purpose of improving prediction of fertilizer requirements 
based on soil tests. It was a joint endeavor involving the 
University of Alberta, Alberta Agriculture, Western 
Co-operative Fertilizer Ltd., Sherritt Gordon Mines Ltd., 
and the Agriculture Canada research stations at Beaverlodge, 
Lacombe, and Lethbridge. Field plots were set out throughout 
the province to study both nitrogen and phosphorus 
fertilizer requirements for a number of years varying with 
cooperator. In most cases, a one-factor-at-a-time 
experimental design was used for both nutrients.' The 
exception was Lethbridge Research Station, which used a 
factorial design. In most cases, there were two test crops, 
barley and rapeseed. 

In addition, data for wheat response to phosphate 
fertilizer in east-central Alberta were included in the 
present study.’ These latter field experiments were 
conducted over the same years as those of the RAYP project 


using a similar one-factor-at-a-time experimental design. 


‘Personal communication with Dr. M. Nyborg. 
2Personal communication with Dr. J. A. Robertson. 
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2. Location of Test Sites 

There were two main objectives of the RAYP project. The 
first was to find the potential yielding ability and the 
fertilizer needs of different textural classes of soils 
within each soil zone. The second objective was to compare 
crop response to fertilizer on stubble and fallowed fields.? 
With these objectives in mind, plot sites were selected by 
each individual cooperator. Legal location of plot sites 
used in this study and their cropping history are provided 
in Appendix A. 

3. Seeding, Fertilizer Application, and Harvesting 

Whenever possible, one crop, namely barley, was common 
to all experiments in the RAYP project. Galt barley was the 
most common variety, but some sites. were Seeded to Betzes or 
Conquest. Where rapeseed was used, Span was the most common 
variety but some sites were seeded to Echo or Torch. In 
general, both barley and rapeseed were grown at a site, but 
some sites had only one test crop. The wheat sites in 
east-central Alberta utilized Thatcher wheat. Crop varieties 
for each site are provided in Appendix A. 

For the phosphorus block of each of the RAYP sites, 
blanket applications of nitrogen, potassium, and sulphur as 
NH,NO;, K,S0., KCl, and Na,SO, were applied either with the 
seed or side banded. The wheat sites had only a blanket 
application of nitrogen as NH,NO,. The phosphate fertilizer, 
in the form of treble superphosphate and/or ammonium 


3Personal Communication with Dr. M. Nyborg. 
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phosphate, was generally placed with the seed and/or banded. 
The number of treatments and phosphate fertilizer rates 
varied among cooperators, ranging from five to nine 
including a check, and the number of replicates also varied. 
Plots were harvested at maturity, air dried, threshed and 
the grain yield recorded. Yield means for each phosphate 
treatment at each site used in this study are provided in 
Appendix F. 
4. Soil Sampling 

Soil cores were generally taken on a site basis prior 
to seeding and divided into subdepths of 0-15 cm, 15-30 cm, 
30-60 cm, and 60-90 cm. Samples were air dried and several 
analyses were conducted by the Alberta Agricultural Soil and 
Feed Testing Laboratory (ASFTL). At the initiation of the 
author's study, 1977, sites were revisited to collect 
additional surface samples for physical analysis and 
determination of organic matter and carbonate content. 
5. Growing Season Precipitation 

The precipitation during the growing season was 
recorded for most sites, but some sites lacked these data. 
For those sites lacking data, approximate values were 
estimated using meteorological data published by Alberta 
Environment and records of neighbouring sites. Precipitation 
values for each site used in the study are presented in 


Appendix C. 


we =~. 3836 ons ees myecuig: 


a 
i 


; i wars: ae 
—— 
~~ bari rere 
7 . 


eve) 





ge 


B. ANALYTICAL METHODS 
1. Soil Physical Analysis 

One composite sample of the 0-15 cm depth of each site 
was ground to 2 mm using a flail grinder. Particle-size 
analysis was performed on these samples by the hydrometer 
procedure (Bouycous, 1951; Toogood and Peters, 1953). 
Particle-size distribution and textural classification for 
each site used in the study are presented in Appendix D. 

2. Soil Chemical Analysis 

Most of the chemical analyses were done by the Alberta 
Agricultural Soil and Feed Testing Laboratory using 
composite site samples for each soil depth. Caco, 
equivalence and organic matter content were determined at 
the University of Alberta on either original samples or 
subsequent samples. 

Soil reaction (pH) was determined on a 1:2 soil:water 
Suspension using a pH meter. Electrical conductivity (E.C.) 
was determined on the same 1:2 soil:water suspension uSing a 
conductivity meter. The conductivity reading was multiplied 
by a factor of 2.063 to express results on a Saturated paste 
extract equivalent. 

Nitrate-nitrogen was determined on a 0.02N CuSO, plus 
0.007N AgNO, extract using the phenoldisulphonic acid method 
(Prince, 1945) as described by Heapy (1971). Extractable 
potassium was determined from a 1:5 soil:ammonium acetate 
extract using a flame photometer. Extractable phosphorus was 


determined using three procedures, the Miller and Axley, and 
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Olsen methods, as described by Alexander et a]. (1972), and 
a modified Miller and Axley procedure. This modified 
procedure utilized a 5 cc (scoop) volume soil sample, 25 ml 
of the 0.03N NH,F in 0.03N H,SO, extracting solution anda 
Shaking period of 10 minutes. After filtration, phosphorus 
in solution was determined on a auto-colorimeter set at 400 
nm using a combined nitric vanadate molybdate procedure 
(Kitson and Mellon, 1947). 

CaCO; equivalence was determined on the 0-15 cm samples 
using the calcimeter method (Bascomb, 1961). Organic matter 
content was obtained for the 0-15 cm sample by measuring 
total carbon by dry combustion using a Leco induction 
furnace, subtracting the portion that was inorganic carbon 
and multiplying by a factor of 1.71. Results of the soil 


chemical analyses are presented in Appendices C and E. 


C. SOIL CLASSIFICATION 

Many of the sites had been classified by some of the 
participants of the project. Those sites which were 
originally unclassified were revisited and classified 
according to the Canadian System of Classification (1978). 
Soil classification of each site used in the study are 


presented in Appendix B. 
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D. DATA ANALYSIS 
1. Response Functions 

For consistency in interpretation, only one type of 
mathematical expression was used for purposes of fitting a 
curve to the yield data. A second order polynomial equation 


of the form 


Yas be + bp4x) + ber ox 2) 


was calculated for each site using the mean yield of each 
treatment as the Y term and the phosphate fertilizer rates 
as the X terms. The coefficient values (b,, b,, and b,) for 
each site-equation are given in Appendix G. The effects of a 
nutrient application on yield are not immediately obvious 
from the yield functions. Therefore, a yield increase value 
was calculated for each site as follows: The maximum yield 
(Y-max) was calculated for each site by taking the first 
derivative of the equation, equating it to zero, solving for 
X (X-max, the fertilizer rate for Y-max) and inserting the 
value of X-max into the original equation to obtain the 
Y-max. Ninety percent of Y-max (Y-90%max) was selected as 
the "optimum yield" for each site since this value would 
likely be in the upper end of the "linear" portion of the 
quadratic curve, meaning that the fertilizer rate to obtain 
this yield should still be providing an economic return 
(Spencer and Glendinning, 1980). The rate of fertilizer 


(X-90%max) required for Y-90%max ("optimum" fertilizer rate) 
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was calculated from the original equation. Yield increase 
was calculated by the difference between Y-90%max and the 
"bo," value or yield at the zero phosphate fertilizer rate. 
Percent yield increase was calculated by dividing the yield 
increase by the Y-90%max and multiplying by 100. Percent 
yield increase was used to remove some of the variation in 
yield caused by environmental conditions and to take into 
account maximum yielding potential differences among sites 
(Colwell, 1978). The percent yield increase was used only 
for the wheat sites. The yield increase or percent yield 
increase was used to characterize the magnitude of the yield 
response to phosphate fertilizer and to provide a common 
yield term for each site for use in subsequent analyses. As 
a result of these procedures, sites could be separated into 
two groups, responsive and unresponsive. Responsive sites 
were those sites having a yield increase greater than 
0.0 100kg/ha (quintal/ha) while unresponsive sites had a 
yield increase equal to or less than 0.0 100kg/ha. Yield 
increase and all intermediate values are presented in 
Appendix H. 
2. Multiple Regression and Least Square Analyses of 

Covariance For Unequal Numbers 

The multiple regression function is a linear 

combination of independent variables that attempts to 
account for the variation in a dependent variable. The 


multiple regression equation is expressed in the form 
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where Y is the predicted value of the dependent variable, b, 
is the intercept of the regression line on the vertical 
axis; b,, b., ...- bn are the partial regression 
Poerlrcrvents tom Yon Pinte, Go) o. os oc kn fare vthe observed 
values of the independent variables. The sequence of 
addition of the independent variables in the regression 
equation is controlled by the proportion of the variation 
(R?) in the dependent variable that is accounted for by the 
independent variable. The greater the contribution to the 
overall R?, the greater the importance that variable has in 
accounting for the variation in the dependent variable. The 
level of significance of the regression equation and the 
individual b values of the independent variables are 
determined from the calculated F ratio. 

In multiple regression analysis, it is necessary to 
code qualitative variables with dummy values. An effect 
coding, which uses a 1, 0, -1 coded values, is one of the 
coding systems used to code qualitative variables. Although 
such systems of coding are valid for equal subclass numbers, 
they are most often used for unequal numbers.* The 
intercept, b., iS an estimate of the grand mean of the 
dependent variable, Y, and each b is an eStimate of the 
treatment effect for the group with which it is associated 
i.e., the deviation of the mean of the group from the grand 


*Personal communication with Dr. R. Hardin. 
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mean, Y. Subsequent to obtaining a significant R?, the mean 
Y value for each qualitative variable is determined by an 
analysis of covariance. The effects of the covariates are 
removed from the analysis and the qualitative variable means 
are adjusted to a common value of the covariates, usually 
the mean of the covariates. This type of covariance analysis 
requires the assumption that the slopes of the regression 
lines are equal among the qualitative variables. Significant 
differences between the qualitative variable means are 
determined by an approximate multiple comparison test. The 
product difference between two means used in computing, now 
accounts for the variance and covariance between the 
qualitative variables and the covariates (Harvey, 1975; 
Mehlenbacher, 1978: Steel and Torrie, 1980). 

In this study, stepwise multiple regression equations 
were computed for the "responsive" sites. Yield increase was 
used as the dependent term to determine the influence of 
soil properties on crop response to phosphate fertilizer 
using multiple regression procedures. Variables considered 
for inclusion were: three soil test procedures, a number of 
quantitative site variables, and qualitative site variables. 
The use of multiple regression is based on the assumption 
that the relationships between the dependent variable and 
the independent variables are linear. To determine whether 
or not this was in fact the case, the dependent variables 
(yield increase or percent yield increase) were plotted 


against each of the soil test phosphorus variables using a 


-_ 


=f 
Stiteetee 7 egw: 





38 


Scattergram program. Visual examination of the scattergrams 
indicated that a nonlinear relationship did exist. The 
natural logarithmic transformation was regarded as best 
approximating a linear relationship where originally a 
nonlinear relationship existed. The effectiveness of the 
transformation was evaluated by the contribution to the 
overall correlation between the dependent variable and the 
best set of independent variables before and after 
transformation. The contribution to the overall correlation 
between the dependent variable and the set of independent 
variables should be greater for transformed independent 
variables than for the non-transformed independent 
variables. The natural logarithmic-transformed variables 
were subsequently used as separate independent variables in 
evaluating classification variables. It should be noted that 
if the non-transformed variable was equal to 0, then O was 
used for the value of the natural logarithmic 
transformation. 

After establishing the best combination of quantitative 
variables, the classification dummy variables were inserted 
into the analysis using an effect coding (Appendix I).5 
Analysis of covariance for unequal numbers was used to 
calculate qualitative treatment means at the means of the 
quantitative covariates. Student-Newman-Keuls' test was used 
for approximate comparison of these means (Steel and Torrie, 


S5Personal communication with Dr. R. Hardin. 
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3. Discriminant Analysis 

In theory, the discriminant function is a linear 
combination of independent variables with a dependent 
variable that represents group membership. With only two 
groups, discriminant function analysis amounts to multiple 
regression analysis with the dependent variable taking the 
values of 1 and 0 (Kerlinger and Pedhazur, 1973). Stepwise 
discriminant analysis begins as a simple one-way analysis of 
variance, based on the highest F value of the variable that 
best discriminates between groups. A second discriminating 
variable is selected as the one best able to improve the 
power of discrimination in combination with the first 
variable. At each step, variables may be removed if they 
reduce discrimination when combined with more recently 
selected variables. Eventually, all variables which 
Significantly contibute to the discriminating power are 
included. (Klecka, 1975; Berg, 1979). 


The discriminant function is expressed in the form 


D = az 3e agz) =f ay 2a Os Pe as 2 dpZp 


for unstandardized data and in the form 
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for standardized data. D represents the score on the 


discriminant function, d, is the constant, the d's are the 
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weighting coefficients, and the Z's are the values of the p 
discriminanting variables used in the analysis. Ideally, the 
discriminant scores (D's) for the cases within a particular 
group will be fairly similar. At any rate, the function is 
formed in such a way as to maximize the separation of the 
groups. The sequential addition of the independent variables 
to the function is dependent upon their discriminating 
power. The greater its ability to separate the groups, the 
greater the importance that variable has in the function. 
The relative importance of each variable is determined from 
the standardized discriminant function coefficients. When 
the sign is ignored, each coefficient represents the 
relative contribution of its associated variable to that 
function. The sign merely denotes whether the variable in 
making a positive or negative contribution (Klecka, 1975). 
The effectiveness of a discriminant function can be 
judged by two meaSurements. The total discriminatory power 
(TDP) which is a measure of the total variability of the 
function attributable to group differences can be calculated 
using the function eigenvalue which is a meaSure of the 
total variance existing in the discriminatory variables 
(Appendix J). A further aid in judging the importance of a 
discriminant function is its associated canonical 
correlation, a measure of how closely the function and the 
"group variable" are related. The canonical correlation 
Squared can be interpreted as the proportion of variance in 


the discriminant function explained by the groups. The 
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Statistical significance of the function can be measured by 
the chi-square statistic. 

The resulting equation indicates to which group each 
member probably "belongs". Thus the function can be used for 
predictive purposes to determine the membership of an 
unknown into one of the groups based on itsS measurement of 
certain properties as defined by the function (Kerlinger and 
Pedhazur, 1973). Classification is achieved through the use 
of a series of classification functions, one for each group. 
Classification equations are derived from the pooled 
within-groups covariance matrix and the centroids for the 
discriminating variables. The resulting classification 
coefficients are to be multiplied by the raw variable 
values, summed together, and added onto a constant. The 


equation for one group would appear as 
Peaeeeigl, er leis Vatit.. eet erpVp tech. 


where Ci is the midoss mace wien scoresfiopaqroupsei athe! ci 7s 
are the classification coefficients with ci, being the 
constant, and the V's are the raw scores on the 
discriminating variables. There is always a separate 
equation for each group and a case would be classified into 
the group with highest score (Klecka, 1975). 

For this study, a stepwise discriminant analysis 
program was used to compare three soil phosphorus test 


procedures, select the optimal set of discriminating 





42 


variables and to compute discriminant functions to separate 
responsive and unresponsive sites. The responsive sites were 
coded as 1 and the unresponsive sites were coded as 0. The 
effect of the natural logarithmic transformation of the 
phosphorus soil tests was determined. The effectiveness of 
this transformation was evaluated by the contribution to the 
canonical correlation between group membership before and 
after transformation. The contribution to the overall 
correlation should be greater for the transformed variable 
than the non-transfromed variable. If this was the case, 
then the natural logarithmic transformed variable was 
subsequently used as a separate independent variable in 
evaluating classification variables. It should again be 
noted that if the non-transformed variable was equal to 0, 
then the natural logarithmic transformation was assigned a 
value of 0. The criterion used to select discriminating 
variables was to maximize the Mahalonobis distance between 
the two groups. The procedure is fairly straight forward for 
data composed of only measured variables, but data 
consisting of both measured and qualitative variables tend 
to be more troublesome. Krzanowski (1980) demonstrated a 
method of discriminant analysis for mixtures of categorical 
and continuous variables using a binary (1, 0) coding of the 
categorical variables. The overall error rate was reduced 
when compared to a weighted coding (0, 1, 2), but still 
remained high. Kendall (1975) stated that there appears as 


yet to be no satisfactory theory to deal with this 
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Situation. He proposed either a dummy coding system or 
alternatively, a separate discriminant function for each 
class of the qualitative classification. Both procedures 
were tried in the present study using the effect coding for 
the dummy system and a separate discriminant analysis for 
each qualitative class containing sufficient members. For 
the dummy system technique, the procedure used was to first 
find the optimum set of independent measured variables using 
the stepwise procedure based on maximizing the Mahalonobis 
distance. Then, uSing a direct computing program option 
which enters all independent variables into the analysis 
concurrently, the dummy variables were included with the 
optimum set of measured variables. This is to insure that 
all dummy variables were included in the discriminant 
function. The stepwise procedure was only used for 
determining a separate discriminant analysis for each class 
of the qualitative classification. 
4, Principal Component Analysis 

In theory, principal component analysis (PCA) or factor 
analysis is a statistical procedure used to interpret within 
the variance-covariance matrix of a multivariate data 
collection (Davies, 1973). Rummel (1967) described the 
working of factor analysis as taking numerous measurements 
and qualitative observations and resolving them into 
distinct patterns of occurrence. No particular assumption 
about the underlying structure of the variables is required. 


The process of principal component analysis can be separated 
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into two steps. First a correlation matrix of the variables 
involved is computed as a meaSure of association. The second 
step is the extraction from the correlation matrix of 
initial components as eigenvalues and eigenvectors such that 
the components are orthogonal or independent of each other 
(Kim, 1975). 

Principal component analysis transforms a given set of 
variables into a new set of composite variables that would 
account for more variance in the data as a whole than any 
other linear combination of variables. The second component 
is defined as the second best linear combination of 
variables, under the condition that the second component is 
orthogonal to the first, and therefore can be defined as the 
linear combination of variables that accounts for the most 
residual variance after the effect of the first component is 
removed from the data. Subsequent components are defined 
Similarly until all the variance in the data is exhausted 
(Kim, 1975). 

In this study, a principal component analysis program 
was used to determine the interrelationships present among 
the independent variables for the responsive sites. The 
variables were standardized such that each variable had a 
mean of zero and a unit variance to ensure a normal 
distribution. This allows one to compare the distribution of 
one variable to that of another when the two variables are 


expressed in different units of measurement (Davies, 1973). 
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5. Selection of Independent Variables 

The variables chosen for the discriminant, multiple 
regression and principal component analyses were those 
considered to have an influence on the availability of soil 
phosphorus and crop response to phosphate fertilizer. These 
included soil chemical and physical properties, and 
qualitative classification variables. A problem that did 
arise during the study was missing data for Olsen-P,and 
Miller and Axley-P soil tests for some sites as a result of 
loss of original samples. Also, because of multicollinearity 
problems, % clay, % silt, and % sand variables could not all 
be used at the same time for the disciminant and principal 
component analyses. As a result, only % sand and/or % clay 
was used. 

The qualitative variables used in this study included 
agro-climatic area, soil zone, and soil order 
classifications. Agro-climatic area classification was 
determined from the Agro-climatic Areas of Alberta map 
(Bowser, 1967) while soil zone classification of each site 
was based on the Soil Zones of Alberta map (Odynsky, 1962) 
and identification of the soil great group using the 
Canadian System of Soil Classification (1976). Soil order 
classification according to the Canadian System of Soil 
Classification (1976) was based on profile examination for 


each site. 





IV. RESULTS AND DISCUSSION 

This chapter is divided into four sections. The first 
three sections deal specifically with the results and 
discussion of the individual crops studied. Each crop 
section is further divided into the three statistical 
procedures used, discriminant analysis, multiple regression 
analysis, and principal component analysis. The final 
section deals with the potential sources of error within 


this study. 


A. Barley 

The results presented in this section represent the 
Statistical analyses of the yield response of barley to 
phosphate fertilizer for 125 sites. A brief summary of the 
chemical and physical characteristics of the field sites 
used are presented in Tables 1 and 2, while frequency 
distribution of the sites in regards to three types of site 
classification and textural classes are presented in 
Table 3. In general, the sites used in this study 
represented a wide variety of site conditions for both the 
responsive and unresponsive groups. The means and standard 
deviations of each independent variable were almost equal 
between the two groups. The most noteworthy difference 
between the two groups was a lower mean value of all three 
“soil test methods, ASFTL-P, M & A-P, and Olsen-P, for both 
depths of the responsive sites. The distribution of the 


responsive and unresponsive sites among the classification 
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Table 1. Mean, Standard Deviation, Maximum, and Minimum 
Values of the Independent Site Variables for 
the Unresponsive Barley Sites 


No. of 
Variables* Mean Std. Dev. Max. Min. Sites 
PH (0-15) 6.91 0.60 8.0 Shae 60 
PH (15-30) Onn 0.74 8.4 By 60 
Esc. ((0= 15) 0.41 0230 Pap) OreZ 60 
Bec. (075-30) 0.48 0.63 4.5 O82 60 
IOeM. £CO-115)) 5 ys2 2255 1:4 7, ee 60 
SEaCOs sO 1:5) Ora? ord 3.6 0.0 60 
% sand’ (0-15) 32223 18.85 78.6 353 60 
% silt (0-15) 36.51 He 8.7 58.0 4.0 60 
¥sclaget0=45) 31.24 16.67 Siva 8.2 60 
_ Pptn. 227529 Fiat 41.9 Tae 60 
ASFTL-P(0-15) = > 4950 393 218.4 455 60 
Ln ASFTL-P(0-15) 3.64 0.73 5.39 1.50 60 
ASFTL-P( 15-30) 225 43755 2 18R9 0.0 60 
Ln ASFTL-P( 15-30) 1-597 1.48 oi Si 0.0 60 
M & A-P(0-15) 5142 39a 22106 Sele 57 
Ln M & A-P(0-15) Sa 0.70 5.40 he. 7424 BZ 
M & A-P( 15-30) 26.1 43.4 ZO7e O0 57 
Ln M & A-P(15-30) ZU TALS Sass 0.0 57 
Olsen-P(0-15) 37 at iitsla si 12260) 1A #2 57 
Ln Olsen-P(0-15) Se 2 0.42 4.81 2ak2 57 
Olsen-P( 15-30) 22.6 TARO 89.6 Seyi 57 
Ln Olsen-P( 15-30) 2393 0..58 4.50 1.91 57 
* Variable Units 
Bec mmhos/cm? 
Pptn. cm 
ASFTL-P kg/ha 
M & A-P kg/ha 


Olsen-P kg/ha 
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Table 2. 


Mean, Standard Deviation, Maximum, and Minimum 


Values of the Independent Site Variables for 
the Responsive Barley Sites 


Variables* 


pH (0-15) 

pH (15-30) 

Bee (O=1/5)) 

E.C. (15-30) 
OM nC O= 1.5) 
Caco, (0=15) 

% sand (0-15) 
Fost le, (O=—15) 

% clay (0-15) 
Pptn. 
ASFTL-P(0-15) 

Ln ASFTL-P(0-15) 
ASFTL-P( 15-30) 

Ln ASFTL-P( 15-30) 
M & A-P(0-15) 

Ln M & A-P(0-15) 
M & A-P(15-30) 

Ln M & A-P( 15-30) 
Olsen-P(0-15) 

Ln Olsen-P(0-15) 
Olsen-P( 15-30) 

Ln Olsen-P(15-30) 


P,0, (90% Max.Yld.) 


* Variable 


eG. 
Pptn. 
ASFTL-P 
M & A-P 
Olsen-P 


P,0, (90% Max. 


Mean 


oOMN000) 07 


21 


26. 
Phe 
6. 
ihe 

203 
Sie 
8. 
ule 

26. 
3 

14, 
PAE 

7a 


276 
<83 
- 36 
AAS 
ath 
siaih 
An 
38. 
26. 
244 


58 
61 
80 


2 
98 
8 
43 
6 
06 
5 
78 
1 
03 
. 
56 
5 


Vatclien) 


Std. Dev. 


0.76 
i. O12 
OS 
Oe 17 
2.26 
2.14 
NBS 
10.85 
10.98 
Tien 
20.6 
0.82 
Wewae 
1.09 
eerenh 
Oren 2 
URS 
0.94 
1h AKe) 
Ooi 
Ue 
0.50 
55:0 


Units 


mmhos/cm? 
cm 

kg/ha 
kg/ha 
kg/ha 
kg/ha 
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Table 3. Frequency Distribution of Unresponsive 
and Responsive Barley Sites per 
Classification Class (Number of Sites) 


Classification Unresponsive Responsive 


Agro-climatic Area 


1 aie 21 
2A 8 8 
2H ifs: 19 
3H 1 8 
3Ha 16 9 
Soil Zone 
Gray 11 iS 
Dark Gray 14 17 
Black 11 20 
Thin Black 14 5 
Dark Brown 10 8 
Soil Order 
Chernozemic 41 38 
Luvisolic 17 24 
Gleysolic 0 2 
Solonetzic 2 1 


Texture (0-15) 
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classes waS approximately equal (Table 3) suggesting that 

Site classification may not be important in the separation 
of responsive sites from unresponsive sites. Finally, the 

general geographical distribution of the sites within the 

province was great enough to represent the major grain 


producing areas of the province. 
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1. Discriminant Analyses 

The objective of the use of discriminant analysis is to 
determine those site variables important for distinguishing 
phosphate responsive and unresponsive sites. Based on the 
Simple examination of the means in Tables 1 and 2, the most 
important variable for separation of the sites would appear 
be the soil test for phosphorus. Therefore, the first step 
was to determine the best soil test procedure for separating 
the sites. There appeared to be little difference among the 
three soil test procedures in their ability to separate 
responsive and unreponsive barley sites when the comparison 
waS made using the same sites (Table 4). The only major 
difference was that for the Olsen-P, both depths were 
important based on the standardized coefficients, while for 
both the ASFTL-P and M & A-P, the 15-30 cm depth did not 
enter the function. Even when all 125 sites were used in the 
analysis the ASFTL-P(15-30) did not enter into the function. 
There was little difference in the total disciminatory power 
(TDP) or the canonical correlation among the functions. 
Because there is little difference among the three methods, 
further discriminant analyses of the barley sites involved 
only the ASFTL-P because the number of sites having this 
information was larger than those sites having the Olsen-P 
OniM 6 -A-P.. 

The overall discriminant function based on the 
quantitative variables only (Table 4) for all 125 barley 


sites indicated that the most important site variable for 





Table 4. Discriminant Analysis for Barley Response to Phosphate 
Fertilizer: (1) Comparison of Soil Tests, and 
(2) Best Overall Function 


Std. 
Variables Coef. 


Unstd. 
Coef. 


1. Comparison of Soil 


(112 sites) 


Ln ASFTL-P(0-15) 1.00 
constant 


Ln M & A-P(0-15) 1.00 
constant 


Ln Olsen-P(0-15) 0.79 

Ln Olsen-P(15-30) 0.33 

constant ; 
(125 sites) 


Ln ASFTL-P(0-15) 1.00 
constant 


2. Best Overall 


Ln ASFTL-P(0-15) Onom 


* clay (0-15) 0.62 
DH G1S-30) 0.53 
Ppotn. 0.33 
Oe Mam Oma.) = Ove sul 
constant 


1.26 
Ste UE! 


1.41 
aaa Uke 


URIBE, 
0.61 
=6.96 


le2s 
=4..23 


Function (125 sites) 


24 
«05 
oe) 
04 
attr) 
a0 


wooo0o— 


** Significant at p Ss 0.01 


Group Centroid 
Unresp. 


Resp. 


Tests 


-0.44 


-0.46 


-0.46 


-0.41 


-0.60 


QO. 


0. 


44 


65 


TDP 


18.4 


20.6 


SIE}. 12) 


Canonical 
Chi-sq. Correl. 


0.40** 


O0.41** 


0.41#* 


0.39%** 


0.53** 


+ 


/ 
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Site separation was the ASFTL-P for the 0-15 cm depth. In 
addition, the other variables which were important included 
% clay, PH (15-30), growing season precipitation, and 

% O.M.. The total discriminatory power for the function was 
low (0.27) indicating a large amount of within group 
variation. The canonical correlation was also low and as a 
result a poor separation of the sites based on this function 
would be expected. 

Inclusion of agro-climatic area, soil zone or soil 
order variables into the discriminant analysis along with 
the quantitative variables did not improve the function's 
ability to separate the sites (Table 5). The low values of 
the standardized coefficients for each of the classification 
variables indicated that these variables were not very 
important in separating the sites. In addition, the total 
discriminatory power and canonical correlation are not 
Significantly improved with the inclusion of these 
classification variables. 

Two possible reasons exist for the lack of improvement 
in the discriminant function with the inclusion of the 
classification variables: either the classification has no 
Significant influence in determining barley response to 
phosphate fertilizer or, the important quantitative 
variables for discrimination differ among the classification 
classes. To check the latter possibility, individual 
discriminant analyses were conducted on each classification 


class having sufficient members. 
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Table 5. Discriminant Analyses for Barley Response to 
Phosphate Fertilizer: Quantitative and 


Classification Variables (125 sites) 


Std. Unstd. Group Centroid Canonical 
Variables Coetw. Cork. meReSspe muUnreSp a TDP @Chiz=sq,., «Correll. 
Agro-climatic Area 
Ln ASFTL-P(0-15) 0.89 Pane 
Zu clay (0-15) 0.69 0.05 
PH (15-30) 0.53 0.59 ’ 
Pptn. 0.28 0.04 
YeOoM. = 15) Seaectsy SS 
Agro-climatic Area 
1 0.18 OieZ26 
2ZAy e022 "-04 48 
2H 0.07 Osa 
3H”) 6-02.49 -0.39 
constant -9.41 -0.62 OF67 0.29 41.6 0.54% 
Soil Zone 
Ln ASFTL-P(0-15) 0.91 te dG 
Si clay ((0-15) 0.59 0.04 
pH (15-30) 0.47 0.53 
Pptn. Qeag 0.05 
nO Me iC 0=15) = (\ia 2 on mae Oat 1/0 
Soil Zone 
Gray <0 06)°=0 2110 
Black Ocoee On a2 
Dark Gray 0.02 0.35 
Dark Brown =cut O29 
constant — G25 =0s63 0.68 0.30 43.1 0.55** 
Soil Order 
Ln ASFTL-P(0-15) O95 ae 
% clay (0-15) 0.61 0.43 
PH (15-30) 0.54 0.60 
Pptn. O7s3 0.04 
Rn OLM. = (O15) = 29 pO le 
Soil Order 
Chernozemic 0.08 0.16 
Luvisolic 0.08 Osis 
Solonetzic 0.05 OaZ6 
constant =954 5-060 OF 65 0.28 24 ls! 0.53** 
** Significant at p < 0.01 
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For all three of these classifications, each individual 
class analyzed varied as to which quantitative variables 
were important for site discrimination (Table 6, 7, and 8). 
The most common variables were ASFTL-P(0-15), 
ASFTL-P(15-30), % clay, and precipitation. With two 
exceptions, the individual class disciminant analysis was 
more effective in distinguishing responsive and unresponsive 
Sites than the effect coded analyses. Both total 
discriminatory power and canonical correlation were improved 
in many instances, and as a result, a high degree of correct 
classification can be anticipated. In general, the classes 
for which it was difficult to separate sites by means of a 
discriminant function were those sites belonging to the 
agro-climatic area '1', the Black soil zone, and the 
Chernozemic soil order. For these sites there were possibly 
other site parameter(s) controlling the response of barley 
to phosphate fertilizer. These could include micronutrient 
deficiency, experimental technique, pests and/or disease. By 
contrast, the best expected prediction of site 
responsiveness would be for those sites belonging to the 
Gray soil zone (Table 7), Luvisolic soil order (Table 8), or 
agro-climatic areas '2H' and '3Ha' (Table 6). 

Examination of the standardized coefficients for each 
of the functions presented in this section indicates a 
number of general trends for the influence of the site 
variables on barley response to phosphate fertilizer. The 


importance of each variable within a given function is shown 
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Table 6. Discriminant Analyses for Barley Response to 
Phosphate Fertilizer for Four Agro-climatic Areas 


Std. Unstd. Group Centroid Canonical 
Variables Coet. scock. srespeesunresp.. (DP Chi-sq. Correl. 


Area '1' (54 sites) 


Ln ASFTL-P(0-15) ata 1263 
Ln ASFTL-P(15-30) -0.82 -0.81 
clays CO=1'5) 0.89 0.08 
% sand (0-15) 0.62 O..05 
ECoCOs me O15) Oss 2 0.98 
constant cls ee) OAS, O.63 0.29 1427 0.56** 


Areas 2A‘ Glo SilLees) 


melee, Ue), the 3 Cum eres 
Zeclay CO=15) 57 0.07 
constant Ines) 0),.76 0276 0.34 6.6 0.63%** 


Area '2H' (34 sites) 


ASFTL-P(0-15) 0.85 0.05 
Poe to= 30) 0.82 0.81 
Pptn. 0.68 0.10 
constant =O 0 y= 07,86 1.09 0.48 pa) ae 0.71%** 


Area '3Ha' (22 sites) 


Ln ASFTL-P(0-15) 0.54 0.58 


Reclay 00.5) (eee: 0.08 
Ln ASFTL-P( 15-30) Wee? 1203 
% sand (0-15) SO 85 Ol Oe 
AOA MOORS Ty, [Ore oe.) ea o 
Pptn. g.32 0.04 
constant SABES Sse: 1.49 OS 25.0 0.88%** 


** Significant at p < 0.01 





Table 7. Discriminant Analyses for Barley Response to 
Phosphate Fertilizer in Five Soil Zones 


Std. Unstd. Group Centroid 


Variables Coef. Coef. Resp. Unresp. TDP Chi-sq. 


Gray Soil Zone (24 sites) 


En ASPTL-PCTS-s0y "1 2¢€5 Peete 


% clay (0-15) 0.90 0.08 
eG. CO—=15)) Ges2 2.50 
constant Areal ko 57 0.64 24.4 


Dark Gray Soil Zone (31 sites) 


ASFTL-P( 15-30) On73 0.03 
% clay (0-15) 1.06 0.08 
pH (0-15) O'S50 0.66 
BGs CO—15) =—74 <Teark 
CaCO. (0-15) 0758 9.24 
ASFTL-P(0-15,) O26 “8 0'S04 
constant = Sloe =O 16 O32 0.40 14.5 


Black Soil Zone (31 sites) 


Pptn. 0.84 Os 
ASFTL-P(0-15) 250 0.03 
Boece (is-30) 0.45 Nie loys} 
constant = Uh eile ete sits! 105 0.36 isc 


Thin Black Soil Zone (19 sites) 


ASFTL-P( 15-30) = vern = (adit 
PH (15-30) -1.32 -3.94 
meee Mee CO 15;) 1205 0.65 
wecaco sw (O=15:) 0.89 1.80 
% sand (0-15) 0.78 0.09 
Bre vO 15:) = (53) ate 1.0 
constant Zee oe 0.61 0250 10.8 


Dark Brown Soil Zone (18 sites) 


Pptn. 0.54 Onar0 
% clay (0-15) 0.93 0.41 
BeGorGlS=30)) Oi) Senn = ieee. s 
ASFTL-P(0-15) 0.67 0.04 
Tt. Mato iol) -0.51 -0.24 
constant =O ela ic. 0.95 Oe: 1220 


** Significant at p $s 0.01 


HY 


Canonical 
Correl. 


0.81** 


0.65** 


0.63** 


0.73% 


0.75** 
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Table 8. Discriminant Analyses for Barley Response to 
Phosphate Fertilizer for Two Soil Orders 


Std. Unstd. Group Centroid Canonical 

Variables €oete)  Coek. Resp. Unreésp. TDP Chi-sq. Correll. 

Chernozemic (79 sites) , 
Ln ASFTL-P(0-15) tia eee? 
PH (15-30) 0.56 OnrT 
clay, (0-15) 0.90 0.06 
Pptn. 0.64 0.08 
ore Mee CO ho) =-0;37 -0.16 
ENeASH LEP (i530) —-0...62 —0s55 
% sand (0-15) 0.64 0.04 

constant =e Ome oS 0.59 0.26 23.4 0.53%** 


Luvisolic (41 .sites) 


Ln ASFTL-P(0-15) 0.55 0.63 
Ln ASFTL-P(15-30) 0.99 OSu7 
mee ava GO= ly) 0.89 0.08 
<7€aGo.) (0-15) 0.26 3.43 
constant sik TS) Salih Pease, 0.59 SHS}. 3" 0.78** 


** Significant at p < 0.01 





ao 


by the magnitude of the standardized coefficient, while the 
Sign, when considered along with the group centroids for 
each group, indicates the behavior of the variable in 
relation to the separation of the sites. Whenever 
ASFTL-P(0-15), % CaCO;, % clay, or growing season 
precipitation appeared in the function, the sign associated 
with the coefficient was consistently positive, while the 
value of the group centroid for the responsive group was 
lower than the unresponsive group. Therefore, as the value 
of these variables increased, the site tended to be 
unresponsive to phosphate fertilizer. Other variables that 
appeared in the various functions, but were not consistent 
in their behavior among the classification classes, included 
Betas ei 5-30) 2 opH( 15-30), B.C. , -%,OsMt hand: %) sand..‘The 
results presented indicate that the critical soil test value 
varies depending upon the soil and environmental properties 
of the site. Therefore, contrary to the approach of Cate and 
Nelson (1965), it would be difficult to use one soil test 


value as the critical soil test value for all soils. 
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2. Multiple Regression Analysis 

The yield increase for each responsive site was 
dependent on the phosphate fertilizer rate calculated to 
attain 90% maximum yield for that site. Therefore, the 
variation of this relationship among all sites would be due 
primarily to the differences of the soil and environmental 
properties. Multiple regression analysis should identify and 
quantify those variables. 

A comparison of the three soil test procedures 
indicated that ASFTL-P was the best soil test procedure in 
accounting for the greatest amount of variation in yield 
increase (Table 9)..The ASFTL-P(0-15) accounted for 15% of 
the yield increase variation with an additional 3% accounted 
by ASFTL-P(15-30). By contrast, for both M & A-P and 
Olsen-P, only the 0-15 cm. depth soil test was significant, 
accounting for only 8% and 4% of the variation respectively. 
For all three methods, the natural logarithmic form of the 
soil test was better than the untransformed soil test 
values. As a result of this comparison, the ASFTL-P was used 
in the succeeding analyses. The best combination of 
quantitative variables accounted for 73% of the variation in 
yield increase (Table 9). These variables, and the 
approximate additional variation each explained, include the 
phosphate fertilizer rate needed to reach 90% of the maximum 
Site yield (47%), ASFTL-P(0-15) (15%), ASFTL-P(15-30) (5%), 
pH (0-15) (3%), the growing season precipitation (2%), and 
the % O.M. (0-15) (2%). 
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Table 9. Stepwise Multiple Regression Analyses for Responsive 
Barley Sites: (1) Comparison of Soil Tests, and 
(2) Best Combination of Quantitative Variables 
for Yield Increase Equation 


Overall 
b Std.Err. F R? Stat Brrr ger 
Variables Value b Value Change Est. Value 
1. Comparison of Soil Tests (55 sites) 
Pro, 090% Max,vid.) 0 -23%45) 0503 54.98 0.54 
Ln ASFTL-P(0-15) -3.56** 0.65 29.91 ets 
Ln ASFTL-P( 15-30) 1.03** 0.49 4.49 0.03 
constant 9.60 SiG 41.62 
Proce (902 "Max Vlad.) 10 (25** 70-04 Bats19 0.54 
Ln M & A-P(0-15) 2 eke Oy Gs 11.34 0.08 
constant Vlas: Si sieht 42.59 
P,0, (90% Max.Yld.) 0.26** 0.04 46.70 0.54 
Ln Olsen-P(0-15) 2 SOs ee NOS, Sissi? 0.04 
constant 7.06 4.00 36% 2 


R?2 


O.7 1s 


0.62** 


0.58** 


2. Best Combination of Quantitative Variables (65 sites) 


P.O, (90% Max.Yld.) 0.19#* 0.03 3223.0 Ole Ac! 


Ln ASFTL-P(0-15) -4,.46** 0.64 48.95 OL 15 
Ln ASFTL-P( 15-30) 1.48%* 0.43 tee Gul 0.05 
pH (0-15) =1.,64%% 0.63 6.85 0.03 
Pptn. =(sttee (0.06 3.87 ‘Ohya 
Peo Ma (0-15) -0.34** 0.18 3.49 0.02 
constant 28.46 329 26.06 


** significant at p < 0.01 


0.73** 
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The influence of each variable can be determined by 
examining the magnitude and sign of its coefficients. As 
expected, the phosphate rate had a positive effect on yield 
increase, but unexpectedly, so did the ASFTL-P(15-30). The 
other variables in the analysis all had negative effects on 
yield increase, that is, the yield increase was depressed as 
the the value of these variables increased. Thus, an 
increase in the ASFTL-P(0-15) reduced the yield increase 
from applied phosphate fertilizer as would be expected if 
the ASFTL-P(0-15) provided a measure of the available 
phosphorus in the soil.-A similar trend can be observed for 
pH, precipitation, and organic matter content, except that 
the magnitude and the yield increase variation accounted by 
these variables was smaller. Even though the influence of 
pH, precipitation, and organic matter content were 
Significant, these variables combined accounted for only an 
additional 6% of the variation in yield increase to 
phosphate fertilizer. 

As noted, an increase in pH appeared to depress barley 
yield increase to phosphate fertilizer. If the same rate of 
phosphate fertilizer was required to reach "optimum" yield, 
the yield increase for an alkaline soil would be less than 
that of an acidic soil. This could possibly be due to a 
greater availability of the phosphate fertilizer under 
acidic soil conditions and/or a lower yielding potential of 
the crop on alkaline soils. Because of the lower yield 


increase on alkaline soils, the "optimum" rate of phosphate 
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fertilizer will be lower. Hallsworth (1969), commenting on 
work by Colwell and Esdaile (1968), makes a similar 
conclusion: "for soils containing say 5 ppm available P, the 
phosphate dressing for most profitable response is twice as 
high on an acid soil (pH 5.5) as it is on an alkaline soil 
(pH 8.0)". This suggests that, as the pH of a soil 
increases, the phosphate fertilizer application rate 
required to obtain the optimum yield for a site should be 
reduced. 

The influence of precipitation on crop response to 
phosphate fertilizer tended to be negative, that is, as 
precipitation increased, the crop response to phosphate 
fertilizer was reduced (Table 9). Strong and Barry (1980) 
found a relatively large crop response to phosphorus under 
arid conditions due to a reduced volume of soil exploited by 
the crop's root system and the relatively high availability 
of phosphorus in the fertilizer band. Thus under arid 
conditions, the crop made more use of the fertilizer 
phosphorus than under non-arid conditions where the crop 
made more use of native soil phosphorus. Hutcheon and Rennie 
(1960) reported a significant decrease in the availability 
of soil phosphorus to wheat as the moisture stress 
increased, and an increase in the relative availability of 
the fertilizer phosphorus banded with the seed. 

Organic matter content of the soil appeared to have a 
negative effect on yield increase (Table 9). As the organic 


matter content increased, yield increase to added phosphate 
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fertilizer was depressed. This could suggest that the crop 
obtained phosphates from organic sources through 
mineralization (Stewart et a], 1980), or that phosphate 
sorption increased with increasing organic matter (Rennie 
and McKercher, 1959; Harter, 1969; Hinga, 1973; 
Lopez-~Hernandez and Burnhan, 1974; Holford and Mattingly, 
foro). 

Having established the effects of the quantitative 
variables, the next step was to assess whether the 
prediction of yield response would be improved by including 
site classification variables. The inclusion of 
agro-climatic area, soil zone, or soil order classification 
variables into the regression procedure resulted in a small 
increase in the equation's overall correlation (Table 10). 
Agro-climatic area accounted for an additional 3% of the 
variation in barley response, with all previously determined 
covariates remaining significant. Soil zone accounted for an 
additional 5% of the variation while soil order accounted 
for an additional 4% variation. However, for both soil zone 
and soil order equations, the organic matter variable became 
nonsignificant and was discarded prior to inclusion of soil 
zone or soil order variables. Thus, the variation accounted 
for by organic matter content was accounted for by the soil 
zone or soil order variables. 

To determine whether a significant difference existed 
among the classes for each classification system, an 


approximate multiple range test was conducted on the 
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Table 10. Yield Increase Equations for 65 Responsive Barley 
Sites with Site Classification Using Stepwise 
Multiple Regression Analysis 


b 
Variables Value 


Agro-climatic Area 


Pee (90% Max. Yid.). O.17#* 


Ln ASFTL-P(0-15) = Oe ee 
Ln ASFTL-P( 15-30) 1.40** 
PH (0-15) -1.57%* 
Sfonejele -0.14*# 
MeO Me CO= hoy -0.29% 
Agro-climatic Area 
1 -0.80 
2A = laa, 
2H =Qieos 
3H 1219 
constant 20 ih 


Soil Zone 


Pow (90% Max.vlas.) 0. 20*+ 


Ln ASFTL-P(0-15) -3.95%* 
Ln ASFTL-P( 15-30) 1.52** 
Pptn. -0.17%** 
pH (0-15) -1.06** 
Soil Zone 
Gray 35S 
Black =O) ey 
Dark Gray =r. Lie 
Dark Brown ee 
constant PM Foe 
Soil Order 
PeO;. (90% Max<Yld.) 0.19% 
Ln ASFTL-P(0-15) =—4 27 £4 
Ln ASFTL-P( 15-30) 1.49% 
Pptn. -0.15** 
pH (0-15) =1.48%# 
Soil Order 
Chernozemic = atin 
Luvisolic 1.83 
Solonetzic deat 2 
constant 25.84 


4 Significantiat ps 0.05 
** Significant at p s 0.01 


Sta. srn. 


b 


F 
Value 


PVE, 
44,15 
10.09 
BeBe, 
593 
Zs 


3925 
42.24 
1391 
8.48 
Bo eve 


367,43 
49.23 
13205 
8.07 
6.21 


R? 
Change 


247 
aio 
05 
oiO:3 
02 


Cc eorea Co 


0.47 
OeS 
0.05 
0.02 
0.03 


Overall 


Std brn. 
BSt. 


F 
Value 


16.96 


Pat eas 


23.45 


R? 


0.76** 


0.78** 


0.77** 
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estimated class means (see Material and Methods). No 
Significant differences were found among the agro-climatic 
area class means, but a significant difference existed 
within both the soil zone and soil order classifications 
(Table 11). Within the soil zone classification, the Gray 
zone had a significantly greater yield increase than the 
other soil zones, with no significant difference among the 
other four zones. Within the soil order classification, 
there was a significant difference between Luvisolic and 
Chernozemic sites with the Luvisolic sites having a greater 
response to phosphate fertilizer. Solonetzic and Gleysolic 
Sites showed no significant difference from either the 
Luvisolic or Chernozemic sites, probably due to the low 
number of sites within each class and the resulting high 
Standard errors for the means. The results of the soil zone 
and soil order were in agreement with each other which might 
be expected since most Luvisolic sites were within the Gray 


soil zone. 
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Table 11. Comparison of Mean Yield Increase 
(100 kg/ha) for Responsive 
Barley Sites in Various Classes 


Classification Mean Std. Error 


Agro-climatic Area 


1 6.34 a Or75 
2A DoS eI lees 
2H 6.Sa 0.81 
3H Ses2u a Wea2 
3Ha 8.97 a eres ia) 
= 7.14 0.44 
Soil Zone 

Gray Qo a. 0.83 
Dark Gray Gasourd 0.78 
Black P SACO de’ 0.70 
Dark Brown Ne Siehe ele: Sasi, 
Thin Black 4.86 b 4.86 
= 6.56 0.44 

Soil Order 
Luvisolic 8.98 a 0.67 
Solonetzic 8.23 ab eh PAG) 
Gleysolic 6.61 ab 2e2o 
Chernozemic 5.6 ib Oc 51 
x 7341 0.99 


Means within a classification having different 
letters are significantly different (P < 0.05) 
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3. Principal Component Analysis 

The results presented in this section are those 
principal components accounting for the majority of the 
variation in the responsive site data. The independent 
variables (eigenvectors) are listed along with their factor 
loadings for each component. These loadings measure the 
degree each variable is involved in each factor pattern, and 
can be interpreted like correlation coefficients (Rummel, 
1967). 

The sum of the five largest eigenvalues explained about 
75% of the total varianée within the data (Table 12). 
Principal component number 1 accounted for about 25% of the 
variation and was heavily loaded, positively, by pH and E.C. 
variables and moderately loaded, but negatively, by 
precipitation and the phosphate fertilizer calculated for 
optimum yield. As the pH and E.C. of a soil increased, the 
phosphate fertilizer rate for "optimum" yield decreased. 
This relationship between pH and phosphate fertilizer was 
noted previously in the results of the multiple regression 
and the discriminant analyses. The relationship between E.C. 
and phosphate fertilizer could represent an effect of the 
soil solution (including NO;,-N) concentration, on crop 
utilization of phosphate fertilizer. This group of variables 
reflect the chemical potential of the soil solution and can 
be labelled the "soil solution component". 

The second principal component accounted for about 17% 


of the variation in the data, and was loaded heavily by 
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Table 12. Principal Component Analysis of Responsive 
Barley Sites: The Five Largest Eigenvalues 


Principal Component No. 1 2 3 4 5 
Eigenvalue 2.796 Uatciely 1.696 I AS 0.834 
(cumulative percentage) 25.4 42.7 58.11 68.3 Io 
Eigenvectors 

pH (0-15) pele SGs 0.3680 k 70.223 0.251 -0.064 
pH (15-30) G2e72ee—. 205° 0. 129 0.149 O,258 
E.G.) (0-15) 0.720 0.144 Oe sas Oait3s ~=0.245 
Bac.” (75-30) 0.608 OF 189 O46 7a On 26 0.245 
we ahs 60795) On SS 0.423 0.264 C2502 0.266 
™ Caco, (0-15) oes 25 ees OrOOn uO 455) m —Oalsiy, 
% clay (0-15) =@2.064 QJ225 Oeste > =0),002") -0.209 
Pptn. SOse20 Posteo: =O. 181 0.644 -0.151 
Ln ASFTL-P(0-15) 0.061 Ono 2 Ome O race One cima Otel 
Ln ASFTL-P( 15-30) 0.087 Coo. 0. oo4) —O. ten 0.378 


P,0, (90% Max.Yld.) jU439 8 =0%335 0.408 -0.088 0.547 
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ASFTL-P(0-15), moderately by ASFTL-P(15-30), carbonate 
content of the soil, organic matter content of the soil, 

pH (0-15) and the phosphate fertilizer rate. Again the 
fertilizer rate displayed a negative loading, that is it had 
an inverse relationship with the ASFTL-P variables. This can 
be interpreted as meaning that as the value of the ASFTL-P 
increased, fertilizer requirement decreased. The inverse 
relationship between organic matter content of soils and the 
fertilizer rate suggests a possible mineralization of 
Organic phosphorus to supply phosphate to the crop. This 
component can be labelled the "available phosphorus 
component”. 

The third principal component accounted for about 15% 
of the variation in the data, and was loaded heavily by the 
clay content of the soil, and moderately loaded by 
ASFTL-P(15-30), E.C.(15-30), and the optimum phosphate 
fertilizer rate. The most important relationship that should 
be noted was between clay content of the soil and phosphate 
fertilizer rate. As clay content increased, the phosphate 
fertilizer rate required for 90% of maximum yield also 
increased suggesting an adsorption reaction between clay and 
phosphate fertilizer. Component 3 could be labelled the 
"phosphate adsorption component". 

Principal component number 4 was composed primarily of 
precipitation, organic matter content of the soil, and 
carbonate content of the soil and accounted for about 10% of 


the variation. Precipitation and organic matter content had 
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positive loadings, while carbonate content of the soil had a 
negative loading. Thus, as precipitation increased, organic 
matter also increased, while carbonate content decreased. 
This relationship reflects the trend seen in the soil and 
climate for the province as one goes from the Brown soil 
zone to the Black soil zone. Therefore, this component can 
be labelled as the "soil zonation component". 

The final component, number 5, accounted for about 7% 
of the variation, and was loaded moderately by the 
ASFTL-P(15-30) and the optimum fertilizer rate, both 
positively. The positive loadings of these two variables in 
this component correspond to the positive coefficients 
observed in the multiple regression analysis. Since the 
majority of the crop roots are found in the top 15 cm of the 
soil, utilization of the available phosphorus in the second 
15 cm could have an effect similar to fertilizer on crop 
response. This component can therefore be labelled the 
"phosphate fertilizer rate component". 

The principal component analysis served to illustrate 
the complex interrelationships among soil variables and in 
particular the relationships between "optimum" phosphate 
fertilizer rate and other soil properties. In order to use 
the soil test as the criterion for predicting phosphate 
fertilizer requirements for optimum yield, the relationships 
between phosphate fertilizer rate for optimum yield and soil 


properties must be taken into account. 
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4. Summary 

The results of the diScriminant analyses of the barley 
Sites indicated that there was very little difference among 
the three soil test procedures. Overall, the most important 
quantitative site variable determining the response of 
barley to phosphate fertilizer was the soil test for 
phosphorus; as this variable increased, the site tended to 
be unresponsive. Other variables which occurred commonly in 
the discriminant functions of the barley sites included 
% clay, % CaCO;, and growing season precipitation. A 
quantitative increase of any of these site variables tended 
to categorize a site as unresponsive to phosphate 
fertilizer. Site classification did enhance the separation 
of the sites, but only when individual classification class 
discriminant functions were determined. 

Multiple regression analysis of the responsive barley 
Sites indicated that the soil test best accounting for the 
variation in yield increase was the ASFTL-P. The best 
combination of quantitative site variables that were 
Significant in accounting for the variation in yield 
increase to phosphate fertilizer was the calculated optimum 
fertilizer rate, ASFPTL-PCO-15),; ASFTL-P(15-30), pH (0-15), 
growing season precipitation, and the organic matter content 
of the soil. Yield increase from phosphate fertilizer was 
depressed by an increase in ASFTL-P(0-15), pH (0-15), 
precipitation, and % organic matter. Inclusion of site 


classification variables into the analysis did improve the 
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prediction ability of the equation. Multiple comparison 
tests of the estimated means indicated no significant 
difference among the agro-climatic areas, but within the 
soil zone and soil order classifications, significant 
differences existed among the class means. The Gray soil 
zone and the Luvisolic soil order were significantly more 
responsive to phosphate fertilizer than the remaining 
classification classes. 

Principal component analysis of the responsive sites 
illustrated the complex interrelationships among the site 
properties, and with the calculated optimum fertilizer rate. 
The site variables measured can be reduced to five 
components representing (1) the soil solution, (2) the 
available phosphorus, (3) phosphate adsorption, (4) soil 
zonation, and (5) the phosphate fertilizer rate. The most 
noteworthy interrelationships were the inverse relations of 
"optimum" phosphate fertilizer rate and each of pH, ASFTL-P, 
and organic matter, and the direct relation between clay 
content and phosphate fertilizer rate within certain 


components. 
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B. Rapeseed 

The results presented in this section represent the 
Statistical analyses of the yield response of rapeseed to 
phosphate fertilizer for 91 sites. A brief summary (means, 
Standard deviations, maximum and minimum values) of the 
chemical and physical characteristics of the field sites is 
presented in Tables 13 and 14, while a frequency 
distribution of the sites in each of three site 
classifications is presented in Table 15. In general, the 
Sites used in this study represented a wide variety of site 
conditions for both responsive and unresponsive groups. The 
most noteworthy difference between the two groups was a 
lower mean value of the soil test levels of phosphorus for 
the responsive sites, especially in the surface depth and a 
higher mean preciptation for the responsive sites. The 
distribution of responsive and unresponsive sites among the 
classification classes was unequal for some classes 
(Table 15) suggesting a greater importance of site 
classification in the separation of rapeseed sites than that 
observed for the barley sites. Finally, the general 
geographic distribution of the sites in the province was 
representative of the major dryland crop producing areas of 


the province as well as the major soil and climatic groups. 
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Table 13. Mean, Standard Deviation, Maximum, and Minimum 
Values of the Independent Variables for the 


Unresponsive Rapeseed Sites 


Variables* 


PH (0-15) 

PH (15-30) 

Ee Ce uG=1s) 

Be Ger 1530) 
%-O.M. (0-15) 

% Caco; (0-15) 

% sand (0-15) 

yw srlt 0-15) 
%.clayi(0-15) 
Pptn. 
ASFTL-P(0-15) 

Ln ASFTL-P(0-15) 
ASFTL-P( 15-30) 

Ln ASFTL-P( 15-30) 
M & A-P(0-15) 

Ln M & A-P(0-15) 
M & A-P( 15-30) 

Ln M & A-P( 15-30) 
Olsen-P(0-15) 

Ln Olsen-P(0-15) 
Olsen-P( 15-30) 

Ln Olsen-P( 15-30) 


* Variable 


Bae. 
Pptn. 
ASFTL-P 
M & A-P 
Olsen-P 


Mean 


6.87 
Vos 
0.40 
O25 
5520 
0.25 
30.82 


poveet 


PALI ONS 
Tei 
Ce 4) 
3.68 
14.8 
1.68 
Dos 
3.80 
Iho 
PES 
41.1 
3.61 
USez/ 
Pais hy) 


Std. Dev. 


DUS 
Deel: 
Denes 
0.68 
2.47 
Oren 
18.00 
12.60 
11.40 
8.03 
40.1 
0.68 


Units 


mmhos/cm? 
cm 

kg/ha 
kg/ha 
kg/ha 


Max. 


_ 
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NP WwW fF OWMM 
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Cie Doc Wat 


Min. 
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Table 14. Mean, Standard Deviation, Maximum, and Minimum 
Values of the Independent Variables for the 
Responsive Rapeseed Sites 


No. of 
Variables+* Mean Std. Dev. Max. Min. Sites 
pH (0-15) 6.69 Vie he 8.1 5.2 52 
PH (15-30) 6.74 0.91 8.2 4.6 52 
Be Ce O— 5) 0.34 0.14 0.8 0.2 52 
Bees Gt5=—3.0) 0.36 0.25 15 Onn 52 
Ta0-MeecO—15) Sig y/ 2.58 1427 2 52 
% Caco, (0-15) 0.19 AR2 Bia 0.0 52 
% sand (0-15) 3/3/1500 19.14 UNS 269 52 
Tost Lte (O—t5) SS Geneon ils 26: 5952 t2e2 52 
% clay (0-15) 20.21 14.10 Tileo Binz 52 
Pptn. 2 ie Od IAI Sine 5.8 52 
’ ASFTL-P(0-15) 28.0 26.8 134.4 0.0 52 
Ln ASFTL-P(0-15) 2.97 0.90 4.90 0.0 52 
ASFTL-P( 15-30) £5. SiS ys! 201.6 0.0 52 
Ln ASFTL-P(15-30) 1.65 ese aim Eh OO 52 
M & A-P(0-15) 29.03 26.0 131.0 4.5 52 
Ln M & A-P(0-15) 3.08 0.78 4.88 1350 50 
M & A-P(15-30) 17.4 34.5 Z00% 5 0.0 50 
Ln M & A-P( 15-30) 2.09 tet 5.30 0.0 50 
Olsen-P(0-15) 24.4 225 G2 thes! 50 
Ln Olsen-P(0-15) 307 0.52 4.14 2.06 50 
Olsen-P( 15-30) IAS 14.9 89.6 4.5 50 
Ln Olsen-P( 15-30) 262 0.65 4.50 1250 50 
P,0,-90% (Max.Yld.) 30.6 46.7 Givin PLEO 52 
* Variable Units 
EaCs mmhos/cm? 
PpEns cm 
ASFTL-P kg/ha 
M & A-P kg/ha 
Olsen-P kg/ha 


P,0, (90% Max.Yld.) kg/ha 





Table 15. Frequency Distribution of Unresponsive 
and Responsive Rapeseed Sites per 
Classification Class (Number of Sites) 


Classification Unresponsive Responsive 


Agro-climatic Area 


1 16 13 
2A 7 0 
2H 12 16 
3H 0 6 
3Ha 4 iy) 
Soil Zone 
Gray 3 14 
Dark , 6 24 
Black 11 1 
Thin Black 141 1 
Dark Brown 5 0 
Brown 3 0 
Soil Order 
Chernozemic 29 27 
Luvisolic 8 23 
Gleysolic 0 1 
Solonetzic 2 1 
Texture (0-15) 
HC 2 3 
( 0 2 
Sic 3 J 
SiC 5 8 
CL 11 14 
SCL 2 0 
SiL 2 4 
L 9 11 
SL 5 8 
LS 0 1 
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1. Discriminant Analyses 

The objective of this series of analyses was to 
determine those site variables important for separating 
responsive sites from unresponsive sites for the phosphate 
fertilizer response of rapeseed. Simple examination of the 
mean values (Tables 13 and 14) indicated that the major 
difference between the two groups was the soil test for 
available phosphorus. 

There was a Slight difference among the three soil test 
procedures in separating responsive and unresponsive sites 
(Tables 16). M & A-P and Olsen-P had the greatest success in 
distinguishing these groups. For all three procedures, both 
depths were important in the function. There was no 
improvement in the correlation of the discriminant function 
when a larger data set was used (Table 16). Since there is a 
close procedural relationship between M & A-P and ASFTL-P, 
but there was a larger sample population having ASFTL-P 
information, ASFTL-P was used to determine the best overall 
function (Table 16). In addition to the ASFTL-P, the other 
variables which were important for site distinction were 
PH (15-30), and growing season precipitation. However, the 
effectiveness of the function to separate sites was still 
poor as indicated by the low total discriminatory power and 
canonical correlation. Inclusion of soil order into the 
analysis with the quantitative variables did not improve the 
ability of the function to separate sites, however, 


inclusion of agro-climatic area or soil zone did improve the 
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Table 16. Discriminant Analyses for Rapeseed Response to Phosphate 
Fertilizer: (1) Comparison of Soil Tests, and 
(2) Best Overall Function 


Std. Unstd. Group Centroid Canonical 
Variables Cock. a CoctemmResp amuUnrespe TDPEEChissq. | Correll. 


1. Comparison of Soil Tests 


(86 sites) 
ASFTL-P(0-15) 81 0.05 
ASFTL-P( 15-30) =1.52 =0.04 
constant ae 4 OS O73 O27 27 6 0.53%** 


Ln M & A-P(0-15) 1.61 2.21 
Ln M & A-P(15-30) -1.24 -1.09 


constant = ne OD 0.90 0.36 38.9 CO. 61ee 
Olsen-P(0-15) 1a25 0.08 
Olsen-P( 15-30) =O Sam O 0.6 
constant = 160. A061 0.85 0.34 B52 0.59** 
(91 sites) 
ASFTL-P(0-15) ene 0.05 
ASFTL-P( 15-30) -1.46 -0.04 
constant =e Ole oO 0.67 O525 26.0 0.51** 


2. Best Overall Function (91 sites) 


Ln ASFTL-P(0-15) leer sie Warez 

ENGASETE=E (15-30) -0 2945 -0474 

pH (15-30) 0258 0.62 

Pptn. =) 300 9-0, 05s. 

constant =e 34 an Ole OD 0.87 0.36 40.0 0.61** 


** significant at p < 0.01 
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function's correlation, with the latter classification 
showing the greatest improvement (Table 17). 

As with the barley sites, the individual classes for 
each classification were analyzed to determine variable 
differences among the classes for site discrimination. 
Results indicated that not only did the quantitative 
variables vary among the classes, but so did their 
importance and behavior (Tables 18, 19, and 20). The most 
common variables were the ASFTL-P(0-15) and ASFTL-P(15-30), 
while the other variables seem to appear at random in the 
functions. In general, the individual class discriminant 
analysis was more effective in separating responsive and 
unresponsive sites than the effect coded analysis. Both 
total discriminatory power and canonical correlation were 
improved in many instances for the individual class analyses 
and as a result, a high degree of correct classification of 
the sites could be expected. The exceptions were 
agro-climatic area '2H' and the black soil zone. Thus, it 
would appear that some additional unmeasured site 
parameter(s) was controlling phosphate response of rapeseed 
in these classification classes. These could include 
deficiency of micronutrients, experimental technique of the 
individual coordinators, pests, and/or disease. 

The behavior of the site variables is determined by 
examination of the standardized coefficients of the 
discriminant function in relation to the group centroids. 


Only the ASFTL-P for both depths, E.C. for both depths and 
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Table 17. 


Discriminant Analyses for Rapeseed Response to 


Phosphate Fertilizer: Quantitative and 
Site Classification Variables (91 sites) 


Variables 


Unresp. TDP Chi-sq. 


Std: 
Coef, 


Unstd. 


Coef. 


Agro-climatic Area 


Ln ASFTL-P(0-15) 
Ln ASFTL-P( 15-30) 
pH (15-30) 
Pptn. 
Agro-climatic Area 
1 
2A 
2H 
3H 
constant 


Soil Zone 


Ln ASFTL-P(0-15) 
Ln ASFTL-P( 15-30) 
pH (15-30) 
Pptn. 
Soil Zone 
Gray 
Black 
Dark Gray 
Dark Brown 
Brown 
constant 


Soil Orde 


Ln ASFTL-P(0-15) 
Ln ASFTL-P( 15-30) 
pH (15-30) 

Pptn. 

Soil Order 
Gleysolic 
Luvisolic 
Solonetzic 


constant 


Tas 
=0 569 
disse 
=Oe22 


0.16 
0.79 
O13 
-0.42 


=0.98 

Oma 
-0.08 
-0.08 


0.83 
0.30 
Oeo7 
=O 59 
=0.67 


1S 


e2o 
ai Over oul 
0.54 
SW) seié: 


= Whose 
0.09 
0.01 


jects, 
=0553 
0.34 
=0203 


0°22 
eos 
eles RE, 
=0°80 
=i 33 


= lereo 
0.54 
-0.08 
-0.01 


1.60 
reo 0 
1.64 
= ees 
Sha tie 
3.04 


12 o8 
=ORoo 
0.58 
0.04 


-0.66 
O02 
0.02 

Ther 


** significant at p < 0.01 


Group Centroid 


Resp. 


Ord 


= Wid AG, Ono2 63—0 


41.2 


Canonical 
Correl. 


0.67** 


0.73%** 


0.62** 
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Table 18. Disciminant Analyses 
Phosphate Fertilizer 


Variables 


Area '1' 
ASFTL-P(0-15) 
Eoce (015) 
am CacOse (0=15) 
% sand (0-15) 
% clay (0-15) 
Ae OsMee (0-15) 
Beco tt5~30) 
ASFTL=P(15=30) 
constant 

Area '2H' 
Ln ASFTL-P(0-15) 
PH (0-15) 
Ln ASFTL-P( 15-30) 
constant 


Area 


Ln ASFTL-P(0-15) 
pH (0-15) 

Ln ASFTL-P( 15-30) 
pH (15-30) 

B.co) (15-30) 
Pptn. 

constant 


'3Ha' 


Std. unsca. 
Coef. Coef. 

(29 sites) 
= 1. 15 SOO 
=teo? MM6Gcs2 
=OC6 2 a oeeo 
Pi ST} 0.20 
2.19 0.34 
=O547 2 = Oe2e 
0.55 O71 
Overs 0.05 
=TOs27 

(28 sites) 
pe2 101883 
Qor7 0.99 
== 68.) —-0-56 
sahil Hee axes 
(21 sites) 
=—390) 94,24 
BB SY 5.66 
2.46 1.48 
SSS) erst, 
0.99 2ES7 
91 2 ae oO eek 
so 


** Significant at p < 0.01 


Group Centroid 


Resp. 


1.98 


-0.47 


1 


03 


for Rapeseed Response to 
for Three Agro-climatic Areas 


Canonical 

Unresp. TDP Chi-sq. Correl. 
=. 0.76 Se | 0.88** 
0.63 0.21 6.8 0.49** 
=47539 OR 82 Zoe 0.91** 


82 





Table 19. Discriminant Analyses for Rapeseed Response to 
Phosphate Fertilizer for Three Soil Zones 


Std. Unstd. Group Centroid Canonical 
Variables Coef. Coef. Resp. Unresp. TDP Chi-sq. Correl. 


Gray Soil Zone (17 sites) 


PH (0-15) -2.77 -4.92 
ASFTL-P(0-15) Ze OS05 
Pptn. bets 0.18 
% sand (0-15) =—ase Ol Oi. 
PH (15-30) 0.74 0.74 
constant 23 5 Se le eal 5.603 OR 87 2 0 0.94% 


Dark Gray Soil Zone (30 sites) 


Ln ASFTL-P(0-15) 1.49 2.05 
Ln ASFTL-P(15-30) -0.69 -0.54 


pH (0-15) 2747.) te. 09 
PH (15-30) =e 4 = oD 
Rec uu \5)) 0.54 3.95 
% sand (0-15) -0.86 -0.05 
Seclay (O- 15) =() O30 0605 
constant 2 Ose =O 2 S329 0.73 a5 0.86** 


Black Soil Zone (24 sites) 


ASFTL-P( 15-30) 1.06 Ostz 
ASFTL-P(0-15) -0.69 -0.04 
% clay (0-15) 0.61 Of. 05 
constant =p ee! Oo Oo Ore 6.0 0.50** 


** significant at p s 0.01 
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Table 20. Discriminant Analyses for Rapeseed Response to 
Phosphate Fertilizer for Two Soil Orders 


Std. Unstd. Group Centroid Canonical 
Variables Coef. Coef. Resp. Unresp. TDP Chi-sq. Correl. 


Chernozemic (56 sites) 


PH (15-30) 0.75 0.94 
ASFTL-P(0-15) 0.94 0.06 
ASFTL-P( 15-30) -0.57 -0.07 
% sand (0-15) (AeA 0.02 
Pptn. -0.28 -0.04 
Pec U0— 115) 0.38 2.49 
%2Caco; (0-15) =O)3 Ole = Ome e 
Fee c15o=30.) = Ole 2 Gian ale ee: 
constant =e) Cee Oe Ost 0.47 SZes 0.69%** 


Luvisolic (31 sites) 


ASFTL-P(0-15) 2200 0,05 
PH (0-15) =(eW6eu-s leo 
ASFTL-P( 15-30) =1.30 —0.02 
Pptn. 0.78 Oreilal 
% sand (0-15) Sas Stal tays: 
EeG wet 5-30) =i Ome ais 
constant 7.34 -0.93 PLS Oval SSIA 0.85% 


** Significant at p < 0.01 
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% clay exhibit a constant behavior when all functions were 
examined. The behavior of the ASFTL-P(0-15) and 
ASFTL-P(15-30) were opposite so that as ASFTL-P(0-15) 
increased, the site tended to be unresponsive, while an 
increase in the ASFTL-P(15-30) tended to result in the site 
being responsive to phosphate fertilizer. E.C. also 
displayed this type of behavior; as the E.C.(0-15) 
increased, the site tended to be unresponsive while an 
increase in the E.C.(15-30) tended to result in the site 
being responsive. The % clay of the soil appeared only ina 
few functions but where it did, an increase resulted in the 
site tending to be responsive to phosphate fertilizer. The 
remaining variables appearing in the functions were 
inconsistant in their behavior probably due to the small 


Sample sizes used in many analyses. 
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2. Multiple Regression Analysis 

As in the case of barley, for each responsive rapeseed 
Site, yield increase was dependent upon the calculated 
phosphate fertilizer rate to produce a yield that was 90% of 
the maximum site yield. Therefore, the variation in this 
relationship from site to site should be due to the 
differences among the sites' soil properties. Multiple 
regression analysis techniques were used to identify those 
variables responsible for this among-site variation in crop 
response. 

A stepwise multiple regression analysis was first 
computed for each of three soil test procedures, after which 
the best combination of quantitative variables to explain 
the variation in yield increase was determined. The 
comparison of the three soil test procedures indicated that 
there was very little difference among them (Table 21). 
ASFTL-P(0-15) accounted for 6% of the yield increase 
variation, M & A-P(0-15) explained 5% of the variation while 
Olsen-P(0-15) accounted for 9% of the variation. For all 
three methods, the natural logarithmic transformation of the 
soil test accounted for a larger portion of the variation 
than the untransformed values. Also, the soil test for the 
15-30 cm depth was nonsignificant for all methods. Since 
there was very little difference among the soil test 
procedures, and because of a larger sample size, the ASFTL-P 
was used in the next step to determine the best combination 


of quantitative variables. The stepwise procedure indicated 
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Table 21 


. Stepwise Multiple Regression Analyses for Responsive 
Rapeseed Sites: (1) Comparison of Soil Tests, and 
(2) Best Combination of Quantitative Variables 
for Yield Increase Equation 


Overall 
b Std.Err. F 13h Stdikrr sar 
Variables Value b Value Change Est. Value 
1. Comparison of Soil Tests (50 sites) 
P,0, (90% Max.Yld.) 0.08** 0.02 29.35 0.47 
Ln ASFTL-P(0-15) -0.68** 0.28 5.94 0.06 
constant 2.60 1.66 26.20 
P.O, (90% Max.Yld.) 0.08** 0.02 MEET 0.47 
Ln M & A-P(0-15) =O. 71¥*. 0.34 4.51 0.05 
constant Zak 1.68 24.86 
P,0O, (90% Max.Yld.) 0.08** 0.02 33.88 0.47 
Ln Olsen-P(0-15) -1.40** 0.47 9.10 0.09 
constant 4.84 1.61 295418 


0.53** 


0.51** 


0.55** 


2. Best Combination of Quantitative Variables (52 sites) 


P,0, (90% Max.Yld.) 0.09#* 0.01 40.18 0.48 


Ln ASFTL-P(0-15) -0.79%** 0.27 8.88 0.06 
%.€aco, (0-15) = 253e« ‘0220 6.85 0.06 
constant 2.88 ewes! PASE Ps 


** significant at p < 0.01 


0.59** 
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that the best combination of quantitative variables 
accounted for 59% of the variation in yield increase of 
rapeseed (Table 21). These variables and the approximate 
amount of additional variation each explained were the 
phosphate rate to attain 90% maximum yield (48%), 
BOpasrtoh-Pi 01715). (6%); °and the % Caco;(0-15)° (6%). 

The coefficients for the quantitative variables 
indicated the specific influence each variable had on 
rapeseed response to phosphate fertilizer. The phosphate 
fertilizer rate for 90% maximum yield had a positive 
influence, so that as the phosphate fertilizer rate 
increased, yield response also increased. Meanwhile, both 
ASFTL-P(0-15) and % CaCO; had negative influences, so that 
as the value of these variables increased, the yield 
increase was depressed. This result suggested that the soil 
test procedure did provide an index of the amount of plant 
available phosphorus present in the soil, and as this 
measure increased, less fertilizer phosphate was required to 
attain the optimum yield. The negative influence of CaCO, 
indicated a possible chemical precipitation and/or 
adsorption reaction of the added phosphates by Caco, 
reducing the availability of the added phosphate (Thomas and 
Peaslee, 1973). Since the presence of carbonates is 
resticted to alkaline soils, the results of the analysis 
would tend to support the statement made by Hallsworth 


(1969) referred to earlier in the chapter. 
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To determine if knowledge of agro-climatic area, soil 
zone, or soil order could improve yield response prediction, 
effect coded variables of these site classifications were 
entered into the multiple regression analysis. The inclusion 
of these classification variables did not significantly 
improve the explanation of the yield response variation of 
rapeseed to phosphate fertilizer (Table 22). As a result, 
estimated means for the classification classes and the 


corresponding multiple range test were not calculated. 







sic b Sabdvg2 en of 
+ , said fiom 19 \aneee 
| s¥ bebos soattey 

sdy cant Canedi 


74685 





Table 22. Yield Increase Equations for 52 Responsive Rapeseed 
Sites with Site Classification Using Stepwise 
Multiple Regression Analysis 


Overall 
b Std.Err. F Re Std serra 
Variables Value b Value Change Est. Value 
Agro-climatic Area 
P.O, (90% Max.Yld.) 0.09#* 0.02 34.91 0.48 
Ln ASFTL-P(0-15) -0.76** 0.29 6.86 0.06 
% Caco, (0-15) -0.51#* O21 Dein 0.06 
Agro-climatic Area 
1 =a 0.41 
2H =Onuto 0.39 
3H 0.43 0.54 
constant 2.87 re 62 se 
Soil Zone 
P,0, (90% Max.Yld.) 0.09** 0.02 325/37 0.48 
Ln ASFTL-P(0-15) -0.82** 0.28 8.68 0.06 
% Caco, (0-15) -0.54** 0.22 6.3 0.06 
Soil Zone 
Gray ORS 2 
Black 0.23 
Dark Gray O97 
constant Zee haRL 11.34 
Soil Order 
P,0, (90% Max.Yld.) 0.09** 0.02 Sabie 0.48 
Ln ASFTL-P(0-15) -0.86** 0.27 9.89 0.06 
% Caco, (0-15) SOSH G21 6.86 0.06 
Soil Order 
Chernozemic Sig per Sat 
Luvisolic 0.63 0.63 
Solonetzic 0.52 eal 


constant PJ oys) 1.60 ene 


** Significant at p < 0.01 


R? 


0.60** 


0.60** 


0.61** 
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3. Principal Component Analysis 

To determine the interrelationships among the 
independent site variables of the responsive rapeseed sites, 
principal component analysis was conducted. The sum of the 
five largest components explained about 82% of the total 
variance of the data (Table 23). 

Principal component number 1 accounted for about 25% of 
the variation and was heavily loaded by pH, E.C., and CaCo,, 
with a moderate loading by % clay and ASFTL-P(0-15). Of 
these variables, only ASFTL-P(0-15) had a negative effect 
while the other variables had positive effects so that there 
waS an inverse relationship between soil test phosphorus and 
the other major variables of this component. As pH 
increased, the availability of soil phosphorus, as measured 
by the soil test, decreased. This could be due to several 
reasons, including the nature of the chemical extracting 
procedure and a lower concentration of readily available 
phosphorus in the soil solution and on the soil colloids. As 
the E.C. increased, indicating a greater ionic concentration 
in the soil solution, the availability of the soil 
phosphorus decreased possibly due to chemical precipitation 
reactions with cations in solution. The inverse relation 
between the soil test phosphorus and carbonates or clay 
content could reflect adsorption equilibrium reations of 
soil phosphates with carbonates and clay particles. This 
component represents the soil solution equilibrium and can 


be labelled the "soil solution component". 
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Table 23. Principal Component Analysis of Responsive 
Rapeseed Sites: The Five Largest Eigenvalues 


Principal Component No. 1 2 3 4 5 
Eigenvalue Zee 2.401 foo 0 Vests 0.923 
(cumulative percentage) 25.6 47.4 61.8 73.8 B22 
Eligenvectors 

pe A O=15) 0.614 Joo ome 6 2 ome ONO m= Ol. 224 
Pee (15-30) 0.649 Oo 0 9052 Smee O22 See — Oi, HOt 
E.G. (0-15) Oe 9S 0.198 Oli Ol2ss 0.214 
BiG. 615=30) 0.644 0.218 0.494 Oo ts 0.398 
2 Mo O- 15) 0.165 0.084 Oe) SO sis) eee 
= Caco, (0-15) OeoZ20 sO. 200 07. 280 On4360e —0)197 
% clay (0-15) Oe 6 5 See 0 5710 0.416 each) Sle itets| 
Pptn. SOROS 05549 On2 72 ea 54) on WA) 
Ln ASFTL-P(0-15) -0.446 0.545 0.247 e570 0.147 
Ln ASFTL-P( 15-30) -0.294 Oe Sita 0s Olt 0.380 0.053 


P,0, (90% Max.Yld.) 0.243 -0.415 -0.342 -0.186 0.648 





= 


The second principal component accounted for about 22% 
of the variation in the independent site data, and was 
heavily loaded by pH, clay content of the soil, 
precipitation, and ASFTL-P, and moderately loaded by the 
phosphate fertilizer rate. The fertilizer rate exhibited an 
inverse relationship with soil test phosphorus, 
precipitation, and pH, and a direct relationship with clay 
content. The soil test apparently provided some measure of 
the amount of soil phosphorus available to the plant since 
with an increase of the soil test, there was a decreased 
need for fertilizer phosphorus, as indicated by their 
inverse relationship in this component. There was an inverse 
relationship between the phosphate fertilizer rate and 
precipitation. This relationship might occur because growing 
season precipitation would tend to increase root development 
of the crop and a greater volume of soil would be utilized 
by the crop to obtain nutrients. As a result, added 
phosphate fertilizer may not have been used as extensively 
by the crop as it would be under arid conditions (Strong and 
Barry, 1980). Hallsworth (1969) suggested a greater need of 
fertilizer phosphate under acidic conditions. A similar 
result appeared in the second component as an inverse 
relationship between pH and the fertilizer rate; that is, as 
pH increased, the optimum fertilizer rate decreased. Since 
phosphate sorption generally increases with clay content, 
one might expect the phosphate fertilizer requirement to 


increase directly with clay content, as was observed in this 
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component. This component can be labelled the "available 
phosphorus component". 

Principal component number 3 accounted for about 14% of 
the independent variable variation and was heavily loaded by 
Organic matter content of the soil, and moderately by a 
number of other variables, the most noteworthy being the 
phosphate fertilizer rate. There was an inverse relationship 
between organic matter and fertilizer rate so that as 
Organic matter increased fertilizer rate decreased 
indicating a possible mineralization of organic phosphate to 
Satisfy crop requirements. Therefore this component can be 
labelled the "soil organic matter component". 

The fourth component accounted for about 12% of the 
variation of the independent variables. The most important 
variables in this component were organic matter and 
carbonate content of the soil, precipitation, and soil test 
phosphorus. Organic matter content, carbonate content, and 
precipitation illustrated a soil zone relationship. In 
general, aS precipitation decreases, organic matter content 
of soils also decreases, while carbonate content of the soil 
increases. It 1S possible that soil test phosphorus may also 
follow a zonal trend. Thus, this component could be labelled 
the "soil zone component". 

The fifth component accounted for about 8% of the 
variation in the independent variables and was controlled 
primarily by the phosphate fertilizer rate for optimum 


yield. Therefore, this component can be labelled the 





3S 


"phosphate fertilizer component”. The loading of phosphate 
fertilizer rate as the only variable in this component would 
tend to suggest that there are other undetermined site 
variables which may influence the optimum phosphate 
fertilizer rate. 

The complex relationships among the independent 
variables for the responsive rapeseed sites were illustrated 
by this analysis. The variation of the phosphate fertilizer 
for optimum crop response was related to many soil 
properties and environmental conditions which control the 


phosphate supply to the crop. 
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4, Summary 

The results of the discriminant analyses of the 
rapeseed sites indicated that there was a slight difference 
among the soil test procedures. The M & A-P appeared to best 
separate the sites, but ASFTL-P had the advantage of a 
larger number of sites available for analysis. Overall, the 
most important quantitative site variable that separated 
responSive and unresponsive sites were the ASFTL-P tests for 
the 0-15 cm and the 15-30 cm depth. As ASFTL-P(0-15) 
increased the site tended to be unresponsive, whereas, 
ASFTL-P(15-30) had the opposite effect. Other variables 
which were consistent in their behavior in the various 
functions included E.C. for both depths, and clay content. 
Site classification did influence the separation of the 
Sites. Inclusion of either agro-climatic area or soil zone 
variables into the function improved the correlation. 
Individual classification class discriminant functions 
provided potentially the most effective means of separating 
sites. 

Multiple regression analysis of the responsive rapeseed 
Sites indicated very little difference among the three soil 
test procedures in accounting for the variation in yield 
increase from phosphate fertilizer. The best combination of 
Significant quantitative variables was ASFTL-P(0-15) and 
% CaCO,. Yield increase was depressed by an increase of 
either or both of these variables. Inclusion of site 


classification variables into the analysis did not improve 
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the equation's prediction ability. 

Principal component analysis of the responsive sites 
illustrated the complex interrelationships among the site 
properties, and with the calculated optimum fertilizer rate. 
The site variables measured can be reduced to five 
components representing (1) the soil solution, (2) the 
available phosphorus, (3) the soil organic matter, (4) soil 
zone, and (5) the phosphate fertilizer rate. The most 
noteworthy relationships were the inverse relationships 
between phosphate fertilizer rate and each of pH, ASFTL-P, 
precipitation, and soil organic matter content, and the 
direct relationship between clay content and phosphate 


fertilizer rate within certain components. 
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C. Wheat 

The wheat sites used in this study were outside the 
RAYP project but were used aS experimental sites during the 
Same period of time in a project having similar objectives. 
This project was designed to determine the response of wheat 
to phosphate fertilizer on Chernozemic and Solonetzic soil 
orders. The results presented in this section represent the 
Statistical analysis of wheat response to phosphate 
fertilizer for 38 sites. A brief description of the site 
charateristics (means, standard deviations, maximum and 
minimum values) are presented in Tables 24 and 25, while 
frequency distribution of the sites according to site 
classification are presented in Table 26. In general, the 
Sites were restricted to acidic pH values and to only a few 
classification classes. In addition, site chemical and 
physical data are available for only the 0-15 cm soil depth. 
The major difference between the unresponsive and responsive 
groups was a higher mean soil test level for phosphorus and 
a higher mean precipitation for the unresponsive sites. 
Finally, the general distribution of the sites was 


restricted to the east-central portion of the province. 
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Table 24. Mean, Standard Deviation, Maximum, and Minimum 
Values of the Independent Variables for the 
Unresponsive Wheat Sites 





No. of 
Variables* Mean Std. Dev. Max. Min. Sites 
PH (0-15) 5.70 Om35 6.3 Sar 13 
E.C. (0-15) 0.43 Ons ORT ess 13 
% O.M. (0-15) 6.65 1.64 Ort Sie iS 
% sand (0-15) 325439 TTA 44,5 20.8 13 
eC silte(0= 15) 40.12 4.21 47.6 34°57 Vs 
% clay (0-15) v dy BBY 5.59 38.6 18.5 13 
Pptn. 30.54 10.61 44.7 aeons 13 
ASFTL-P(0-15) 6s5 352.0 116.5 22.4 ils! 
Ln ASFTL-P(0-15) 4.00 0.59 4.76 Siac 1S! 
M & A-P(0-15) 68.7 34.5 128.8 25.8 13 
Ln M & A-P(0-15) Ave) 0).52 4.86 325 13 
Olsen-P(0-15) ase 6.01 TA Ae) 22.4 iz! 
Ln Olsen-P(0-15) SO Ge s7 an27 Srelit 13 
* Variable Units 

E.G. mmhos/cm? 

Pptn. cm 

ASFTL-P kg/ha 

M & A-P kg/ha 


Olsen-P kg/ha 








Table 25. Mean, Standard Deviation, 
Values of the Independent Variables for the 


Responsive Wheat Sites 


Variables 


PH (0-15) 

E.Co (0-15) 

O.M. (0=15) 
sand (0-15) 
silt (0-15) 
clay (0-15) 
Pptn. 
ASFTL-P(0-15) 

Ln ASFTL-P(0-15) 
M & A-P(0-15) 

Ln M & A-P(0-15) 
Olsen-P(0-15) 

Ln Olsen-P(0-15) 
P,0, (90% Max.Yld.) 


se sf af ae 


* Variable 


Bas 

Pptn. 

ASFTL-P 

M & A-P 
Olsen-P 

P,0, (90% Max. 


Mean 


5.63 
0.40 
Devt 
40.29 
34.63 


vas) ), 


Std. Dev. Max. 


0.28 é.3 
0.14 0.7 
1.60 8.0 
10.04 Seo) 
6.12 41.1 


Sioek: 37.4 
7.49 40.4 


Gres 70.6 
0.58 4,26 
16.53 84.0 
0.42 4.43 
S58 49.3 
0.34 3290 


Quel 106.4 


Units 


mmhos/cm? 
cm 

kg/ha 
kg/ha 
kg/ha 
kg/ha 


Min. 
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Maximum, and Minimum 


No. of 
Sites 


oo 
20 
PHS 
25 
ao 
75. 
25 
25 
ao 
715) 
25 
25 
78 
Pho) 


100 








Table 26. Frequency Distribution of Responsive 
and Unresponsive Wheat Sites per 
Classification Class (Number of Sites) 


Classification 


Agro-climatic Area 
2A 


Soil Zone 
Black 
Thin Black 
Dark Brown 


Soil Order 
Chernozemic 
Solonetzic 


Texture (0-15) 
CL 
SCL 
L 
SL 


Unresponsive 


ou -Of 


ono~ 


Responsive 
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1. Discriminant Analyses 

The objective of this series of analyses was to 
determine those site variables important for distinguishing 
responsive and unresponsive wheat sites to phosphate 
fertilizer. There appeared to be very little difference 
among the soil test procedures for purposes of separating 
responsive and unresponsive sites (Table 27). In general, 
separation was very “a. The best overall functions were 
determined using Olsen-P and ASFTL-P in separate functions 
(Table 27). Again there was very little difference between 
the two procedures for separating sites into unresponsive 
and responsive. The function using Olsen-P also included 
% O.M. and E.C., while the function using ASFTL-P had only 
% O.M. aS an additional variable important for 
discrimination. Since the function using ASFTL-P only 
required one additional variable to obtain the same degree 
of separation as that for the function using Olsen-P, it was 
much eaSier to use. Therefore, comparison of site 
classification was made using the ASFTL-P function. 

Inclusion of the site classification variables into the 
function did not improve the function correlation 
(Table 28). As a result, individual class discriminant 
analyses were determined. These functions varied as to the 
number and types of variables important for site separation 
(Table 29 and 30). Even the best phosphorus soil test 
procedure varied among the classes. For separating the 


Sites, the best functions were within the soil zone and soil 
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Table 27. Discriminant Analyses for Wheat Response to Phosphate 
Fertilizer: (1) Comparison of Soil Tests, and 
(2) Best Overall Function 


Std. 
Variables Coef. 


Unstd. 
Coef. 


1. Comparison of Soil 


(38 sites) 
ASFTL-P(0-15) 1.00 
constant 
M & A-P(0-15) 1.00 
constant 
Olsen-P(0-15) 1.00 
constant 


2. Best Overall Functions (38 sites) 


Ln Olsen-P(0-15) 1.07 


2 O.M. (0-15) 0.54 
E.C. (0-15) =0.39 
constant 
ASFTL-P(0-15) 1.00 
% O.M. (0-15) O77 
constant 


** significant at p < 0. 


0.04 
=1.09 


0.04 
PERE: 


0.08 
=2.80 


2.89 
0.34 
SPIE, 
SOc, 


0.04 


0.48 
aaa 


01 


Resp. 


Tests 


=Ono7 


=0.36 


sO; 41 


= Oia 


=O toe, 


Group Centroid 
Unresp. TDP Chi-sq. 


0.71 


0.70 


O79 


0.91 


1.00 


hae) 


Canonical 


Correl. 


0.46** 


0.46** 


0.50** 


0.56** 


0.60** 
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Table 28. Discriminant Analyses for Wheat Response to 
Phosphate Fertilizer: Quantitative and 
Site Classification Variables (38 sites) 


Std. Unstd. Group Centroid Canonical 
Variables Coef. Coef. Resp. Unresp. TDP Chi-sq. Correl. 


Agro-climatie Area 


ASFTL-P(0-15) 0.99 0.04 
SosM. (0-15) 0.74 -0.46 
Argo-climatic Area 
1 0.04 0.04 
constant SEaS7) (Wig sv 1.00 OeSs 1Syen1 0.60** 


Soil Zone 


ASFTL-P(0-15) 0.90 0.04 
mOuM cme GO= 1.5.) 0.69 0.43 
Soil Zone 
Black = Qyen Son Ole tS 
Dark Brown = 0) Sana Orne 
constant —d 0-0. 55 12102 0.34 15s 0.60** 
Soil Order 
ASFTL-P(0-15) te 110 0.05 
20M. V0—15) Neng AS} 0.49 
Soil Order 
Chernozemic Ons 0.30 
constant =5.00 -0.54 1.04 0.35 Gre 0.61** 


** Significant at p s 0.01 
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Table 29. Discriminant Analyses for Wheat Response to Phosphate 
Fertilizer for One Agro-climatic Area and Two Soil Zones 


Std. Unstd. Group Centroid Canonical 
Variables Coef. Coef. Resp. Unresp. TDP Chi-sq. Correl. 


Agro-climatic Area '1' (26 sites) 


ASFTL-P(0-15) Vee 0.05 
MOeM. (0-15) 0.82 0.65 
constant Nl) SG 0.88 E35 12:50 0.62** 


Black Soil Zone (15 sites) 


Pptn. 0.83 0.10 
Zoclay (0-15) 3.96 027s 
% sand (0-15) SSO 0.45 
Ln M & A-P(0-15) 0.93 1.95 
EG-n0—45)) 0.83 5.29 
constant -47.28 On 6er=2. 10 (Sm ow) 0.81%** 


Thin Black Soil Zone (13 sites) 


Pptn. Vos 0.36 
Zz clay, (0-15) 220 0.38 
% sand (0-15) a1 OeS7 
Ln Olsen-P(0-15) Sheath Ah 
SEO.Me. CO 15) 2.39 641 
Peco CO—i.5)) eles Si aes! 
PH (0-15) = (O/cur Oine oO) 
constant eafotb TWP ha 4.51 0.97 2036 0.99%* 


** Significant at p < 0.01 





106 


Table 30. Discriminant Analyses for Wheat Response to 
Phosphate Fertilizer for Two Soil Orders 


Std. Unstd. Group Centroid Canonical 
Variables Coef. Coef. Resp. Unresp. TDP Chi-sq. Correl. 


Chernozemic Sites (18 sites) 


Ln Olsen-P(0-15) -2.04 -6.16 


% O.M. (0-15) cease, SUNS 
Pptn. Aw, Os19 
ReGen cO= 15) litem 3294 
PH (0-15) Onde 2.29 
% sand (0-15) 1.91 Oueaiia 
mee Lays (O= 15) ie 6 0.19 
constant 110.42 e022 104 Onze 17.4 0.87** 


Solonetzic Sites (20 sites) 


% sand (0-15) 2.90 0.41 
PH (0-15) Penh Shy 17. 
t7O.M. (0O—15) hes Tae 
ReClay | 0= 1:5) 221 0.45 
M & A-P(0-15) -0.91 -0.04 
constant 13787 Wezel 0.60 1557 0.80*#* 


ee ae a a ee ee 


** significant at p < 0.01 
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order classes, while the function for agro-climatic area '1' 
had a poor ability to separate sites as indicated by the 
relatively low total discriminatory power and canonical 
correlation. However, care must be exercised when examining 
these functions because of the small sample size which may 
have resulted in a general inconsistant behavior of the site 


variables among the functions presented. 
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2. Multiple Regression Analysis 

As in the case of the barley and rapeseed sites, the 
calculation of wheat response to phosphate fertilizer meant 
that yield increase was dependent upon a calculated 
fertilizer rate. The variation of this relationship among 
all sites should be due to variation in the site properties, 
and multiple regression procedures could be used to identify 
those site variables responsible for this variation. 

It should be noted that the dependent variable used in 
these analyses was percent yield increase (see Material and 
Methods). This was done because of the very large variation 
in yield increase that could not be explained by the 
independent site variables other than the phosphate 
fertilizer rate. Percent yield increase was used in an 
attempt to remove some of the unmeasured environmental 
factors which may have influenced the variation in crop 
response. No comparison of the soil test procedures was 
necessary since only the Olsen-P proved to be significant in 
accounting for variation of percent yield increase. The best 
combination of quantitative variables as determined by a 
stepwise multiple regression analysis, and the approximate 
additional variation each explained, included: the phosphate 
fertilizer rate for optimum yield (31%), E.c.(0-15) (10%), 
and Ln Olsen-P(0-15) (20%) (Table 31). An increase in E.C. 
tended to enhance the percent yield increase, while an 
increase in Olsen-P depressed the percent yield increase. 


Altogether, this function was able to explain 61% of the 
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Table 31. Percent Yield Increase Equations for 25 Responsive 
Wheat Sites: (1) Best Combination of Quantitative 
Variables, and (2) Site Classification Using Stepwise 
Multiple Regression Anaylsis 


Overall 
b Std.Err. F R? SyeclMo sas MS 
Variables Value b Value Change Est. Value 


(1) Best Combination of Quantitative Variables 


P.O; (90% Max.Yld.) 0.46** 0.09 PASV SEN) OST 


B.C. (0-15) 56.87** 14.58 15.21 O10 
In Olsen-P(0-15) -20.14%* 6.17 10.65 0.20 
constant 48.01 7.56 10.85 


(2) Site Classification 
Agro-climatic Area 


P.O, (90% Max.Yld.) 0.39%* 0.09 21.14 Oran 
Boe 0-15) 63.04** 13.68 Bile es ORO 
Ln Olsen-P(0-15) -21.14** 5.69 US gle 0.20 
Agro-climatic Area 
1 =3376 1.70 
constant 51.86 6.95 10.86 


Soil Zone 


P.O, (90% Max.Yld.) 0.40** 0.09 2159 esi | 
Bee 15) (SV ASE ES eh 21.54 0.10 
En Olsen-P(0=-15) =20.69%* 5.71 Ass tie 0:20 
Soil Zone 


Dark Brown 5ee2 FAH 
Black -0.85 2205 
constant 48.26 6.96 Bec 
Soil Order 


P.O, (90% Max.Yld.) 0.46** 0.09 24.49 0.30 
Be Gem. O— 05) 55.58** 18.11 9.42 OSG 
Ln Olsen-P(0-15) -20.56** 7.17 B22 0.20 
Soil Order 
Chernozemic (0) 5 Sis) 2.81 
constant 50.04 Wa fhe U AUS 


** significant at p s 0.01 


R? 


0.61** 


0.69** 


0.70** 


0.61** 
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variation in percent yield increase. 

To determine if inclusion of site classification would 
improve percent yield increase prediction, effect coded 
classification variables were forced into the function 
(Table 31). The inclusion of agro-climatic area or soil zone 
accounted respectively for an additional 8% and 9% of the 
percent yield increase variation. Honea: inclusion of soil 
order variables did not improve the regression correlation. 
To determine if a significant difference existed among the 
agro-climatic areas or soil zones, an approximate multiple 
range test was used on the estimated class means (see 
Material and Methods). The results indicated no significant 
difference among the means within either agro-climatic area 
or soil zone classifications (Table 32), even though a 
relatively large percentage of the variation in percent 
yield increase was accounted for by these variables. This 
was probably due to the large variation in the estimated 


means aS indicated by the high standard errors. 





= 
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Table 32. Comparison of Mean Percent Yield 
Increase for Responsive Wheat 
Sites in Various Classes 


Classification Mean Sta. Error 


Agro-climatic Area 


1 hei SREP ot 1.96 

2A 16.32 a 2.61 

= 15.82 1.63 

Soil Zone 

Black iWon PAS! al 2.44 

Thin Black To 5aa SG 

Dark Brown 1Gsst) a 2.68 

x 15.70 7% 


Means within a classification having different 
letters are significantly different (P < 0.05) 





3. Principal Component Analysis 

The interrelationships among the measured independent 
Site variables of the responsive wheat sites were determined 
using principal compnent analysis. The sum of the four 
largest components explained about 85% of the total variance 
of the data (Table 33) 

Principal component number 1 accounted for about 34% of 
the variation. It was heavily loaded by pH, % organic 
matter, and the calculated optimum fertilizer rate and 
moderately loaded by % clay, precipitation, and Olsen-P. The 
phosphate fertilizer rate had an inverse relationship with 
PH, % Organic matter, % clay, and precipitation, and a 
direct relationship with Olsen-P. This suggested that as pH, 
% organic matter and/or precipitation increased, the optimum 
fertilizer rate decreased. This would imply a mineralization 
process or a phosphate sorption mechanism by the soil 
Organic matter, a greater importance of fertilizer 
phosphorus under arid conditions, plus a greater need for 
phosphate fertilizer by wheat as soil pH decreased. The 
direct relationship of the fertilizer phosphate requirement 
with the soil test for phosphorus (Olsen-P) is contrary to 
the definition of a soil test, and suggests that the soil 
test did not provide a measure of the available phosphorus 
in the soil. The inverse relationship between % clay and the 
phosphate fertilizer rate is again contrary to that found in 
the literature. This component could be labelled as the 


"phosphate fertilizer component". 
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Table 33. Principal Component Analysis of Responsive 
Wheat Sites: The Four Largest Eigenvalues 


Principal Component No. 1 
Eigenvalue 2.412 
(cumulative percentage) 34.5 
Eigenvectors 

pH (0-15) 0.855 
Bec. (0-15) ; = Oren 
R2O.M. oe t0=15) 0.675 
% clay (0-15) On537 
Pptn. 0.367 
Ln Olsen-P(0-15) sa ee de 
P.O, (90% Max.Yld.) = (eho 


1.741 
Bee 


=O 310.26 
Oa7 71 
0.407 

=U. 70 
0.564 
0.783 

=O Ss 


0.996 
HENS 


=0% 207 
-0.442 
0.248 
0.423 
0.485 
Ose 
Vas 


tg 
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The second principal component accounted for about 25% 
of the variation and was heavily loaded by E.C. and Olsen-P, 
and moderately loaded by % organic matter and precipitation. 
The direct relationship between organic matter and 
precipitation suggested a soil zone trend, but the 
domination of the component by E.C. and Olsen-P suggested a 
minor role of the zone trend. Both E.C. and Olsen-P are an 
indication of the ionic potential of the soil solution, E.C. 
for ionic concentation and Olsen-P for solution and adsorbed 
phosphorus. Therefore, this component was labelled the "soil 
solution component". 

Principal component number 3 accounted for about 14% of 
the variation and was loaded moderately by E.C., % clay, 
precipitation and the optimum phosphate fertilizer rate. A 
sorption relationship was indicated by this component, i.e., 
as clay content increased, the salt content of the soil 
solution (E.C.) decreased and the phosphate fertilizer rate 
needed for optimum growth increased to overcome phosphate 
sorption by the clay. Therefore, this component was labelled 
the "clay sorption component". 

Principal component number 4 accounted for about 11% of 
the variation with the important variables being E.C., 

% clay and precipitation. No explanation for the 
relationship of these three variables can be offered. 

Principal component analysis was meant for data 
reduction of large data sets and not for small data sets as 


was the case here. A number of unidentifiable or contrary 
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relationships was found which may be due to the relatively 
_ small size of the data matrix. 


- 1 ria .@ rir 
: » 


bbe . 
= Rebs 3 
re 
a eis 
ve 
a 6 _ 
nas 
a, 
; 7 pers: ° 7 fend 1 ¢ «4 ices nee sta 2 a 
a ; a - I 


_ a ; 
poae.) Orier anes .¥ on Thcehaee 


Val lL sGae 
3 _ 


aka - 
_ he =p 
” = 






- ! 


ih. re th 


oe 


} oy aginaned 2siae | 


. «| fa 
sf iv. osed Ding 





116 


4. Summary 

The results of the discriminant analyses of the wheat 
Sites indicated that there was very little difference among 
the soil test procedures for separation of responsive and 
unresponSive sites. In addition to the soil test for 
phosphorus, organic matter content of soils was an important 
discriminating variable, as was E.C., depending upon the 
soil test procedure used in the analysis. A high soil test 
for phosphorus and/or organic matter tended to allocate a 
Site into the unresponsive group, whereas, a high E.C. 
tended to allocate a site into the responsive group. 
Individual classification class discriminant analyses 
resulted in more highly correlated functions than the 
function using the effect coded variables. This was due to 
the difference in the list of discriminant variables and 
their importance and behavior among the classes. 

Multiple regression analysis of the responsive wheat 
Sites indicated that the Olsen-P was the only soil test 
procedure able to significantly account for variation in the 
percent yield increase of wheat to phosphate fertilizer. The 
only other measured quantitative variable which was 
Significant was E.C.. Inclusion of either agro-climatic area 
or soil zone into the analysis increased the correlation 
coefficients of the equations. Soil order did not have the 
same effect. However, even with the improved correlation, 
there was no significant difference among the classification 


class means. 


Principal component analysis of the responsive wheat 
Sites revealed some recognizable relationships among site 
properties, and with the calculated optimum fertilizer rate. 
The site variables measured can be reduced to three 
components representing (1) phosphate fertilizer rate, 

(2) soil solution, and (3) clay adsorption. The most 
noteworthy relationships were the inverse relation between 
phosphate fertilizer rate and each of pH, % organic matter, 
and precipitation. Contradictory results were also noted, 


possibly being due to small sample size. 
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D. Sources of Variation 


Several potential sources of variation exist in the 


Study of crop response to fertilizer. These are discussed 


with reference to the present study. 


1. 


The general field designs used in this study varied 
among the cooperators and were quite unique when 
compared to those found in the literature. As a result 
of careful examination, the procedure outlined in the 
Materials and Methods appeared to be the only route open 
to satisfy the objectives. Some of the problems 
encountered included: 

(a) In the original design of the project, a basic 
assumption was made concerning the relationship between 
cropping history of a site and nitrogen levels in the 
soil. It was assumed that fallowed sites would contain 
more plant available nitrogen than previously cropped 
sites, and as a result, blanket rates of nitrogen 
fertilizer differed depending on cropping history. Sites 
cropped the previous year received more nitrogen 
fertilizer than sites fallowed the previous year. This 
was compounded by use of different blanket nitrogen 
fertilizer rates among the cooperators. Therefore, 
cropping history aS a site variable became related to 
nitrogen fertilizer rates. Separation of these variables 
was not possible and a combined variable was used. 
Analysis of covariance using an effect coding indicated 


no significant effect of this combined variable on yield 
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response to phosphate fertilizer for all three crops. 

(b)The plot design, number of treatments and 
replication varied not only among cooperators, but also 
from year to year for a particular cooperator. 

(c) In a number of cases, the highest phosphate 
fertilizer treatment was not great enough to establish a 
true maximum yield for a site. For these sites, 
calculation of 90% maximum yield was based on the 
highest fertilizer rate and not on an extrapolation of 
the response function. 

(d) The design of most of the experimental sites 
provided no information on possible interactions of 
plant nutrients. 

The type of equation used for calculating the response 
function for each site was chosen based on visual 
examination of the plotted yield data for each site, 
ease of calculation, and ease of mathematical 
manipulation. Only one type of function (second order 
polynomial) was used, and in some cases the equation was 
forced to fit the data such that the fit was poor. Poor 
fits were due primarily to insufficient number of 
treatments to adequately define the response curve and 
to possible lack of uniformity within the plot site. 

The lack of precipitation data for some sites forced the 
use of estimated values based on the nearest 
meteorological station. These estimated values may not 


have reflected the actual rainfall for the plot sites in 
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question. The influence of the distribution of 
precipitation over the growing season and the initial 
Soil moisture conditions were not determined due to a 
lack of data. 

Incomplete data for the M & A-P and Olsen-P procedures 
for some sites, due to the loss of original soil 
samples, forced-comparison of soil test procedures being 
made on a reduced data set. 

Site soil analyses were based on a composite soil sample 
for the site and not on individual treatments and/or 
replicates. This resulted in the assumption that the 
soil samples were representative of the plot site, and 
that the plot was uniform-in terms of soil properties. 
One or more of the variables investigated may have been 
truly unrelated to crop response but remained correlated 
due to chance. In the present study, attempts were made 
to give plausible explanations for significant 
correlations between independent and dependent 
variables. Definite causal relationships were however, 
difficult to determine. The validity of certain factors 
should be checked by analysis of new data. 

The correlation between dependent and independent 
variables may have been nonlinear. This source of error 
was minimized in the present study by making 
scattergrams of dependent versus independent variables 
as described in the Material and Methods chapter, and 


applying the appropriate transformation to the 
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independent variable to approximate a linear 
relationship. 

The number of sites, or sample population, in many 
analyses was quite small, thus possibly influencing the 
reliability of the results. 

Multicollinearity exists when any independent variable 
is correlated with another independent variable or with 
a linear combination of other independent variables. 
Multicollinearity is common and even inevitable in much 
of the data in soil science. Correlation among the 
independent variables causes three main problems: 

(i) the standard errors of the regression coefficients 
are increased, (ii) as the extreme case of (i) is 
approached, computational difficulties arise, (iii) the 
omission of variables may result in biased estimators 
for the regression parameters of the remaining variables 
if the missing variables are correlated with those 
remaining. In general, there is little that can be done 
about multicollinearity except to take a larger sample, 
preferably in a way that decreases multicollinearity 
(Wesolowsky, 1976). 

The basic difficulty with data derived from a series of 
fertilizer experiments is that the sources of variation 
differ between and within experiments. If these are not 
recognized, it is easy to obtain invalid tests of 
significance by using inappropriate eStimates of error 


variance. The source of error affecting between site 
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relationships are primarily due to factors varying in an 
unknown or unidentified manner throughout the region. 
Since these error effects vary with both location and 
time, the error variance cannot be estimated by 
replication. Rather it must be estimated indirectly, as 
by the residual mean square of an appropriate regression 
analysis of variance.(Colwell, 1978). 

The selection of 90% maximum yield for a site as the 
optimum yield may not have been valid. This selection 
was based on the examination of a general response curve 
which indicated that potential yield values near the 
maximum yield for the site changed very little, 
depending on the partial regression coefficients for the 
Site, while the fertilizer rate could change quite 
dramatically. To provide a standard procedure, 90% of 
the maximum was arbitrarily selected as a yield that 
approximated an economic optimum as well as a biological 
optimum. 

In this study, a simple separation of sites into 2 
categories, responsive and unresponsive, was used. This 
Separation did not take into account different levels of 
responsiveness (high, medium, and low). 

With a few exceptions, analyses using effect coded 
variables (discriminant and multiple regression) were 
unable to indicate differences among the classes of a 
classification. This could be due to the assumption that 


the slope of the regression lines are equal among the 
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classes when effect coding is used. If this assumption 
was not correct, then a weighted coding may have been 
necessary. 

The results of the statistical analyses in this study 


were not verified with data external to this study. 
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V. SUMMARY AND CONCLUSIONS 

The aim of this study was to determine the influence of 
various soil properties and site classifications on the crop 
response to phosphate fertilizer in Alberta. As noted, 
rather poor correlations exist between the soil test for 
phosphorus and percent yield from combined field experiments 
in Alberta (R? = 0.53).* Good correlations between yield and 
the nitrogen and phosphate fertilizer rates were found for 
individual site-years by Heapy (1971) when soil tests for 
available nitrogen and phosphorus were included in the 
response function. However, when the individual site-years 
were combined, correlations were poor. Greenhouse studies 
have shown high correlations between yield response and soil 
test phosphorus (Robertson, 1962). Significant differences 
among cereal crops with respect to crop response to 
phosphate fertilizer was noted by Robertson et a] (1968). In 
addition, numerous studies have noted the influence of 
various soil properties on the chemical reactions and 
availability of phosphate fertilizer within soils. Since 
(i) the correlations from greenhouse studies have been 
considerably better than those for field studies, and 
(ii) the correlations from individual site-year field 
experiments were better than those for which site-years were 
combined, there would appear to be an influence of the site 
environment (soil and climate) on the crop response to 
phosphate fertilizer. Therefore, rather than attempt to 


‘Personal communication with Dr. J. A. Robertson. 
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develop and/or test new soil test procedures, the influence 
of soil and climatic properties on the yield response to 
phosphate fertilizer was examined. 

The analyses of crop response in this study were broken 
down into two fundamental questions based on the purpose of 
a soil test: (1) Will a crop respond to phosphate fertilizer 
application at a particular site? and (2) If the answer to 
(1) is yes, then what is the magnitude of the response? To 
answer these questions, this study attempted to determine 
the influence of various site properties uSing two separate 
but related analyses: (1) discriminant analysis to separate 
Sites into responsive and unresponsive categories, and 
(2) multiple regression analysis to account for the 
variation in yield increase of the responsive sites. In 
addition, principal component analysis was used to determine 
the interrelationships among site variables of the 
responsive sites. The results of these statistical 
techniques were used to try to understand the variation in 
site response to phosphate fertilizer application. 

Results of the analyses of the barley sites indicated 
that the most important site property influencing both site 
response and yield increase to phosphate fertilizer was the 
soil test (ASFTL-P). Other site variables that were 
important for site separation included clay and Caco; 
content of the soil, and growing season precipitation while, 
Soil pH, growing season precipitation, and organic matter 


content of soils significantly accounted for variation in 
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yield increase of the responsive sites. Site classification 
improved the correlation coefficients of both the 
discriminant and multiple regression analyses, indicating 
Significant differences in crop response to phosphate 
fertilizer among some classes, particularly those sites in 
the gray soil zone or members of the Luvisolic soil order. 
Principal component analysis indicated that the required 
phosphate fertilizer rate for "optimum" yield response was 
inversely related to ASFTL-P, soil pH, and the organic 
matter content of soils. Thus for barley, the phosphate 
fertilizer rates should be reduced as ASFTL-P, pH, and/or 
% organic matter increase. 

Results of the analyses of the rapeseed sites suggested 
that the crop response to phosphate fertilizer was 
influenced by site properties different from those for the 
barley sites. Again, the most important site parameter 
influencing crop response to phosphate fertilizer was the 
soil test for phosphorus (ASFTL-P). The other site variables 
that significantly influenced site separation were E.C. and 
clay content of the soil, while CaCO; content of the soil 
was the only other site parameter that accounted for 
variation in yield increase of the responsive sites. Site 
classification was important for site Separation but not for 
explaining variation in yield increase. Principal component 
analysis of the responsive sites indicated trends similar to 
those found for the responsive barley sites. The required 


phosphate fertilizer rate for "optimum" yield was inversely 
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related to ASFTL-P, soil pH, and soil organic matter 
content, but also, to growing season precipitation. 
Therefore, phosphate fertilizer rates for "optimum" yield 
response of rapeseed should be reduced as ASFTL-P, soil pH, 
Organic matter content of soils and/or growing season 
precipitation increase. 

mreereddaleiys the locations of the wheat sites 
differed considerably from those of either barley or 
rapeseed, making crop comparisons almost impossible. The 
results of the analyses of the wheat sites indicate that a 
soil test for phosphorus was the most important site 
variable influencing crop response. The other site 
properties influencing site separation were organic matter 
content of soils and soil E.C., while only the additions of 
soil E.C. explained variation in percent yield increase of 
the responsive sites. Site classification had a variable 
influence. For site separation, site classfication appeared 
to be important, especially for individual class functions 
and for determining the best soil test procedure for 
phosphorus. For the variation in percent yield increase of 
the responsive sites, inclusion of site classification 
resulted in a large improvement in the correlation 
coefficient of the percent yield increase equation, but 
there was no significant difference among the class means 
when compared. Principal component analysis showed a number 
of the same trends as observed for the responsive barley and 


rapeseed sites, that is, the required phosphate fertilizer 










ve. thy low SOs oF Hema 
cas] a ara sud ,2nSa7T 
a 


aN " ‘> Stai(oC2zcne .3)92073m 


al Vetteiat 10 Fata 
24982 5 ’ : 
es 
1 <7 
tld 
& ar. >a 
4 
e 


jars e~ eatin. @ sO St 
eo Aas 


. a . - e 
‘4 7". San ox 139 congels 


>: seni sen is Ou eet 


ip 


4; ans eect 19% k 
tan 







a6 ob4p7sa. Sisinv eed nz eErrer a “von & 





, il oA An : CARE hn! 


128 


rate for "optimum" yield was inversely related to soil pH, 
soil organic matter, and growing season precipitation. 
However, the relationship between the soil test for 
phosphorus (Olsen-P) and phosphate fertilizer rate was 
contrary to the barley and rapeseed results, and to the 
commonly expected relationship. This contradiction could be 
due to either the small number of sites or to the inability 
of the soil test to provide an indication of the available 
phosphorus status for these sites, especially those classed 
as Solonetzic. 

In conclusion, the soil test for phosphorus does not, 
by itself, provide a satisfactory measurement for separation 
of responsive and unresponsive sites, nor for the variation 
in yield increase of the responsive sites. The inclusion, of 
other site properties did improve the correlation 
coefficients, but their contribution to the overall function 
R? was generally smaller than that of the soil test. Site 
classification using either effect coding or analysis of 
individual classes did improve on the correlations, with the 
individual analyses having the better results for site 
separation. It would be preferred that the coded function 
was more successful because of the difficulty in using 
individual class functions. The results of this study cannot 
be considered as conclusive and they need to be verified 
with external data. They do suggest that the phosphate 
fertilizer rate for "optimum" yield response should be 


reduced as ASFTL-P, soil pH, organic matter content, and/or 
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growing season precipitation increase. The properties 
identified as influencing crop response to phosphate 
fertilizer can be used in further modelling designed to 
derive more specific calibration curves for predicting 
phosphate fertilizer requirements. Additional work is needed 
to determine the influence of meteorological variation, 
cropping history, soil and fertilizer nitrogen levels and 
micronutrient levels on the crop response to phosphate 
fertilizer. Alternative approaches to measuring the 
phosphorus fertility status of soils may also have to be 


examined. 
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APPENDICES 


Site Identification Code 
Code Cooperator 


Agriculture Canada, Beaverlodge 

Alberta Agriculture, Edmonton 

Agriculture Canada, Lacombe 

Western Co-operative Fertilizer Ltd., Calgary 
Agriculture Canada, Lethbridge 

Dr. J.A. Robertson, University of Alberta 
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APPENDIX A 


Experimental Year, Crop Variety, Cropping History, and 


Legal Location of Experimental Sites 


NA indicates data were Not Available 





Table A-1. 


Site 


BO 1 
B02 
B03 


B04 
B05 
BO6 


BO7 
B08 
BOS 


Year 


nor 
1O71 
1371 


ia 
1971 
1971 


oh 
tort 
age 
aoe 
i i 
O72 


ee 


372 
122 
1972 


Pde 
1272 


T9372 
$372 
'o73 


1373 
$273 
1973 
1973 
1973 
i 23 


1973 
je WES 
RE 
ake 
hotdes 
1974 
1974 
1974 
1974 
1974 
1974 
1974 
1974 


Barley Sites 


Variety 


Galt 
Galt 
Galt 


Galt 
Galt 
Galt 


Galt 
Galt 
Galt 
Galt 
Galt 
Galt 


Galt 


Galt 
Galt 
Galt 


Galt 
Galt 


Galt 
Galt 
Galt 


Galt 
Galt 
Galt 
Galt 
Galt 
Galt 


Galt 
Galt 
Galt 
Galt 
Galt 
Galt 
Galt 
Galt 
Galt 
Galt 
Galt 
Galt 
Galt 


Past Cropping History 


1970-Fallow 
1970-Fallow 
1970-Oats 
1969-Wheat 
1970-Cropped 
1970-Barley 
1970-Barley 
1969-Barley 
1970-Wheat 
1970-Fallow 
1970-Fallow 
1970-Fallow 
1971-Fallow 
1971-Fallow 
1970-Barley 
1971-Fallowed Fescue 
1970-Fescue 
1971-Fallow 
1971-Fallow 
1971-Barley 
1970-Rapeseed 
1971-Rapeseed 
1971-Barley 
1970-Fallow 
1971-Rapeseed 
1971-Barley 
1972-Barley 
1971-Barley 
1972-Fallow 
1972-Fallow 
1972-Fallow 
1972-Fallow 
1972-Fallow 
1972-Wheat 
1971-Fallow 
1972-Partial Fallow 
1972-Cropped 
1972-Cropped 
1972-Cropped 
1972-Cropped 
1973-Fallow 
1973-Fallow 
1973=Fallow 
1973-Fallow 
1973-Cropped 
1973-Cropped 
1973-Cropped 
1973-Cropped 


Legal Location 


CSD -11-265073= 10-W6 


NW 
NE 


NE 
LSD 


13-081-02-W6 
35-108-12-W5 


{hes By aS tah 


06-26-073-10-W6 
05=376-07/1-09-Wwe6 


O2=0:78-20-WS 
PAs OLEaES-Wo 
6-072 ia WG 
23-063-01-W6 
09-0 70=10-W6 
16-072-11-W6 


23-078-10-W6 
Ot=709-12=W5 
Ot=108=13-W5 
07-072-07-=W6 


S2>012R biEwG 
+2-078-19-W6 


24 110=19-W5 
O2=h 08st 3-W5 


08-07-072-07=W6 


16-36-083-24=W5 
02—)g-07= 15=nS 


09-097 we 
02-108=33-W5 
05-109-07-W5 


03-22-080-02-W6 


09-26-073-10-W6 


01-108-13-W5 
16-4170 9-Wo 
16—107—1o-Wes 
24-107=13=N5 
14-108-13-W5 
04-107-12-W5 
$6—107= 1 5-Wo 
01-108-13-W5 
F710 7-1 3"0o 
Oi=108=13-45 
08-108-17-W5 
16-1.10-19-WS 
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Table A-1. 


Site Year 


E01 


he a 


1971 
1375 


Eh a 
197d 
igs. 
STs 
19714 
19 9H 
hats 


‘Siz 
O72 


Lo 2 
1972 
1972 
Lore 
ty 2 
LRH 
epee 
no73 
Bee 


(Ae 
tod 


ios 


1973 
y Rey ee 


ues 


Barley Sites (cont.) 


Variety 
Galt 


Galt 
Galt 


Galt 
Galt 
Galt 
Galt 
Galt 
Galt 
Galt 


Galt 
Galt 


Galt 
Galt 
Galt 
Galt 
Galt 
Galt 
Galt 
Galt 
Galt 


Galt 
Galt 


Galt 


Galt 
Galt 


Galt 


Past Cropping History 


1970-Oats and Barley 
1969-Sod-Breaking 
1970-Barley 
1970-Fallow 
1969-Oats 
1970-Fallow 
1970-Barley 
1970-Barley and Oats 
1970-Barley 
1970-Barley 
1970-Fallow 
1970-Fallow 
1969-Sweet-Clover 
1971-Barley 
1971-Fallow 
1970-Hay-Sod 
1971-Barley 
1971-Oats and Barley 
1971-Barley 
1971-Barley 
1971-Wheat 
1972-Barley 
1972-Barley 
1972-Barley 
1972-Wheat 
1971-Rapeseed 
1972-Barley 
1972-Fallow 
1971-Rapeseed 
1972-Wheat 
1971-Wheat 
1972-Barley 
1972-Barley 
1971-Barley 
1972-Barley 
1971-Rapeseed 


Legal Location 


SW 


18-055-—23-We 


259-56" 27-We 
24-055-24-W4 


18-062-26-W4 
25-032-04-W5 
25-046-27-W4 
16-058-25-Ww4 
29=053=U1-wo 
16-063-26-W4 
27-046-25-W4 


04-049-27-W4 
01-049-22-W4 


05-049-19-W4 
04-049-19-W4 
30=032-02-Ws 
26-033-01-Wo 
O3-U3s3-47-We 
25-U59=20-44 
ty -O59=21-We 
13-054-24-Ww4 
O520o9=1 3-W4 


07-057-24-W4 
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31-058-21-W4 


20-060-17-W4 


20 Sh =e iawe 
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24-059-17-W4 
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Table A-1. 
Site Year 
LO 1 1971 
LO3 1971 
LO4 1971 
LO5 1971 
LO6 1971 
LO8 1971 
H10 606 1972 
BVT 1972 
L12 1972 
mis 1972 
m1e° 1972 
trios 1972 


Barley Sites (cont.) 


Variety 


Galt 
Galt 
Galt 
Galt 
Galt 
Galt 
Galt 
Galt 
Galt 
Galt 
Galt 
Galt 


Past Cropping History 


1970-Barley 
1970-Barley 
1970-Rapeseed 
1970-Wheat 
1970-Fallow 
1970-Barley 

197i- Cats 
1971-Sweet Clover 
1971-Rapeseed 

197 ?-Oats 
1971-Sweet Clover 
1971-Sweet Clover 


Legal Location 


02-055-08-W4 
08-054-11-W4 
U6-03-26-Ww4 
12-054-09-W4 
08-054-11-W4 
35-052-08-wW4 
36-053- biew4 
14-053-08-Ww4 
34-054-07-Ww4 
33-053-11-W4 
327054-13-W4 
36-053-11-W4 
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Table A-1. Barley Sites (cont.) 


Site Year Variety Past Cropping History Legal Location 


WO1 1971 Conquest 1970-Fallow SW 14-026-23-W4 
WO2 1971 Betzes 1970-Fallow SW 29-027-19-W4 
WO3 1971 Conquest 1970-Fallow NW 07-024-26-W4 
WO4 1971 Betzes 1970-Fallow NE 20-023-23-W4 
WO5 1971 Betzes 1970-Fallow NE 06-031-21-W4 
WO6 1971 Conquest 1970-Fallow SW b-022-25-w4 
WO7 1972 Betzes 1971-Fallow NE 33-023-28-W4 
WO8 1972 Betzes 1971-Fallow SW 02-023-27-W4 
wWOS 1972 Galt 1971-Barley NE 02-030-01-W5 
1970-Cereal 
Wid 1972 Betzes 1971-Barley SW 34-031-27-W4 
1970-Cereal 
Miz. «1972 Betzes 1971-Fallow SW 09-024-21-W4 
Wis). 1972 Betzes 1971-Fallow SW 26-01:1-27-—w4 
W14 1972 Galt 1971-Fallow SE 18-027-28-W4 
Wino 1972 Betzes 1971-Fallow SE 04-027-21-W4 
W116 1972 Betzes 1971-Fallow NE 24-025-23-W4 
why 1972 Betzes 1971-Cereal SW 17-032-07-Wd 
W18 1972 Betzes 1971-Fallow NW 14-017-02-W5 
W1i9° 1972 Galt 1971-Cereal NE 14-03 1-C2-Ww5 
W220) 1972 Betzes 1971-Grazed Crop Cover SE 16-018-29-We 
1970-Cereal 
W22 1972 Betzes 1971-Fallow NE 26-024-27-W4 
2 PS a Be WAP Betzes 1972-Fallow SW 26-011-27-W4 
W24 1973 Betzes 1972-Fallow SW 01-024-28-WwW4 
wW25 1973 NA 1972-Fallow SW 07-032-23-W4 
W26 1973 NA 1972-Fallow SE 15-033-25-W4 
Wer. t1h973 Betzes 1972-Fallow SW 05-028-22-W4 
W28 1973 Galt 1972-Fallow SW 18-027-21-W4 
Wed 31973 Betzes 1972-Barley NW 26-024-27-W4 
Woleee toro Betzes 1972-Rapeseed NW 09-026-23-W4 
W34 1973 Betzes 1972-Wheat SW 26-011-27-W4 
1971-Fallow 
W36) 1973 Galt 1972-Barley SE 30-029-25-W4 
1971-Barley 
W37 1973 Galt 1972-Barley SE 36-029-29-W4 
1971-Cropped 
W38 1973 Galt 1972-Barley NE 17-03 2-00-W5 
1971-Barley 
W41 1973 Betzes 1973-Barley SE 13-032-24> We 
W42 1974 Galt 1973-Wheat SW 15-033-25-W4 
W43 1974 Galt 1973-Oats and Barley SE 23-029-01-W5 
W44 1974 Galt 1973-Cropped NW 12-034-01-W5 
W46 1974 Galt 1973-Barley NW 22-038-28-W4 


W47 1974 Galt 1973-Barley NE 4$1-038-01-W5 
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Table A-1. 
Site Year 
Tot 1973 
TO02 1973 
mO3 1973 
TO7 1974 
TO8 1974 
TO9 1974 
10° ~1975 
mi2 1975 


Barley Sites (cont.) 


Variety 


Galt 
Galt 
Galt 
Galt 
Galt 
Galt 
Galt 
Galt 


Past Cropping History 


1972-Rapeseed 
1972-Fallow 
1972-Barley 
1973-Fallow 
1973-Wheat 
1973-Fallow 
1974-Fallow 
1974-Fallow 


Legal Location 


06-021-23-W4 
14-022-26-W4 
02-023-28-W4 
06-021-23-W4 
02-023-28-W4 
14-022-26-W4 
06-021-23-W4 
14-022-26-W4 
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Table A-2. Rapeseed Sites 
Site Year Variety Past Cropping History Legal Location 
B41 1971 Span 1970-Fallow SE 23-085-0.1-W6 
Ba2. 1971 Span 1970-Oats NE 35-108-s2-W5 
1969-Wheat 
B43 1971 Span 1970-Cropped NE 12-24 0=13-W5 
B44 1971 Span 1970-Barley LSD 06-26-073-10-W6 
B45 1971 Span 1970-Barley SW 34-071-09-W6 
1969-Barley 
B46 1971 Span 1970-Wheat SW 02-09 8=20-W5S 
B47 1971 Span 1970-Fallow USD 11=26-0/3-10-We 
B48 1971 Span 1970-Fallow SE 17-078-19-W6 
Bao (1971 Span 1970-Fallow NW f3-US9-02=Wo 
13/1 a he Wy ps Span 197 1-Rapeseed NW Po Font INE a WS eal he 
bon -1972 Span 1971-Barley NE 02-7056-13-W5 
Boe 1972 Span 197 1-Paliow NE 09-0970-10-Wo 
pos i972 Span 1971-Fallow NW VI=708-—1S=W5 
B54 1972 Span 1971-Fallow 23-078-40-=W6 
Boo 1972 Span 1971-Barley SE 17-U 78-3 9-WG6 


1970-Fallow 
1971-Barley 
1970-Rapeseed 


Bao 4972 Span C2 =072=09 =WG 


Bo? 661972 Span 1971-Fallow NW fo=072=11-WsS 
BSS 1972 Span 1971-Rapeseed SW 32-072-11-W6 
1970-Volunteer Barley 
Bode to TS Span 1972-Barley LSD 08-07-072-07-W6 

1971-Barley 
HoOU;) 1973 Span 1972-Barley LSD l6=36-U83—24—-W5 
Bo 1973 Span 1972-Fallow ESD, 02-=49—108—15-Ws 
Hoe, 21973 Span 1972-Fallow NW O90 dea ys re 
Hos 1973 Span 1972-Fallow NE 02-208= 13-5 
B64 1973 Span 1972-Fallow NW 05-109-07-W5 
Boo 1373 Span 1972-Wheat LSD 03=22-0s0-=02-NG 


1971-Fallow 
1972-Partial Fallow LSD 09-26-073-10-W6 
1971-Fescue 


B66." 1973 Span 


Boi 1973 Span 1972-Cropped NW 01-108 =16=N5 
B68 1973 Span 1972-Cropped NW [eter SALE Peet Baw Pah be 
BOSS 97S Span 1972-Cropped SW 16-107—15-=W5 
Bye 1973 Span 1972-Cropped SW 24-—107-135-W5 


B71 1974 Span 1973-Fallow NE 14-108=13-W5 
B72 1974 Span 1973-Fallow SW 04-107-12-W5 
B73 1974 Span 1973-Fallow NW 16-10 7=15-Wo 
B74 1974 Span 1973-Fallow NE 01-108-13-W5 
Bio? 21974 Span 1973-Cropped SE ii ooh 9 tO Foaety He half 
B76 1974 Span 1973-Cropped NW W106 13-05 
BIT 3974 Span 1973-Cropped NW 08-108-17-W5 
B78 1974 Span 1973-Cropped NW 16-110-19-W6 
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Table A-2. Rapeseed Sites (cont.) 


Site Year Variety Past Cropping History Legal Location 


E37 1971 Span 1970-Barley NW 18-049-27-W4 
Ese. 1971 Span 1970-Barley NW 16-058-25-W4 
moa 1971 Span 1970-Barley SE 29-033-O07-W5 
E40 1971 Span 1970-Barley and Oats SE 24-046-27-W4 
E41 1971 Span 1970-Fallow NE 09-049-26-W4 
E42 1971 Span 1970-Barley SW 25-032-04-W5 
E44 1971 Span 1970-Fallow SW 24-055-24-W4 
E45 1971 Span 1970-Oat and Barley SW 18-055-23-W4 
1969-Sod-Breaking 
E46 1972 Span 1971-Barley NW 04-049-27-W4 
E47 1972 Span 1971-Fallow NW 01-049-22-W4 
1970-Hay (Sod) 
E48 1972 Span 1971-Barley SE 05-049-19-w4 
E49 1972 Span 1971-Oats and Barley NE 04-049-19-W4 
TASS a he Bi Span 1971-Cereal SE S0-032-02-W5 
E54 1972 Span 1971-Cereal NW 26-043-0)—W5 
Beto 7 2 Span 1971-Wheat SE 03-033=27-W4 
Ba? 1973 Span 1972-Fallow NE silasU50—40 WS 
moo (933 Span 1972-Barley NW 13-054-24-Ww4 
E60 1973 Span 1972-Barley SE 17-059-21-Ww4 
E61 1973 Span 1972-Barley SW 09=059—168—w4 
1971-Barley 
EG2 973 Span 1972-Barley NE 07-057-24-WwW4 
1971-Rapeseed 
Boe. 11973 Span 1972-Barley NE 24-059-17-W4 
1971-Rapeseed 
E64 1973 Span 1972-Wheat SE 20-060-17-W4 
1971-Wheat 
BPouneiors Span 1972-Barley SW 25-059-20-W4 
E66 1973 Span 1972-Barley SW 21-058-24-Wa 
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Table A-2. Rapeseed Sites (cont.) 


Site Year Variety Past Cropping History Legal Location 


L48 1971 Span 1970-Rapeseed SE 06-037-28-W4 
L49 1971 Span 1970-Barley SE 02-055-08-W4 
L50 1971 Span 1970-Wheat NW 12-054-09-Ww4 
L51 1972 Span 1971-Sweet Clover SE 39=053—4)1-We 
Boe 1972 Span 1971-Oats SE 36-053-71-W4 
Bos) 1972 Span 1971-Rapeseed NE 34-054-07-W4 


B54 1972 Span 1971-Sweet Clover SE 32-054-14-W4 
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Table A-2. Rapeseed Sites (cont.) 


Site Year Variety Past Cropping History Legal Location 


W49 1972 Echo 1971-Fallow NW 09-026-23-W4 
W50 1973 Span 1972-Barley NW 26-024-27-W4 
W52 1973 Span 1972-Barley SE 18-033-23-W4 
W5s 1973 Span 1972-Barley SE 36-029-01-W5 


W54 1973 Span 1972-Barley NE 117-032-095 —Wo 
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Table A-2. Rapeseed Sites (cont.) 


Site Year Variety Past Cropping History Legal Location 


T19 1971 Span 1970-Fallow SE 34-018-24-W4 
veo. 1971 Span 1970-Fallow SW 33-016-27-W4 
Tep 197.1 Span 1970-Fallow SE 33-005-=27-Wwe 
wee 19%) Span 1970-Fallow NE 27-002-14-W4 
Bin 1972 Span 1971-Fallow NW 26-002-14-W4 
mao 1972 Span 1971-Fallow SW 337005-27-="4 
mis 1973 Torch 1972-Fallow SE 14-022-26-W4 
T14 1973 Span 1972-Fallow NE 05-008-01-W5 
iS 2 1973 Span 1972-Fallow SW 33-005-27-W4 
eis. i973 Torch- 197Z2-Barley NE 02-023-28-W4 
T26 1974 Span 1973-Fallow SW 06-021-23-W4 
T27 1974 Span 1973-Fallow NE 28-002-14-W4 
T28 1974 Span 1973-Fallow SW 33-005-27-W4 
T29 1974 Span 1973-Wheat SE 02-023-28-W4 
T30 1974 Span 1973-Fallow NW 14-022-26-W4 
way 91975 Span 1974-Fallow SE 14-022-26-W4 
mse 1975 Span 1974-Fallow SW 06-021-23-W4 
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Table A-3. Wheat Sites 
Site Year Variety Past Cropping History Legal Location 
J01l 1969 Thatcher 1968-Fallow NW 03-046-17-W4 
J02 1969 Thatcher 1968-Fallow SW 25-032-17-Wa 
J03 1969 Thatcher 1968-Fallow NW 36-038-14-W4 
J04 1969 Thatcher 1968-Fallow NW 21-039-28-W4 
J06 1970 Thatcher 1969-Fallow SE 15-039-19-W4 
J07 1970 Thatcher 1969-Fallow NW 03-046-17-W4 
J08 1970 Thatcher 1969-Fallow SW 25-032-17-W4 
J09 1970 Thatcher 1969-Fallow NW 36-038-14-W4 
J10 1970 Thatcher 1969-Fallow NW 21-039-18-W4 
J11 1971 Thatcher 1970-Fallow NE 07~031-17-W4 
J1i2 1971 Thatcher 1970-Fallow NW 21-039-18-W4 
ois 1971 SThatcher 1970-Fallow NW 12-032-18-Ww4 
Ji14 1971 Thatcher 1970-Fallow NW 36-038-14-W4 
O15 1974) "Thatcher 1970-Fallow SW 25-032-17-W4 
J16 1971 Thatcher 1970-Fallow SE 18-039-18-W4 
tive 19/4 SThatcher -1970-Falblow NW 15-047-17-W4 
J18 1972 Thatcher 1971-Fallow W 29-050-19-W4 
uae 1972 "Thatcher 1971-Fallow SE 31-050-19-W4 
o2u.) 1972 “Thatcher 1971-Cropped NE 36-050-20-Ww4 
J22: 1972 Thatcher 1971-Greenfeed NW 21-039-18-W4 
1970-Cropped 
eo) 19/2 SThatcher 197 1-Fallow SE 15-039-19-w4 
O24 1972 “Thatcher 1971-Fallow NW 36-038-14-W4 
we) 19/72 eThatcher 2197 1-Cropped SW 13-049-17-W4 
J26 1972 Thatcher 1971-Cropped SE 25-050 Li -W4 
1970-Fallow 
J27 1972 Thatcher 1971-Cropped NE 09-050-17-W4 
1970-Cropped 
J28 1972 Thatcher 1971-Fallow NW 24-049-17-W4 
J29 1973 Thatcher 1972-Cropped SE 18-049-17-w4 
J30 1973 Thatcher 1972-Fallow NW 09-050-19-Ww4 
ese Gls Thatcher 1972-Fal low SE 29-050-19-W4 
J32 1973 Thatcher 1972-Fallow SE 31-050-19-W4 
uJ33. 1973 Thatcher 1972-Fallow SE 18-039-18-W4 
J34 1973 Thatcher 1972-Barley NW 03-046-17-W4 
1971-Fallow 
J35 1973 Thatcher 1972-Fallow NW 36-038-14-W4 
J36 1973 Thatcher 1972-Greenfeed SE 20-029-18-W4 
J37 1973 Thatcher 1972-Fallow SW 25-050-17-W4 
J38 1973 Thatcher 1972-Cropped NW 09-050-17-w4 
J39 1973 Thatcher 1972-Fallow SE 15-050-19-wé 
J40 1973 Thatcher 1972-Fallow SE 26-049-17-W4 
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APPENDIX B 


Classification of Experimental Sites 


List of Classification Abbreviations 











G. Gray 
D.G. Dark Gray 
BL. Black 
TBL. Thin Black 
rat. Dark Brown 
Br Brown 


Parent Material 


Lac Lacustrine 
cfs Come aura ab Lacustro Fill 
ryid ib idl Ae 

Fl Fluvial 

Aeo Aeolian 

Resid Residual 

SL Sandy Loam 

S Sand 

Sorted Till Sorted Till 


Soil Classification 


abbreviations follow Canadian System 
of Soil Classifcaton (1978) 
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Table B-1. Barley Sites 


Soil Agro-climatic Parent Sav 
Site zone Area Material Classification 
B01 D.Gi 2H Pacer ll SZ.DG 
BO2 BG. 2H PayTill D.GL 
B03 Gia 3Ha F1l/Aeo O.GL 
BO4 G. 3Ha Lac SZ.GL 
BO5 D.G. 2H Cacwril!] SZ.DG 
BO6 D.G 2H oie 111 D.GL 
BO7 G. 2H Lac SZ.GL 
B08 G. 3H Fl 0.GL 
BO9 G. 3H Race Til) SZ.GL 
B10 5 3H Fl O.GL 
B11 D.G 2H pacer SZ.DG 
B12 GC. 3H Lac Till SZ.GL 
B13 Gy 3H Lac Till S$Z.GL 
B14 G. 3Ha F1/Aeo O.GL 
B15 D.G 3Ha F1/Aeo D.GL 
B16 D.G 2H bee. Till SZ.DG 
B17 G. 3H Pacer. | SZ.GL 
Bi8 G. 3H Fl O.GL 
B19 G. 3Ha Lac SZ.GL 
B20 D.G 3Ha F1l/Aeo D.GL 
B21 DEG. 2H Lac Till SZ.DG 
B22 D.G. 2H Fl D.GL 
B23 G. 3Ha F1l/Aeo 0O.GL 
B24 D.G 3Ha F1/Aeo D.GL 
B25 D.G. 3Ha F1l/Aeo D.GL 
B26 Bliss 3Ha Fi GO. BL 
Bay D.G 2H Lac DG.SO 
B28 D.G 2H GacatTia) SZ.DG 
B29 G. 3Ha F1l/Aeo O.GL 
B30 G. 3Ha Lac SZ.GL 
B31 G. 3Ha F1l/Aeo O.GL 
B32 G. 3Ha Fl O.GL 
B33 D.G. 3Ha F1l/Aeo D.GL 
B34 D.G. 3Ha Fl GL.DG 
B35 G. 3Ha F1/Aeo O.GL 
B36 G. 3Ha F1l/Aeo O.GL 
B37 G. 3Ha F1l/Aeo O.GL 
B38 D.G 3Ha F1/Aeo D.GL 
B39 BLG. 3Ha F1l/Aeo D.GL 
B40 G. 3Ha Lac wee 


—~ 
» @—- - 
osAVit 


“ 
Aap 


Se 
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Table B-1. Barley Sites (cont.) 


Soil Agro-climatic Parent Soil 
Site zone Area Material Classification 
E01 BL. ] Lac E.BL 
E02 DG. 2H Fl O.DG 
E03 BL. 1 Lac E.BL 
E04 BL. 2H Fl GLE.BL 
E06 DG: 3H Sg O.DG 
E07 BL. 1 BETTI Ly O.BL 
E08 BL. 1 PL/Til? E.BL 
E09 BL. 2H Tai} O.BL 
E10 G. 2H Fl O.GL 
E11 Bu. i Fl E.BL 
E13 TG 1 rr O.DG 
E14 BL. 1 PUsSTi 1 E.BL 
E15 BL. ] Resid BL.SS 
E17 Bise 1 Til BLA.SZ 
E20 BL. 2H a el O.BL 
E21 BL. 1 Til O.BL 
E22 BL. 1 eye) O.BL 
E23 LG. 1 Tibi D.GL 
E24 D.G. 2H ite O.DG 
E25 BL. 1 Lac E.BL 
E26 G. 1 Tea O.GL 
B2/ Bie. 1 Tey E.BL 
E29 eG. 1 TPL D.GL 
E30 G. 2H Wea alk Gis 
E32 eG. | Pyle HU.LG 
E33 DG. 1 Ti) } O.DG 
E34 DG. 2H Fl GL.DG 








) ee9i¢-ysiget ii | olde ‘ it 
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Table B-1. Barley Sites (cont.) 


Soil Agro-climatic Parent Soil 
Site zone Area Material Classification 
L0O1 BL. 2H Fl O.BL 
L003 D.G. 2H ER/Tili O.DG 
L04 BL. 2H Lac/Till,Lac O.BL 
LO5 BL. 2H PRsALULT O.BL 
LO6 Du. 2H Tay) O.DG 
LO8 D.G. 2H TH O.DG 
L10 BL. 2H PU/TrLI O.BL 
L11 Gu 2H Ti) 1 O.GL 
G2 BL. 2H F1,F1/Till O.BL 
L13 G. 2H er O.GL 
L14 BL. 1 FU/Spr1l/Till O.BL 
L15 BL 2H ly fae a O2BG 


2ngs 


_— 
iis 
unizo-oO7TeR oe 

ort > oe 4 






s2i2 centage 1 te 


T anes 9 ea 
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Table B-1. Barley Sites (cont.) 


Soil Agro-climatic Parent Soil 
Site zone Area Material Classification 
WO1 Ss 2A Lac O.DB 
WO02 Desa. 2A Lac O.DB 
w03 TBE: 1 je EP a a O.TBL 
W04 L28. 2A Mili Till /Resia O.DB 
WO05 BoB. 1 Lac O.DB 
WwO06 Dei. 2A Fl R.DB 
WO7 TBL. 1 F1/Till O.TBL 
wos BL. 1 Ply Tit O.TBL 
wo9g BL. 1 ma bl O.BL 
W10 TBL. 1 F1/SL E.TBL 
W12 Det. 2A Lac/S RDB 
W13 D8. 2A Fl O.DB 
wWi4 THis 1 Till/Resid,Till Oo roh 
W15 DE. 2A Lac R.DB. 
W16 Don. 2A pL Sea be O. DB. 
W17 BL. 1 rey O.BL 
wi8 BL. 3H rad O.BL 
W19 BL. 2H Lac O.BL 
W20 Ble. 2H aval ia’ O.TBL 
W22 TBL. 1 Tele Ply Til O.TBL 
W23 D.B. 2A isi O.DB 
W24 Huahine 1 Lac Ca ol 
W25 4 1 Taeiel O.TBL 
W26 TBL. 1 ed EF a O.TBL 
W27 Dis. ZA Lac O.DB 
W28 D2B: 2A Lac Sa.U05 
w29 in aving 1 juraal On Bk 
W31 D.B. 1 Fl O.DB 
W34 as 2A Fl R.DB 
W36 TRL. 1 Bi O..cb BG 
W37 BL. 1 Payot O.BL 
W38 De 5. 1 Fl O.DB 
W4 4 TBL. 1 Lac R.TBL 
W42 TBL. 1 Lac pLac/Till O~. TBL 
W43 Dias 1 ULL ON O.BL 
w44 BL. 1 Fl R.HG 
W46 BL. 2H Lac/Resid E.BL 
2H Lac, Lacy Tia D.GL 
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Table B-1. Barley Sites (cont.) 


Soil Agro-climatic Parent Soil 
Site zone Area Material Classification 
TO1 Do. 2A Fl O.DB 
TO2 TEL. ] PP erweril i OTB 
T03 TBL. 1 Fl,F1/Till,Till O.TBL 
RO Dao. 2A Fl O.DB 
TO8 duishing 1 Brest ceri? er aL 
TOS TBL. 1 FE/Ti1l],F1 O.TBL 
T10 Pe5 < 2A Fl OTDB 


Biz TBL. 1 LF MS it Re OM ape a ©. TBL 


mf 52 


“7 ha iq 
< 7 
o 
Wires 
7) 


a 






aa | 
* 3) oe 
Joyo) esri2cyeised) Shee eee 


Table B-2. 


B78 


Soil 
zone 


e 


tte e 
QAA+ AMM QQ 


QQ 


° e °. 


QQ 


NQYOIYOAQNAATDOIADAMNDOIOGOWOOMQDVIOIAHNDGVIANDOGOITDGGGHDINGVGOHNADA 


Rapeseed Sites 


Agro-climatic 
Area 


Parent 


Material 


Fl 
F1l/Aeo 
Lac 

Baca tik 
Titre bac 


Biyori. 
Lac 
F1l/Aeo 
Lac Till 
F1l/Aeo 
Lac Till 
Racer. 
acm il 1 
bac Tall 
Bac Till 
Bacarz.1 


sei 


Soil 


Classification 


O.GL 
O.GL 
$Z.GL 
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Table B-2. Rapeseed Sites (cont.) 


Soil Agro-climatic Parent Soil 
Site zone Area Material Classification 
E37 D.G 1 Lac O.DG 
E38 BL 1 RU/sTiyi E.BL 
E39 BL 2H eri? O.BL 
E40 BL 1 Sry Tell O.BL 
E41 BL 1 F1/Till E.BL 
E42 D.G 3H Tile O.DG 
E44 BL. 1 Lac E.BL 
E45 BL 1 Lac E.BL 
E46 D.G 1 Try O.DG 
E47 BL 1 PPATL LY E.BL 
E48 BL. 1 Resid BL.SS 
E49 BL 1 daw its BLA.SZ 
E53 BL. 2H 58 BE O.BL 
E54 BL. 1 gh cae O.BL 
E55 BL. 1 iL T O.BL 
E57 BL Ce 1 Tia D.GL 
E58 BL. 1 Lac E.BL 
E60 eG. 2H die at O.DG 
E61 DG. 1 Aube O.DG 
E62 BL. 1 Ti) E.BL 
E63 H.G, 2H Fl GL.DG 
E64 G. 2H es ed O.GL 
E65 D.G. 1 eal D.GL 
E66 D.G. 1 Weta ie HU.LG 
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Table B-2. Rapeseed Sites (cont.) 


Soil Agro-climatic . Parent Soil 
Site zone Area Material Classification 
L48 BL. 2H Lac/Till,Lac O.BL 
L49 BL. 2H Fl O.BL 
L50 BL. 2H Ri Tak) O.BL 
G51 BL. 2H PE/Ti 1a O.BL 
L52 Bi. 2H i Ae BI O.BL 
L53 Hie 2H F1,F1/Till O.BL 
L54 BG. 1 BPLySer l/s Till O.BL 





Table B-2. 
Soil 
Site zone 
w49 DzBs 
wW50 THES 
W52 TBL 
W53 D.B 
w54 BL. 


Agro-climatic 
Area 


ah awit at ot fA) 


A 


Rapeseed Sites (cont.) 


Parent 
Material 


Fl 
ee 
Fl 
Fl 
Tele 


Sorl 
Classification 


O.DB 
O.TBL 
O.TBL 
O.DB 
O.BL 
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Table B-2. Rapeseed Sites (cont.) 


Soil Agro-climatic Parent Soil 
Site Zone Area Material Classification 
T19 orks 2A FLi/Ti Vi O.DB 
T20 TB. 1 PET il O.TBL 
T2'1 TBL. 2H Lac,Lac/Till R.TBL 
Pa2 BL. 2A FIAT o.5 
TT? BL. 2A Raat bl OE 
T18 TBL. 2H Lac,Lac/Till R.TBL 
Tvs TBL. 1 ne ee Ge O.TBL 
T14 Blis 2H Lac RoBE 
TtT5 TBL. 2H Lac, Lac/Till R.TBL 
T16 TBL. 1 EISFU/TITI Till = §«©69O. TRL 
T26 Des. 2A Fl O.DB 
m2 Bis 2A BEATE LI O.5 
T28 TBL. 2H Bac sbac/Tili R.TBL 
T29 TBL. 1 Herel /Pi Tl O.TBL 
T30 THLie 1 BM/TLLY OSTBL 
aot (Mae 1 EPP ber Peril 1 OFTBE 


32 D.B. 2A Fl O.DB 





Table B-3. 


Wheat Sites 





Soil Agro-climatic Parent Soil 
Site Zone Area Material Classification 
J01 DLs 1 Pi Peels Til) TEL.SS 
J02 Doo. 2A Fl O.DB 
J03 be es 2A igh op Pat DB.SS 
J04 TBL. 1 Tick} O.TBL 
J06 TBL 1 Lac TBL.SS 
J07 TBL 1 Fiytrida, Fl $Z.TSBL 
J08 D.B 2A F1/Till O.DB 
Jog D.B. 2A sp ae a DB.SZ 
J10 TBL 1 FLZTill O.TBL 
J11 yeu. 2A Lac ,Lac/Till O.DB 
J12 TRL. 1 F1/S 2 lod AVE 
ais D8 2A Fl OLDE 
J14 DSB. 2A rit DB.SZ 
J15 Do8. 2A he a! S$Z.DB 
J16 TBL 1 Lac bac/Till, Till TBLCSS 
47 TBL 1 Till,Lac/Till SZ.TBL 
J18 BL 1 gieial E.BL 
J19 BL 1 Sorted Till GLE.BL 
J20 BL 1 an! O.BL 
J22 TBL 1 Ti) SZ. bu 
J23 PBL « 1 Pre bacevTil ts Thies 
J24 o.5 2A upwind DB.SZ 
J25 BL 1 Pog Tit) /Resia BL.SO 
J26 BL 1 Bieri T11 t BL.SO 
J27 BL 1 et Bi..o2 
J28 BL 1 Pill, F1/Til1 BL.SZ 
J29 BL. 1 a Gehl Oe ee E.BL 
J30 BL 1 Vi ls Wh EA Re) a el a E.BL 
J31 BL 1 iy E.BL 
a2 BL 1 Sorted Till E.BL 
J33 TBL 1 TiLisLacs hee TBL.SS 
J34 TBL 1 Past 2 Ble 
J35 Dea 2A died bil DB.SZ 
J36 TBL 1 see a TSLos.e 
J37 BL ] dia Ge Ohler ol Wide ya sa] BL.SZ 
J38 BL 1 F1/Till BL.SZ 
J39 BL. 1 ipa al BL.SS 
J40 BL. 1 Ta bd (ele ree BL.SS 





APPENDIX C 


Growing Season Precipitation, and Soil Chemical 


Analyses of Experimental Sites 


Analysis Units 

Bec. mmhos/cm? 
NO;-N kg/ha 

K kg/ha 

& indicates estimated value 


NA indicates that data were Not Available 
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Table C-1. Barley Sites 


Site Pptn. pH pH BoC, EeGarecalO. RO0M. NOs K 
(em) 0-15 15-30 0-15 15-30 0-15 C315 ee a ENO ail 
BO 1 Sed = 4.6 3 eae Cie 0.0 so Ze 405 
BO2 258 S46 al 0.4 Os 0.0 hes 39 a205 
BO3 te..5 G29 a6 fee | Os 0.0 meas 8 249 
BO4 20 ee: Le a he 3 Oe Ge0 3.6 10 365 
BO5 aes eS 4.6 es Or2 Gel Pall 2 306 
BO6 eee. Ges es Oe3 Pez 0.0 Meet 10 TRS 
BO7 300 566 4.5 Urs Cas O20 poi ie. 407 
B08 AS Be 7.4 *.9 0.4 Cus Os a4 43 149 
BO9 oes Gi 4.6 at Gel 0.0 Zed 20 330 
B10 teu Gand 4.9 ed Cus 0.0 aah 29 va ie | 
B11 f2n% ez 620 Os 0.6 0.0 ec 32 358 
Bi2 Lon o 6.0 Reo Ose a 0.0 2<9 19 448 
B13 Piad 6.4 Sac Oe O..1 0.0 4,4 27 370 
B14 16.0 7.4 7.4 0.3 O23 One 3.4 oa 459 
B15 Gres Gus Gee ee i eae 0 5.4 20 336 
B16 2G 5.4 tse 023 Se 0.0 Set sie 403 
B17 Teac Graz 29 O53 Oa 0.0 4.0 é 543 
B18 25.4 Beal 8.0 eZ 0.4 Tey a eye 8 146 
B19 Beeg Gaz ae 0.4 ed O00 Sn} 50 739 
B20 e.5 Gad 6.6 Vas as 0°. 4.8 18 15H 
B21 7.4 ac G3 Ons Oak 0.0 hes 49 448 
B22 lest mat Greet O..2 O.2 One oP: 9 403 
B23 24.1 Lat 7.8 0.4 Ges ez Set q; 370 
B24 Sad 7.4 Meo Gres 0.4 0.4 4.3 SE 336 
B25 38.1 Bad fi he vee OZ Greet 3.4 19 230 
B26 thas 6.7 6.8 O78 vas 0.0 4.0 90 1378 
B27 24.1 6.8 6.4 Ore GC xz 0.0 Zed 1 O07 
B28 18.3 6.0 er Ois2 Oa 0.0 wal 3 543 
B29 Sts Scat decer Ups! Oas ee) 4,3 9 ac 
B30 3 age yy: 6.9 0.6 G25 OV Grid 56 851 
B31 24.1 ties Sime ae Vege On) 1.4 2 Bal ot 
B32 366.5 tao Vee 0.6 ted OZ 14.7 18 409 
B33 Pap BS LAr dies 0.4 Oe 0.0 S00 ee 347 
B34 20.8* 6.4 G.2 0.4 Oras 0.0 4,4 20 969 
B35 CRS ae Sod ga 6.9 0.4 ee 0.0 ted ie 398 
B36 aap ls 6.6 a Oe2 ese 0.0 4.5 19 252 
B37 To Oe 6676 Gar 0.2 Oe 0.0 120 10 Sts 
B38 2e6 649 J oi O03 O62 0.0 acer 12 347 
B39 20.0% 7.5 be 0.5 aap! 0.6 SiS) 37 319 
B40 Dies hel G3 0.6 0.4 0.0 7.4 27 840 





Table €-1. 


Site 


E01 
E02 
E03 
E04 
E06 
E07 
E08 
E09 
E10 
Ei1 
E13 
E14 
E15 
E17 
E20 
E21 
E22 
E23 
E24 
E25 
E26 
E27 
E29 
E30 
E32 
1 £ | 
E34 


PpEn. 


(cm) 


2246 
ooo 
2euG 
41.9 
269 
Zuge 
A Waid 
£28 


Barley Sites (cont.) 
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Table C-1. Barley Sites (cont.) 


Site Pptn. pH PH B.C. E.C. %CaCO; %O.M. NO;-N K 
Kom) ~U-15 15-30 0-15 15-30) 0-15 =o 1S = 4 
LO1 3346 6.9 pe 0.3 Wes: 0.0 4.4 2 222 
L03 24.9 Gan 6.8 Upc Oz 0.0 3.8 7 377 
L04 eae 6.8 ew te is 0.0 11.4 18 431 
LO5 24.6 6.6 Gat On 2 Owe 0,0 wel se 1142 
LO6 24.9 6.6 6.8 aes. ie 0.0 Sao 29 236 
L08 23.4 Gao 6.5 Ure ge eu eeu Ze 310 
L10 $£5.0* 6.5 6.8 On2 Coe 0.0 re 8 224 
L11 Pat eas 6.4 6.6 Deo les wet 24 24 S20 
L12 Zpeue 6.7 6.8 Use aes 0.0 Sa 8 683 
L13 ho 0 6.8 Ot G3 O53 Gu 4.1 28 403 
L14 14.7 G2 6.4 Ue Oo2 Ca ee 30 of A) 
L15 to,0e 6.3 6.6 Ov C26 0.0 cy Zea 302 






2yslyae 92143 aie uy 
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Table C-1. Barley Sites 


Site Pptn. pH pH Bets E.C. %CaCO,; %O.M. NO;-N K 
ome 0-45 “15-309 0-46) J5-30) G15 0-415 0-15 0-15 
W001 i3.0 6.4 ell 0.4 0.8 0.0 Ls PR: 48 1434 
WO2 12 6.6 6.9 0.4 0.4 CeO oF, a5 1109 
wO03 8.9 6.4 Te3 0.3 0.4 0.0 shag! 19 831 
wo4 LP. 2 Soe es 0.3 0.4 0.0 Soe 21 506 
WO5 16.3 YR Te | Oe 6.8 =) 36 1413 
WO06 14.0 1 ae 7.8 0.4 0.4 16.3 a.6 65 Sty 
WO7 2.9 6.6 7.4 025 0.4 0.0 S08 78 1014 
wos 2. e.. -7-..0 ares: 0.8 Sys: Oca Gn: 78 1831 
wo9 qec0s. 7% 2 2 et 0.6 O52 Or had 29 986 
W10 28.4* 6.6 7.6 0.3 Os 0.0 626 43 504 
W12 hes ‘SIO Gee 0:6 a. 6 a0 wes VAS 1378 
W13 tOe3 vie G 8.0 0.5 0.4 0.0 2.9 Pad 1159 
wi4 eee We Sel 06 0.4 0.4 9.4 94 246 
W15 18.0% 7.6 8.0 55 0.6 0.0 4.9 54 1394 
W16 22s 6.4 dxz 0.4 0.4 0.0 Tivee, 34 728 
W17 False, Ek TAT 7.4 0.4 O53 0.0 fips 18 325 
wi18 hee Use 9.7. 4 el 0.4 OZ 0.0 Dee 20 801 
w19 28.4% 8.0 8.4 0.4 0.4 0.0 8.0 11 364 
W20 2308% 7A 13 0.3 0.4 0.0 4.8 9 560 
W22 20-51 Ga we 0.4 0.4 0.0 5.0 43 963 
W23 be el YA he Ue os 0.4 Ort 8.4 46 1042 
W24 | 0 ae, ia pS 5 OFS Bel head 29 806 
W25 2a59% 7.1 leet 0.4 Co 0.1 3 ot 35 801 
W26 20.3 6.4 6.8 Pye wes 0 ate “3 638 
W27 Siiece 7.5 oh 0.6 0.4 0 Sls 20 1635 
W28 23.46 8.0 8.0 0.6 Oe 0.0 4.1 10 1406 
W29 12.4 eg ee Urs Ov2 Oe 6 6.4 7 896 
W31 10.9 a0 hed Ge2 Oe2 Ors.2 io 3 Te2 
W34 13.0 7.4 Wie 1s 0.4 Os 1 8.4 6 ns) 
W36 oa fede, Whos 0.3 0.4 Ord is fel 10 905 
W37 16.0 6.6 Ged 0.4 or 0.0 8.9 18 946 
W38 24.6 en 720 Ce 0.4 OO 2,0 15 560 
w41 14.2 6.6 ede os Os 0.0 8.0 67 1170 
W42 14.5 Gea 6.9 Oz D hege) 0.0 cy 45 605 
W43 9.9 ight Le 4 0.6 0.4 Od OZ tz 907 
w44 18.23 o,0 3 eege 5 O55 0.4 9.0 134 246 
W46 $oe2 Bye 8.4 0.4 0.4 Oe Sti 45 291 
W47 16.0 6.9 6.9 Oe2 ee 0.0 Se 34 482 
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Table C-1. Barley Sites 


Site §Pptn.- pH PH B.C. E.€. %CaCO; %O.M. NO;-N K 
con. 80-15 15-30). 0-12 "i523 G=15 c's Dog Cg ee Mae k 
TO1 eat BED 8.0 O25 0.4 eo aac sn oo 
TO2 teed pel 7.4 Oe Cut Gel co a 69 1086 
TOS 22e3 Fak Fea 0.7 0.8 oa 4.0 36 746 
TO7 Get wel 2 a 0.4 0.4 Pee 360 36 1385 
TO08 363 6e3 BES 0.6 0.6 0.0 52d 45 1000 
TO9 6.9 6.6 lage 0.5 G68 0.0 a. 60 1280 
T10 G39 See 8.4 Oe3 0.4 eo ao 25 716 
T12 Vet ee Teg 0.4 Oe5 0.6 4.2 56 Tae 
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Table C-2. Rapeseed Sites 


Site Pptn. pH pH eG. Becvewmcaco, *0,M/-NO,-N K 
fem) “0-15 -15-30 0-15 15-30 0-15 0-15 b= 15a (C= 15 
B41 °F G.. 4.9 oo 0.3 0.0 eeu 29 277 
B42 fc. 5 ooo 5.6 0.3 Ce 0.0 e235 8 249 
B43 19.0 5.4 Se Owe Oo3 0.0 328 10 365 
B44 ao Se 4.6 ao3 One 0.0 fioee 4 306 
B45 ce pe G.2 5. 5 i.3 Oca 0.0 Best 10 703 
B46 30.0 Se ae se Os 0.0 EE 15 407 
B47 SE sya 4.6 v.34 1 pee 0.0 Sree. 27 405 
B48 Pha ee 7.4 1.9 0.4 O.3 0.3 5.4 43 149 
B49 25.7 SaaS aed 0.4 3 0.0 Cees: 34 3e2 
B50 ae: Ge sap 0.4 O22 0 4.0 34 672 
B51 16.5 625 6.6 O23 5 0.0 4.8 19 174 
B52 te. Bae 6.0 0.3 0.6 0.0 o.o 32 358 
BS3 16.5 Gna 6.5 Oe 3 os oma 5.4 22 347 
B54 1 a 6.4 36 q.2 0.4 0.0 4.4 2. 370 
B55 25.4 he 8.0 G2 0.4 Re See 8 146 
B56 19.6 5.4 6.8 ies O.7 0.0 Baal 3% 403 
B57 5 6.0 5.3 0.2 GO. 0.0 263 19 448 
B58 26.3 6.2 Bad OS Dad 0.0 4.0 4 543 
B59 7.4 a, 6 Cres res 0.4 co as 49 448 
B60 7 fe ae 6.2 Ong Oz aye) a5 9 403 
B61 24.1 yee: aoc 0.4 O53 2 a..2 fi 370 
B62 an. 2 7.4 Ce Ons 0.4 On 4.3 31 336 
B63 38.1 a3 aes o.2 0.2 O.4 3.4 19 230 
B64 whooe Git 6.8 0.8 Tao 0.0 4.0 90 1378 
B65 24.1 6.8 6.4 Die Q.2 0.0 Piece 1 907 
B66 oes Gru eat Oe ae, OE oe 3 543 
B67 SGrat ee Rieck Ons te ann 4.3 9 325 
B68 aes Lae 6.9 06 0.5 Oe Gue 56 851 
B69 24.1 tile Le: O53 Orie Ouse 1.4 2 515 
B70 3o°.5 tie G ie 0.6 rae! Osi 14.7 18 409 
B71 22.6 1 ey oe 0.4 3 0.0 SU ae 347 
B72 20.8% ~6.4 6.2 0.4 seas Ou 4,4 20 969 
B73 (oo eee Fs 6.9 0.4 ee OO ie 13 398 
B74 22.6 aaG or Oz Oi 0.0 ee. 19 the IF 
B75 19.0* 6.6 Gu 0.2 Vet 0.0 rye! 10 ey tS 
B76 22.26 SS deg Os Oe 0.0 aes 12 347 
B77 2008*) 7.5 dice US ae. 0.0 eel 2 BY 319 
B78 PHY et Ge 0.6 0.4 0.0 7.4 23 840 
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Table C-2. Rapeseed Sites (cont.) 


Site Pptn. pH pH EC. E.C. *£CaCO;, %0.M. NO,-N K 
(cm) 0-15 15-30 0-15 15-30 0-15 0-15 0-15 0-15 
E37 6.2 6.7 0.3 f.2 0.0 y ae. 3 343 
E38 257 5.9 6.0 0.2 0.2 0.0 8.6 16 223 
E39 19.8 6% ue 0.3 0.3 0.0 6.4 13 272 
E40 26.2 6.8 1: 0.3 Die: 0.0 10.7 20 159 
E41 23.9 6.4 6.5 0.6 tees 0.0 8.6 83 252 
E42 2725 6.1 6.3 0.2 0-2 0.0 9.4 1 392 
E44 22.6 5.8 5.9 G22 O73 0.0 2 52 577 
E45 233.6 6.0 6.5 3 0.4 G20 ie? 31 613 
E46 29.5 6.5 6.6 Ge2 G2 0.0 6.5 11 246 
E47 23.6 6.0 6.0 0.5 G3 aged 7.9 90 498 
E48 2129 6.5 6.7 0.5 0.6 0.0 6.0 2 370 
E49 23.9 6.3 fo 0.9 4.5 0.0 9.2 25 465 
E53 35.8 LEG 8.4 0.5 0.4 0.0 8.9 20 510 
E54 2724 738 Sa 0.6 0.5 0.0 8.6 Ws 280 
E55 267 Go 7x O22 Ong Os0 8.8 8 543 
E57 23.9 6.8 =) 0.3 0.3 0.0 4.8 44 426 
E58 241 6.0 6.5 0.4 o:.3 0.0 sibieg:! 31 560 
E60 2a-2 6.5 6.8 O02 0.2 0.0 4.0 7 342 
E61 29.0 6.6 Te? abe. 0.5 Pac Ge 12 291 
E62 24:8 6.9 G.3 0x2 Oe2 0.0 6.2 7 622 
E63 26.7 5 ae 8.0 0.5 O74 23 6.4 7 325 
E64 28.4 6.7 6.5 ee be? 620 S06 4 342 
E65 22.4 8.0 a2 0.4 Ons GeU che 6 437 
E66 24.9 8.0 eal 0.5 0.5 One Te? 4 291 
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Table C-2. Rapeseed Sites (cont.) 


Site Pptn. pH PH Be. aCe seawos %O.M.. NO,-N K 
tem)* 0-15 15-30 0-15 15-30 0-75 0-15 (15) G-15 
L48 eae! 6.8 Ger ae os O.0 10.4 18 431 
L49 S256 Chad 6.8 1m Vas sd 4,4 2 249 
L50 24.6 6.6 Gu eo 0.3 0.0 #0 18 1018 
L51 5.0% 6.3 6.6 G3 Cus 0.0 Sievt 22 S02 
L52 3,05 6.5 6.8 G22 ae 0.0 58 8 224 
L53 ee lts, 6.7 Gad ers a3 0.0 ade’ 8 683 
L54 Perede 6.2 6.4 oie Gaz 0.0 Sia 30 370 
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Table C-2. Rapeseed Sites (cont.) 


Site Pptn. pH CHU GE.C,  &.Ge *CacO, %0.M: NO,-N K 
(cm) 0-15 15-30 0-15 15-30 0-15 0-15 0-15 0-15 
wad 22.4% 6.8 456 (2°. £54 0.0 3.9 43 683 
Wee 12.8 7.4 756 0.2 013 0.0 6.4 9 879 
W52 Bee 35” ECE 4. 0.4 p20 5.2 11 756 
meee 26.G 2 .t . GR e2 7 ae SO 0.0 7.0 11 638 
wecmreed 27.9 <9 i ae 0.0 7.8 21 398 
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Table C-2. Rapeseed Sites (cont.) 


Site Pptn. pH pH a eee E.C. %CaCO, %O.M. NO,-N K 
femie O91 s —'oei0rea- >) 15-30) 0-95 ° 0-15 0-15 0-15 
T19 10.9 G4 6.6 a ar OFZ re Sat 36 a70 
T20 10.3 how Pes rez Whe 0 4.9 63 1464 
Te! Pld God Mee Ore Ov? 0.0 4.5 20 1262 
T22 LAs G22 Gas 0.7 ORS eeu Pe 40 987 
UZ Ux Py. Veh Font 0.4 0.4 eo ae) 24 889 
T18 Se: Pe Te Jeg tis o.0 ers: 11 1243 
Tis ts.8 DA 7.4 1p Ota Ocal aae 69 1086 
T14 5.8 Gel Sa J ipet 0.6 oe aed os. 1281 
gig ike 8.6 e0 hee 0.4 0.6 fede eS ae 1437 
T16 2a<9 7.4 Tea Oat 0.8 ied 4.0 36 746 
T26 he 8.0 Se 0.4 0.4 HG sia Q 36 1385 
Pai ea, Sel Sas 0.4 0.4 2.6 AeU 18 a fish: 
T28 a rer: Bee Fels: ete 0.6 0.0 4.1 ne 1576 
T29 aed Gut Lae 0.6 0.6 0.0 Site 45 1000 
T30 6.9 6.6 735 Tits: 0.8 0.0 jo FAY. 60 1280 
oS | 9.7 hee 7.9 0.4 O55 0.6 4.2 56 fhe 
Tse 9.9 Gee 8.4 te 0.4 348 Leo 25 716 
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Table -C-3. Wheat Sites 


Site Pptn. pH pH EVCe ECV G2Cac0s0° 20.M.aNOj=N K 
(cm) 0-15 15-30 0-15 15-30 0-15 0-15 0-15 0-15 
J01 14.5 5.7 Gt 0.4 NA 0.0 7.5 29 426 
J02 1425¢. 5.5 ST 0.3 NA 0.0 i ae 9 482 
J03 25.4 5.0 5.5 6c NA 0.0 5.6 40 358 
J04 44,7 5.1 5.2 0.5 NA 0.0 10.1 46 504 
J06 S6.c% 86.3 NA 0.4 NA 0.0 9.7 56 1098 
07 214% 559 NA 0.4 NA 0.0 Tes 48 737 
J08 2212-95 .6 NA Orne NA 0.0 S54 16 Tae 
J09 2023" 6528 NA agen’ NA oro 5.2 29 502 
310 36.8% 6.2 NA 0.8 NA 0.0 G4 21 1611 
53 by 18.9% 25.2 NA 0.8 NA 0.0 S47 18 533 
J12 18.392 2506 NA 0.4 NA 0.0 6.5 48 538 
13 TRC RRsES -7 NA aire NA 0.0 22 29 629 
314 TO°8e> 65.6 NA 0.5 NA 0.0 5.8 4] 414 
J15 18.8% 5.4 NA G43 NA 0.0 29 16 563 
J16 18.3* 5.5 NA 0.4 NA 0.0 5.9 36 482 
wiz 20-8%°5:6.0 NA O77 NA aig 5.0 41 346 
J18 2279428529 NA (ies NA 0.0 6.8 48 409 
J19 27294. 65.8 NA Ovs NA 0.0 722 46 265 
J20 22.9% 6.0 NA O22 NA g.0 Boel 8 316 
J22 34.3%: 5.4 NA a Bets NA 0.0 5.6 18 815 
u23 32.5 5.3 NA O25 NA 0.0 6.0 26 603 
J24 251s). 5.4 NA 0.5 NA 0.0 4.5 49 513 
J25 ‘7.3 5.5 NA 0.5 NA 0.0 sia: he 233 
J26 a2 oes 5.8 NA 0.4 NA 0.0 6.1 15 321 
J27 22.9 57 NA 0.5 NA 0.0 6.9 22 181 
J28 298-554 NA 0.7 NA 0.0 6.0 TZ 329 
J29 FOULs) 6.3 NA Fee NA 0.0 7.8 10 329 
J30 {257 5 7 NA 0.4 NA 0.0 735 67 327 
734 14.7% 6.0 NA 0.4 NA 0.0 a4 41 491 
J32 35.8 sen NA 0.5 NA 0.0 8.0 74 448 
333 40.6 5.4 NA 0.6 NA 0.0 Teas 91 702 
J34 39.1 5.8 NA G23 NA 0.0 6.3 13 764 
J35 36.8* 5.3 NA 0.5 NA 0.0 54 56 804 
J36 34.5 sere! NA 0.4 NA 0.0 Be 20 673 
137 aod) 5.7 NA 0.6 NA 0.0 6.5 77 427 
338 33.0* 5.6 NA 0.5 NA 0.0 5.4 12 2o4 
J39 1457 5.8 NA G7 NA 0.0 6.2 83 344 
J40 40.4 5.6 NA 0.4 NA 0.0 5.8 57 287 
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APPENDIX D 


Soil Physical Analysis and Textural Classification 


for the 0-15 cm Depth of the Experimental Sites 


List of Textural Class Abbreviations 


High Clay 

Clay 

Strity Clay 
Sandy Clay 
Silty Clay Loam 
Clay Loam 

Sandy Clay Loam 
Loam 

Sandy Loam 
Loamy Sand 

Sand 

Silt 
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Table D-1. Barley Sites 


Textural 
Site % Sand ye Site % Clay Class 
BO1 202 a9.6 34.2 CL 
BO2 24.2 Saou east L 
BO3 7S Be 56.0 16.3 SiL 
B04 i3.0 51.6 Bone SiCL 
BO5 26.6 SG 34.7 CL, 
BO6 Wide a 42.2 50.0 SiCL 
BOQ7 £6..0 Song oO ea SiCL 
BO8 in.2 53.0 34.0 SiCL 
BOS PAST 38.4 40.8 C 
B10 29.07 42.3 28.0 CL 
B11 21.4 41.3 eS. Gib: 
B12 20.9 39.9 a2 CL 
B13 2e<1 3552 = Pe CL 
B14 20,0) 58.0 16.0 SiL 
B15 45.1 40.4 14.5 L 
B16 tS PES Bho ao. 2 C 
Bi? rs a 40.8 aO'6 i CL 
B18 16.6 54.9 26.5 Sich 
B19 io<9 Bo a 24.9 So 
B20 25.0 Bas in Sib 
B21 fo. 2 46.9 34.9 Sich 
B22 20.0 47.4 2.0 StG 
B23 She 45.8 Ta aad L 
B24 68.6 Zhe 10.2 SL 
B25 3p ae 20.4 10.2 SL 
B26 16.9 54.7 28.4 SiCL 
B27 i.8 43.4 44.8 Sic 
B28 24.2 40.9 34.9 CL 
B29 ee EM fe are 14.8 L 
B30 1a) D2a.0 op eee | SiCL 
B31 si Are S2ne 12,3 SL 
B32 AN) 35.4 Saal C 
B33 oi bg eds. Ee ae Si 
B34 Bt 24.6 7136 HC 
B35 Git 24.1 Bee SL 
B36 Sr. 362.6 12.4 L 
B37 69.7 4.0 2003 SCL 
B38 66.5 24.1 9.4 SCL 
B39 61.4 222 16.4 SL 
B40 Eas Sar 25,7 SiCL 
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Table D-1. Barley Sites (cont.) 


Textural 
Site % Sand % SELLE % Clay Class 
E01 14.6 43.4 42.0 Sic 
E02 fa pes 50.1 24.1 SiL 
E03 '.6 eur a ae Sich 
E04 43.6 Suet ayes L 
E06 Soa 54.4 Zo63 SiL 
E07 ce a 30.1 19.6 L 
E08 PAS oP ee cee Gh 
E09 a1.9 39.4 26.1 CL 
E10 25.7 Sb6 oe ae | Sit, 
E11 74.9 14.8 U3 SL 
E13 S35 39°55 27.0 L 
E14 40.9 34.1 250K) L 
E15 3351 a 27.4 L 
E17 hay 4 43.2 34176 CL 
E20 18.0 13.3 Ores Sach 
E21 Pad Poste 43.8 2667 CL 
E22 20.0 49.5 Fo be L 
E23 S60! 36.9 FAAS L 
E24 42.7 32.4 24.9 L 
E25 16.6 40.6 42.8 Sac 
E26 so a State, 14.4 SL 
E27 PASS 43.3 30.8 CL 
E29 Seno oe ie Zid L 
E30 pe PG Show PS cet SL 
E32 ba. 9 20.0 EG SL 
E33 32. 39.4 2056 on 
E34 Sta 2 34.0 Zi CL 
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Table D-1. Barley Sites (cont.) 


Textural 
Site % Sand y gap a ts % Clay Class 
L01 fa.G 2.2 Oe LS 
L03 a ees. 29.8 12¢9 SL 
LO4 9.0 Ee eat 239.9 SiCL 
LO5 20.63 18.3 11.4 SL 
LO6 G5 oe 2348 thes SL 
L08 cy ae. 330 24.9 L 
L10 74.4 14.1 1 a SL 
L11 AS<s S5.5 ree L 
Li2 40.4 29.9 29.8 CL 
L13 54.8 20 au Lo. Si 
L14 44.6 27.9 alg (oh L 
L15 74.4 14.4 Heo SL 
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Table D-1. Barley Sites (cont.) 


Textural 

Site % Sand % Si Lt % Clay Class 
WO01 wei 29.4 6255 HC 
Ww02 9.9 Sic 58.4 c 
w03 20.6 44,3 See CL 
w04 = is 42.0 24.7 L 
wo5 Ses 1520 Sitter HC 
w06 50.0 B53 Tee5 L 
WO7 31.0 S6.2 a8 CL 
w08 20.4 47.5 32.71 CE 
wo9 Seno 45.4 Ziel L 
wW10 42.1 Soe 19.7 L 
W12 10.6 hegre Tals HC 
W13 Boece enue 18.6 SL 
W1i4 = fe lag 45.1 2320 if 
W15 f2,6 Zila 60.0 HG 
W16 2001 48.8 CPi rea Cr 
W17 240 42.1 24.6 L 
wi8 2133 $37 45.0 &: 
wi9 13.9 48.3 32.8 SiCL 
W20 S607 Saye h Han de CL 
W22 24.2 44.0 313.6 CL 
W23 oe eee | PT aa | Ti. SL 
W24 29,1 40.1 aed CL 
W25 i i a Sho aay) 47.1 e 
W26 Bled © ote pals chi ts CL 
W27 aS 24.6 iiss HC 
W28 5.4 Fa oer: 73.4 HC 
W29 22.0 46.1 31.4 CL 
W31 SGnG ae, 26.5 SGi 
W34 59.7 PNR elec SL 
W36 55.0 PRS io Tes SL 
W37 730A T 46.6 2d] L 
W38 Pisa 43.5 OG CE 
W4 1 29.4 Sietode. SHS CL 
W42 41.6 Pn da 36.9 CL 
w43 7 et | 47.4 30.9 CL 
w44 20. BO 20.6 SiL 
W46 34.2 B 2 aod 
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Table D-1. Barley Sites (cont.) 


Textural 
Site % Sand % Silt % Clay Class 
TO1 ret 47.1 Sige L 
TO2 Pe Pah 45.6 3h CL 
703 28.4 42.8 28.8 CL 
TO7 S262 47.7 266.1 L 
TO8 2e05 47.0 S027 CL 
TO9 28.4 43.3 2602 CL 
T10 S164 48.2 202 L 
T12 26.9 41.1 30.0 CL 





Table D-2. 
Site % Sand 
B41 29.7 
B42 P9 
B43 PaO 
B44 26.0 
B45 12.5 
B46 EG .0 
B47 Auece 
B48 $552 
B49 24 62 
B50 1 6 
B51 29.0 
B52 20.4 
B53 45.1 
B54 hie: 
B55 £646 
B56 18.3 
B57 20.9 
B58 vie oem 
B59 TE 2 
B60 260 
B61 i We 
B62 68.6 
B63 69.7 
B64 16.9 
B65 tc 6 
B66 24.2 
B67 Ct ae 
B68 TG.) 
B69 55.5 
B70 ae 
B71 oh ee 
B72 me 
B73 Gihes 
B74 Sec 
B75 69.7 
B76 66.5 
B77 61.4 
B78 16.3 


Rapeseed Sites 


SLL 


ee 
56.0 
bea 


% Clay 


28.0 
Gis 
Ko ier 
34.7 
oS Gets, 
aoa 
34.2 
a1.6 
PE a 
24.9 
7 ae. 
oY Pe 
12.5 
Es 
20.5 
43.2 
S22 
36.1 
Je en 
2h. 
14.7 
10.2 
00 a2 
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Textural 
Class 


7 ; V5 
a - 7 
7 : 








Lé 
a] os i Grd]. 
* 1 


192 


Table D-2. Rapeseed Sites (cont.) 


Textural 
Site % Sand % Svtt % Clay Class 
E37 1 is 48.5 33,7 SiCL 
E38 2507 42.1 342 CL 
E39 31.9 39.4 FA: Pe CL 
E40 50.3 30.1 19.6 L 
E41 35.5 3h.9 20.0 L 
E42 22.3 54.4 25.4 SiL 
E44 10.6 > i Poy Mee SiCL 
E45 14.6 43.4 42.0 Sic 
E46 1 oe Wie. 39.5 27.0 L 
E47 40.9 34.1 Sr apae L 
E48 oS S400 27.4 L 
E49 te ee 43.2 ano CL 
E53 18.0 43.9 Ss SiCL 
E54 ae 43.8 25.7 CL 
E55 ZOwo 49.5 2859 L 
E57 S289 BRIS rat BPS s L 
E58 16.6 40.6 42.8 Sic 
E60 42.7 32.4 24.9 L 
E61 a24) 39.4 20.6 CL 
E62 Fee athe 435.3 40,6 CL 
E63 36.2 34.0 21.0 CL 
E64 he ge Bee ton SL 
E65 36.1 36.9 27-50 L 


E66 See 20.8 16.73 SL 
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Table D-2. Rapeseed Sites (cont.) 


Textural 
Site % Sand % Silt % Clay Class 
L48 19.0 BAe 29.9 SiCL 
L49 78.6 t2.2 Gee LS 
L50 1.3 1G.3 124 SL 
L51 74.4 14.1 iE ehe SL 
L52 74.4 14.1 Lae: SL 
L53 40.4 29.8 29.8 CL 
L54 44.6 Zhe Zhe L 
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Table D-2. 


Site % 


w49 
W50 
W52 
W53 
W54 


Rapeseed Sites (cont.) 


Sand 
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Textural 
Class 


Taal) 
: i: x 


t1G>) 25712 Sebesgen), st ena 










Lee on4az 7 —~eat 
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Table D-2. Rapeseed Sites (cont.) 


Textural 
Site % Sand oe it % Clay Class 
T19 30.4 go.5 oT of CL 
T20 34.1 39.0 2679 L 
r24 10.4 26.9 62.9 HC 
T22 S732 S729 24.6 L 
ated pot oe 41.4 246.7 L 
T18 9.5 46.3 44,4 Sic 
£43 eon 45.6 a ieo CL 
T14 Pe a Pe, 6972 HC 
e15 12.4 26.6 617.0 HC 
T16 28.4 42.8 fo 00 CL 
T26 32.2 47.7 20.41 L 
E27 Se as) 24.1 SiL: 
T28 me.8 20 6:1 63.1 HC 
729 o2a3 47.0 30.7 CL 
30 28.4 43.3 28.3 CL 
631 28.9 41.1 Sto0 CL 
re ae 48.2 2057 L 


i 
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Table D-3. Wheat Sites 


Textural 
Site % Sand £°Silt % Clay Class 
J01 38.3 a9. 22.0 L 
J02 69.5 1622 14.2 SL 
J03 30.7 40.9 20.4 L 
J04 29.6 47.6 2250 CL 
J06 20.8 40.6 Gish CL 
J07 39.6 40.6 19.8 L 
J08 58.9 24.7 1625 SL 
J09g9 42.9 34.7 22.4 L 
J10 3052 41.5 28 CL 
J11 36.4 34.9 20<7 CL 
J12 48.7 30.8 2035 L 
J13 Bo.7 apie S060 CL 
J14 44.8 ee Fle pete L 
J15 se fees Zeek 22.4 SCL 
J16 25.9 39.1 ch eal CL 
J17 40.6 38.9 20.5 L 
J18 oN ed | 34.1 34 1 CL 
J19 37.8 Set oe Let CL 
J20 2.3 e525 37.4 CL 
J22 So ees! Skeet 2655 L 
J23 Pa bees 42.3 S60 is 
J24 40.7 Spe eS ae le L 
J25 46.8 206 2025 L 
J26 40.6 34.8 24.6 L 
J27 ae2.3 3951 24.7 L 
J28 30.4 41.1 28.8 Ci 
J29 35.9 Sy oe 26.9 L 
J30 44.5 2 ENG, ioe L 
ait ZG ii 42.4 28.9 CL 
J32 26.3 40.5 ase CL 
J33 24.8 46.4 oO CL 
J34 Bee adie 257 L 
J35 36.4 39.0 24.7 L 
J36 eeu sdeeu S059 CL 
J37 44.6 34.9 20.20 L 
J38 40.7 36.8 Cea L 
J39 te | 43.2 the ret | L 
J40 39.4 osu 22.6 L 


924 cnemd Gkee ahha? Ss | 


© s Le ' 
Ps 2% pa! e 
; = te 
= aut ve 
’ } 7 . 
— 
g.¢ rq 
> 
& j — } 
mal 
‘ a § 
= t 
’ — 
=~ ad 7 
‘ ; 
| in 
ie 
j ; 
a8 3S 
: Foe i,o° 2¢ 
: aa 
6.02 a j 
Ear Tz 





197 


APPENDIX E 


Available Phosphorus (kg/ha) Analyses of Experimental 


Sites 


NA indicates data were Not Available 
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Table E-1. Barley Sites 


Site ASFTL-P M & A-P Olsen-P 

00> La)---UTSe30) (0=95)) Ane=30) LO=T5)> (15-30) 
B01 34.0 0.0 305% 4.5 2300 20.2 
BO2 2000 0.0 31,4 4.5 24.6 2042 
B03 49.3 23.5 7 Paes 23 eh 29.4 eae 
BO4 Ct S.9 gaa 0 14.6 Va Pla 
BO5 13.4 0.0 14.6 0.0 340 4.5 
B06 o2; 0.0 68.3 Bini 40.3 6.7 
BO7 37.0 0.0 3255 222 A Wa2 
BO8 Oe 0 0.0 6.7 1.1 5 Pa 9.0 
BO9 AZo 0.0 16.8 ane 16.8 23a 
B10 the 0.0 ae RE 9.6 v4 ro 3352 
Bij 295) aia ee 26.0 a. 6 ee | 14.6 
B12 vee Zink ao Zee o20 12.0 
B13 13.4 a We 3.4 sa Pay 12.8 
B14 11,2 66.1 153.4 i026 12254 65.0 
B15 SEAL Zee ry ea 12.3 rare, 
B16 31.4 3 ot GP ars 9.0 S20) Les, 
B17 24.6 | 28.0 4.5 35.8 toe 
B18 2 Pel Guat 0.0 2052 10.49 
B19 ier eee 16.8 rao 13 Hise 
B20 ee) 13.4 (Ae fate de As 10.1 
B21 49.3 aa0 54.9 14.6 2560 RS 
B22 68.3 J ee oe ee tase 50.4 13.4 
B23 115.4 vo ee laZeo Oda 2 59.4 40.3 
B24 26.9 L203 25<0 12.6 os.0 0 
B25 14.6 Vou LOS 16.8 1 Pe 9.0 
B26 30.2 Ti Ea 14:4 3854 31.4 
B27 hs Ive 44.8 ees 5O a0 22% 4 
B28 Bak Peat Pees oie 2 iZa3 eat 
B29 1s 5.6 nou ares 90 4.5 
B30 haw 14.6 16.8 10.41 16.0 9.0 
B31 108.6 2026 Oia 20075 Soh ace, oon6 
B32 Papas 5.6 2a 4.5 25.8 134 
B33 Se BS toes CS egre. 11926 AD ae Sexe 
B34 a0 0.0 14.6 fins 26.9 13.4 
B35 134.4 125.4 13129 144.5 62.7 98.2 
B36 teas 3.4 20.2 isn Ue, 4.5 
B37 218.4 2.1ge3 22406 207.2 Lone tit prs} 
B38 2.4 3259 S320 31.4 22.4 ies 
B39 Ure ara | ue 134.4 1.8 POG 2 5.0 
B40 19.0 ie PLES 7.8 26.9 9.0 
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Table E-1. Barley Sites (cont.) 


Site ASFTL-P M & A-P Olsen-P 
CO=157),52C 15230) (8-15)  Cis=30) co= 1S “i5-30) 
E01 Sage 22.0 78.4 7.8 7a62 Rae 
E02 78.4 dens. eae 3.4 40.3 ie ay 
E03 40.3 4.5 3.8 Le 31.4 Hered 
E04 34.79 15.7 23.5 Meet 24.6 S022 
E06 40.3 1am 38.1 bi ae sae 2hed 
E07 26.9 14.6 24.6 at a ad 13.3 
E08 26.9 ou 14.6 a. 6 ae 6.7 
E09 28.0 3.4 22.4 14.6 oP 16.8 
E10 26.9 12.3 2052 3.4 12.0 14.6 
E11 20.3 7.8 Leo 10% $5.67 13.4 
E13 3564 4.5 40.3 220 2840 1 giave 
E14 aoe Veal 65.0 Pp as] 40.3 9.0 
E15 40.3 ba 2 8 fe 0 rae ee 13.4 
E17 48.2 Meet 49.3 3.4 aa60 Le ee 
E20 1 Pe 0.0 14.6 0.0 Pam Sire, Wd a2 
E21 Sed 0.0 1.9 ae 1S20 14.6 
E22 20.0 2.6 22.4 eehe Pikes bh e 
E23 3.4 1,1 ee hall 7.8 4.5 
E24 oa Bee Lon 2.0 Lond Jeu 
E25 10 tac 1220 7.8 pe ae ol 
E26 40.3 flies as Bee 2652 ora 26.9 
E27 a 24.6 53.6 Z0eG 20.0 aes 
E29 S520 eZee Lae, 34.7 31.4 42.6 
E30 eae Pipes 54.9 20.9 2350 22.4 
E32 mec 0.0 a: ia 14.6 oe 
E33 ted oye 14.6 Gal Mise g.0 
E34 io Pa 2 wil 240 4.5 Rise 9.0 
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Table E-1. Barley Sites (cont.) 


Site ASFTL-P M & A-P Olsen-P 
LOrte se 15-30) (Hots) .tees0) (O-d5)> ohiSes07 
L01 Te «2 82,9 66.1 68.3 38.4 31.4 
L03 98.2 Je=d 61.6 38.1 302 22.4 
L04 28.3 ie ae, 31.4 14.6 38.9 2620 
L05 44.8 22.4 40.3 16.8 a 14.6 
L06 a8.8 Piso 34.7 Pee, (Wy ie) 14.6 
L08 26.9 28.0 3.4 24.6 pape ee 19.0 
L10 14.6 +e NA NA 14.6 NA 
L11 59,0 0 NA NA 16.8 NA 
L12 29.2 1g NA NA 29.3 NA 
L13 40.3 6.7 NA NA 35.8 NA 
L14 13.4 13.4 NA NA 14.6 NA 
L15 hex 2A NA NA 13.4 NA 
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Table E-1. Barley Sites (cont.) 


Site ASFTL-P M & A-P Olsen-P 
(0—15).--( 15-30) CG—7S)) = C1530) (O15)? ChS—3G) 
wo1 FA) Grtk 40.3 Th 44.8 2755 
wo2 2955 D2 31.4 4.5 33.6 14.6 
w03 2103 0.0 28.0 5.6 26.9 24.6 
wo4 20.2 Broek 26.9 12.3 23.5 14.6 
wos 14.6 304 13.4 7.8 33°.6 NA 
W06 4.5 0.0 5.6 0.0 13.4 13.4 
WO7 29.1 Ze 33.6 9.0 39.2 yon. 
w08 69.4 262 hie 33.6 S54 £307 
wo9g 40.3 Gx 44.8 ee: 48.2 29.1 
W10 ce 0.0 33.6 Ge 35.8 14.6 
W12 4.5 0.0 5.6 6'.7 24.6 2a 
W13 aS ae Co 23.5 Rig! 24.6 15.7 
wW14 $59 4.5 16.8 fea 29.1 35.8 
W15 to 0.0 26.0 are ey, Tae 
W16 44.8 202 26.9 4 Sled 33.6 F253 
W17 ee.” raiey 15.7 9.0 24.6 20 22 
w18 tad 0.0 hiee 5.6 20.2 13.4 
wi9 14.6 0.0 1 aes 0.0 224 19.0 
w20 101 0.0 13.4 Ong 22.4 Eu 
W22 22.4 7.2 26.9 13 S022 thc 
W23 7a: bay ada, 28.0 a4 224 9.0 
W24 s0n2 Pha 35.8 9.0 33.6 13.4 
W25 a4 4.5 9235 Aad 16.8 134 
W26 35.8 4.5 59.4 202 33.6 17.9 
W27 a | Gas 3720 14.6 41.4 26.9 
W28 19.0 a, 22.4 ent 35.8 iat 
w29 41.4 22 48.2 9.0 29.41 13.4 
W31 47.0 10el 51.5 15.7 Beso (ed 
W34 20.2 7.8 19.0 17.9 157 228 
W36 69.4 15.7 57.1 Dias 35.8 228 
W37 45.9 6.7 67.2 21,3 42.6 22% 
wW38 65.0 24.6 75.0 336 69.4 44.8 
W41 61.6 piste y: NA NA NA NA 
W42 wae 112 NA NA NA NA 
W43 G2 16.8 NA NA NA NA 
w44 14.6 5.6 NA NA NA NA 
w46 728 3.4 NA NA NA NA 
W47 Ps ew | 19.0 NA NA NA NA 
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eo7l2 geisha 1° elgeiee 
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Table E-1. Barley Sites (cont.) 


Site ASFTL-P M & A-P Olsen-P 
.0-15). °° (15-30) (0-15) (15-30) (0-15) 9115-30) 
T01 22-24 6.7 24.6 16.8 31.4 Leo 
TO2 63.8 5.6 86.2 7.8 67.2 9.0 
T03 25.8 5.6 35.8 10.4 35.8 17.9 
TO7 56.0 1422 33.6 4.5 35.8 9.0 
TO8 63.8 14.6 2447 t2s3 62.7 13.4 
TO9 sist (ein 3.4 70.6 5.6 58.2 9.0 
T10 40.3 4.5 40.3 14.6 49.3 13.4 
T12 26.9 3.4 3226 17.9 35.8 13.4 


1-8 afdet oe 


) Sed 
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Table E-2. Rapeseed Sites 


Site ASFTL-P M & A-P Olsen-P 
Pueior? (15=30) (Os 15,.2 345730) LO hog) ste OC) 
B41 Toaeu 0.0 eae 5.6 Ya Soe 
B42 S23 2300 oe ee 20 34 2 ail eee 
B43 34.7 6.7 37.0 14.6 22.4 1672 
B44 13.4 0.0 14.6 0.0 a. 0 4.5 
B45 62.7 0.0 68.3 ene 40.3 Ga 
B46 37.0 0.0 is ge Ante 2560 Bey 
B47 Sf.0 0.0 3.4 4.5 23.0 anes 
B48 a) 0.0 6.7 Lat te00 wa 
B49 26.9 0.0 31.4 4.5 24.6 20.2 
B50 13.4 3.4 16.8 a. 6 24.6 Ue eae 
B51 5.6 0.0 6.7 3.4 13.4 9.0 
B52 29 cal 22 2850 Sg 34.7 14.6 
B53 =, at ihe eh f2.3 T2203 
B54 13.4 Tie 12.6 3.4 Lbs Hy? ee, 
B55 9.0 Tt 6.7 0.0 20.2 10.2) 
B56 oh 3.2 3220 2...0 g2.0 oY | 
B57 20 2a coe ene oe eae. 
B58 24.6 3.4 28.0 #5 3060 1S 
B59 49.3 2.0 ei tee, 14.6 Zone Wing 
B60 66.3 1S 62.2 Poe 50.4 13.4 
B61 TiS. 4 Tikian Lizec 67.2 59.4 40.3 
B62 26.9 ae) 2550 14.6 i708 a.0 
B63 14.6 ow 1654 16.8 1152 g 0) 
B64 3052 ey’. S255 13.4 36a oe 
B65 Nae, Si. 44.8 tg) DOU 22.4 
B66 2.0 Lait Zac Ode? Tana o20 
B67 ea 5.6 a0 4.5 chau t) 4.5 
B68 pens, 14.6 16.8 10 cI 16.8 S20 
B69 108.6 201.6 101.49 200.5 5155 $976 
B70 BR apie, sir de Pile le} 4.5 2060 13.4 
B71 Gs Goh ata bY PX i Oodosces Bios Tone oT 0 bees 
B72 a0 O08 14.6 iS, 26% 13.4 
B73 134.4 12 ae UES PPA: 144.5 open | oe7e 
B74 2 Rea: 3.4 7h ae Pas Lies 4.5 
B75 218.4 eM BOF | 22a ae A0T a2 Diese Pideeas 
B76 31.4 ate aye ens 31.4 22.4 thee 
B77 Lehaad 3.4 134.4 fi Te oe aa 
B78 19.0 Bak othe wit 2600 ou 
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Table E-2. Rapeseed Sites (cont.) 


Site ASFTL-P M & A-P Olsen-P 

CO= 1552" 15-40) VOetS)) 2085" 30) COTE Pe ea 0y) 
E37 ge 0.0 2350 3.4 a 16.8 
E38 26.9 pe! 14.6 5.6 awe 6.7 
E39 28.0 3.4 22.4 14.6 2349 16.8 
E40 20.2 14.6 24.6 950 $052 13 
E41 33..5 9.0 390 Tae 29.6 16.8 
E42 40.3 Laat 30m 13.4 S062 Zales 
B44 40.3 £.2 43.7 7.8 31.4 1223 
E45 tae dee! 70.2 (hye Io.ae (ets 
E46 3041 4.5 40.3 sis See tate 
E47 oi Wes 1.1 in) 7.8 40.3 site 
E48 40.3 fae 3022 ae Zito 13.4 
E49 48.2 Vee) 49.3 3.4 33.6 dite 2 
E53 i. 0.0 14.6 O70 2350 Love 
E54 6.7 0.0 TERS: ieee tou 14.6 
E55 ZO. 5.6 ace S a 0 FP eis] 
E57 35.8 22.0 Ao. 34.7 Sl. = 42.6 
E58 ioe Hires phe ley & 126 To rat 
E60 Lond Ces 13.0 9.6 16.8 ei 
E61 Tees. Zee 14.6 Ga bei eee t, 
E62 oe ey) 24.6 ee i= 2220 25.8 fs Ae, 
E63 92.0 242 OU 4.5 Nee 9.0 
E64 oo. 2 Zana 54.9 26.9 Pa eRe) 22.4 
E65 3.4 trast 4,9 1.1 120 2. 
E66 7.0 0.0 L253 ms. 14.6 as, 
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Table E-2. Rapeseed Sites (cont.) 


Site ASFTL-P M & A-P Olsen-P 
(0-15): (76=30) 0-16) 415-30) LO-ial= <15-30) 
L48 38.4 To Sie4 14.6 c1s5 2an5 
L49 49.3 38.1 42.6 24.6 25 <5 Zee 
L50 45.9 14.6 34.7 33 S26 720 
L51 ore, Zins NA NA 13.4 NA 
L52 14.6 neo NA NA 14.6 NA 
L53 aoee 1S NA NA Pe NA 


L54 13.4 13.4 NA NA 14.6 NA 





a 
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Table E-2. Rapeseed Sites (cont.) 


Site ASFTL<P M & A-P Olsen-P 

VO 19) ./( Te=30) CO=tS) (15330) CO)" Aa 5e30) 
W49 26.9 Siok NA NA NA NA 
W50 40.3 250 45.9 13.4 356 13.4 
W52 69.4 fied PACee 2620 47.0 Zen 
W53 do6U 2620 80.6 a7.0 a 40.3 
W54 10.1 1054 13.4 Ghee 22.4 22.4 
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Table E-2. Rapeseed Sites (cont.) 


Site ASFTL-P M & AP Olsen-P 
F0215)-* ¢ 15-30) (O=tS)e GS-30) fO=15)" (15-30) 
Lis BE.2 2.8 28.0 10.1 ie Pe 4.5 
T20 42.6 bat 59.4 24.6 44.8 iqee 
ed 42.6 4.5 47.0 CRS) 44,8 920 
ae 23.6 Gare 3.4 1a) 22.4 oe 
mye 22.4 4.5 fleets 3.6 oeeo meu 
T18 24.6 ane 24.6 9.0 Saat 132 
abs 63.8 aa0 86.2 Fars, Gfae ele |) 
T14 x<G al 4.5 0.0 40.3 AES 
2t5 A a0 aaa! 14.6 40.3 VES 
T16 2a60 aa6 a0 TOrnt cL ete Teg 
T26 oie eat Thaz a20 o5 Chel s B30 
TA? ho. 4 Bitte Wiel rt Cae 4.5 
T28 43.7 eee 42.6 3.4 49.3 4.5 
729 Gad 14.6 Ta Pe: Satan 13.4 
T30 pein | 3.4 710,..8 oe 6 Sh e2 gD 
St 2640 6 22.0 19 3586 13.4 
32 40.3 4.5 40.3 14.6 49.3 134 
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Table E-3. Wheat Sites 


Site ASFTL-P M & A-P Olsen-P 
(0-15)57°(15=30) (0-15) e415 —30) CO-157)5 (15-30) 

J01 ie 29 34.7 NA 26.9 NA 
J02 15.7 4.5 R770 NA 2208 NA 
J03 52.6 3245 84.0 NA 49.3 NA 
J04 28.0 104 51.5 NA 35.8 NA 
J06 38.1 NA 4740 NA 44,8 NA 
J07 S720 NA 39.2 NA 22a NA 
J08 35.8 NA a4 27 NA iv. NA 
J09 le Say NA A730 NA 26.9 NA 
J10 116.5 NA 128.8 NA 62.7 NA 
ont 43,7 NA 42.6 NA pars ae NA 
a2 ard NA 33.6 NA 22% NA 
J13 525 NA 33.6 NA 22.4 NA 
J14 38.1 NA 56.0 NA 31.4 NA 
J15 51.5 NA 50.4 NA 26.9 NA 
J16 56.0 NA 56.0 NA 40.3 NA 
O17 52.6 NA 58.2 NA 203 NA 
J18 i bay: NA 16.8 NA 13.4 NA 
J19 14.6 NA 2072 NA 17.9 NA 
J20 ee NA 13.4 NA 12.4 NA 
J22 6227 NA 63.8 NA 35.8 NA 
J23 107.5 NA 12120 NA ied! NA 
J24 43.7 NA 50.4 NA chia NA 
J25 aad NA 4-7..4 NA 26.9 NA 
J26 69.4 NA 73.9 NA 40.3 NA 
i) 38.1 NA 42.6 NA chee NA 
J28 3720 NA ie NA 26.9 NA 
329 28.0 NA 42.6 NA 26.9 NA 
J30 24.6 NA Sit NA 26.9 NA 
re 22° 4 NA ENG NA 22:58 NA 
332 26.9 NA 49.3 NA 31.4 NA 
J33 104.2 NA 98.6 NA 67.2 NA 
J34 54.9 NA 2e8 NA 40.3 NA 
J35 70.6 NA 63.8 NA 44.8 NA 
J36 11220 NA 107.5 NA 53.8 NA 
337 49.3 NA 35.8 NA 35.8 NA 
338 44.8 NA 25.8 NA S44 NA 
339 66.1 NA 28.0 NA 40.3 NA 


J40 41.4 NA a6 NA 31.4 NA 
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APPENDIX F 


Mean Yields (100 kg/ha) of the Phosphate Treatments 


for Experimental Sites 


check indicates site yield for nil fertilizer 
treatment 
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Table F-1 Barley Sites 
P,0; Site 
Rate 
kg/ha BO1 B02 B03 BO4 BO5 B06 BO7 B08 BO9 B10 
eneck 23.4 23.9 oo, we ee 20 ao SLO eee sete cee eet 
0 eae 6 21.5 Se eee Molen io oe. GeO 24a) 2286 
11 
17 ieneleh 2. 1S. Se eea ete. t Stee SU. Oe Se eo chon eae ot 
22 
28 
34 Saree 60s hs ie eames  aae, 939,55," 47,9: 936.5 7.33.5 
45 
50 Pome oo) toca eee ie a toed , 35.46 52.9). 20.6735 20 
56 
67 png sf of 4 Deo fa. 8 4245.9 95325 983657 3763 
84 
90 
101 Beep de OS 2 ee eo Ale SF ee. 56.6" S72. 34.6 
134 
Reps 3 3 3 3 3 3 3 3 3 3 
P30; Site 
Rate 
kg/ha B11 B12 B13 B14 Bis B16 B17 B18 B19 B20 
Peectecuao. 15.3 “2608 (2927595562 2766. 15.9 4059. 25. 58 
0 Sao beh kee ee eS ee Oreo tS ee SS 7il er ees eee 
11 
Av S4,0..2/.4 39.2 (32.6) 044.0 .35.2 (39.8 S203) 2 5. Ome ee 
ee 
28 
34 A361 28.3: 43.5" 13099893 042, 8 eso, be Sei ee 
45 
50 46.9 33.7% 4401 “2S eto 5) 4, 7 2G oe. ae. ce 
56 
67 45:5 35.8) 490675 9055 4) 40-7 i eo eo ee 
84 
90 
101 A722) 34.09 49, Sees A Ol) TO 2 ey 3S. S/.6 257 
134 
Reps. 3 S 3 3 3 3 3 3 3 3 
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Table F-1. Barley Sites (cont.) 


P,0; Site 


kg/ha B21 B22 B23 B24 B25 B26 B27 B28 B29 B30 








Reps. 3 3 Ss 3 3 2 3 3 3 a 
P,0; Site 
Rate 





kg/ha B31 B32 B33 B34 B35 B36 B37 B38 Bao B40 


check 11.3 au (Af sie ( Sao oe tee 20.6. Te Tao: 9458) al 
0 wee 49.5 “SSO Cae, oes S79 ° 18.5 4098 3G 2 Sto 


+ 2509 27.6 ~ST2 Se.G 346.6 S7.0 20.5 3h. 0s ee en 


34 22-0 19.6 49.3 S24 SY Sa eI S99) eee 
50 2¥.3. 23.9 53.9 “Sige ao6. 53a 2a 08 ao) ee eee 


67 16.5 7.6. S2V2e2aG. 2 52.0 505° 25.37 41.) AG. elas 


101 22.5 -24.9 54.0 soar SO.h S89 aie 4052 SO.0 40.50 


Va 
es 
ellis < on 
en oe 
~ —e 






orn 10.08 .C.08 c&.f+  fivtd 22 ¢€.28. B.4f C4 


ee 2 ee 
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Table F-1. Barley Sites (cont.) 


P.O, Site 





kg/ha E01 E02 E03 E04 E06 E07 E08 E09 E10 E11 





merireo.G 46.2  Oi.c" gece du.0 98,5 19.4 25.8 11.1 18.9 
0 Bee, 3065) Soe oo eee ou 22.8) 135.2 20.8 20,4 


17 Soro, 939.5, Ss0nee Sone Goono 27.0 30.8 -36.4 35.2. 26.4 


0) 
~) 
Ww 
~ 
i 
WwW 
wo 
. 

ine) 
WwW 
ine) 
co 
WwW 
~~] 
a 
WwW 
oo 
e 

(as) 
Nh 
~) 
O07 
W 
\O 
e 

co 
i 
—, 
\O 
Ww 
~] 
~] 
N 
\oO 
0) 


101 BeeGe fone 028 “8955 5a. 29.0 47.6 43.7 36.3 26.3 
134 


Reps. 3 3 5 3 3 3 3 3 S: 3 
PlO; Site 
Rate 





kg/ha E13 E14 E15 E17 E20 E21 E22 E23 E24 E25 


eneck 34.6. 41.7 Deh me DO MOMeeeuEr oe. d 2b. 0 Shiy o TO lias 
0 okt Si. a Oe ee eon oe sae 2). dae oe Ones ae 
11 Ai OA Ae he ieee, 26250 S10 127i ee oe eee 
17 BBG: 39.35 O27 SG 32257 628.4) 30.) 374) (28k sae ee 
be S80) 444° (26225 25.50 2922 3254) (34.7 (28 eee eee 
28 
34 4325 38.4 -2907) 133508 29.7. -33.6:0939 308 (305225 Oe 
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F-1. Barley Sites (cont.) 
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Table F-1. Barley Sites (cont.) 


P.O. Site 


kg/ha WO1 WO02 W03 w04 WO05 W06 W07 WO8 wOo9 W10 


Beegeroce 36.2 20,6 20.0) S166 49.8 .35:3 30.0 14.9 19.4 
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Table F-1. Barley Sites (cont.) 


P,0; Site 
Rate 
kg/ha W23 W24 W25 W26 W27 
meee 80.8 . 27. 4.95.4. 20 ae 
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Table F-1. Barley Sites (cont.) 


P.O; Site 





kg/ha TO1 TO2 TO3 T07 TO8 TOS T10 Tt2 
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Table F-2. Rapeseed Sites (cont.) 
BO. Site 
Rate 
kg/ha B61 B62 B63 B64 B65 B66 B67 B68 B69 B70 
check 2.7 4,3 5.5 5.4 eat 4.7 U.7 Seo Peo ead 
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Table F-2. Rapeseed Sites (cont.) 
Paves Site 
Rate 
kg/ha E37 E38 E39 E40 E41 E42 E44 E45 E46 E47 
check 1.6 5.0 8.5 7s tO.) 81079 4 Nees $2.97 e208 
0 10.6 1. Sed Sem cell 14.2 met Weise 15.0 19.4 
Be 20-5," 2057 
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Table F-2. Rapeseed Sites (cont.) 
P,0, Site 
Rate 
kg/ha E63 E64 E65 E66 
check 4.5 G25 oS eo 
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F-2. Rapeseed Sites (cont.) 
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Table F-3. Wheat Sites 














PyO; Site 
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check Grae cote. le. o. «2a. -26.2> 910.9 
0 Tere Lh. Oo) he eee eee oes 25.56 6166.6) 628.766 634.0. 14,0 
11 
17 liane 
Ze 
28 awe) ae ee 9.9) U2 sei 2 
34 18.9 
45 0 paar Oe Wipes a lly IMM ten A | Se A. 
50 2054 
56 Ammer U, UO. ehesboo . oilad) eno 
67 poooe 21.6 1400+ 25.6 20.4 
84 Be Al Los ee oat) Oe e 
90 
101 
134 BO G2 2208 ohSsalo- 2405 
Reps 3 3 3 3 3 3 3 3 ¢ 4 
Pau: Site 
Rate 


kg/ha J12 J13 J14 J15 J16 he? J18 J19 J20 J22 





Bue eee. 20.4) 20s | are estat: pee. 18,6 “22 12.7) 1338 
0 Bea 262 de Lose MO eee Oe ko.) 2003) 9 Ol ee ee 


17 2h 0 26.0 24.0 boeUee2578, 25.6. 52959 . 20 ee ee 


34 2528 31.2) 25.48.0009 e20.6 26.56 35.4 22 Cee ee 
50 25.3: 31.5) 26.4 (2a eb 8 26. el ee ee 


67 25.3°. 3052 “22a ae, 0 2A BSG. oo ee care oie 






d ‘ 
< = 
f 
‘ 
~ 
' 
~ ~~ 
/ 
- ea 
> 
- P ene 
> . e ¥- a 
ro 


206 O85 o.8e Mee Kae es 2, eS | SE oe 
| 5. BS. 






- * : = 7a 
a q i aa tux 


AS MOS a2 At AEE Bo 


225 





Table F-3 Wheat Sites (cont.) 
POs Site 
Rate 
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APPENDIX G 


Second Order Polynomial Coefficients 
for Experimental Sites 
Note Coefficients calculated on the basis 
of mean treatment yields. 
R? indicates goodness of fit. 
* Tstonificanteae 1m. = 0205% 
#40. Sugnificanteatep < 0.07; 


Lack of significance is due to the 
number of treatments at the site. 


Units for b, are 100 kg/ha. 
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Table G-1. Barley Sites 





Site Be By be R? 
BO 1 26.92 0.176 =. 00055 0.95% 
B02 A ae OC wt74 = 700126 O..16 
B03 toaro. -O.008 0.00020 210 
B04 20.46 —-0.048 0.00068 0.14 
B05 32.205 Oat = -0G083 0.91% 
B06 10752 0.100 =@,0009% 0.42 
B07 Cy Pees. 0.249 =. 00192 0.54 
BO8 res i 0.951** -0.00659%**x 0.98*x 
BOS 24.94 256 =0.00TS6 Org: 
B10 25-02 0.349 =0 60276 0.86 
Bi1 32.40 Oe 3S -0.00214 0.88* 


B12 17.48 0.472** -0.00346* 0.95% 
Bt ae 0.347** -0.00205* 0.97% 


B14 zo.de —~0,006 0.00092 0.66 
Bt 33430 0.561** -0.00432x 0.96** 
B16 32.64 0.253 -0:.00192 0.80 
Bre See?) -Do1e9 0.00196 0.49 
B18 2a abe 0.415% =O0,0029% 0.91% 
B19 19.04 0.427 -0.00287 0.67 
B20 14.42 0.492 =0. 00232 0.97 ** 
B21 ae ae 0.068 -0.00064 Oas2 
B22 40.46 -0.040 0.00018 ead 
B23 Zoeko OUTS 0.00161 0.35 
B24 Pas rat O.325 =O. 99544 0.47 
B25 25.04 ean -0.00280 ong 
B26 8.14 0.008 -0.00034 0.36 
B27 a4<62 ~0.104 O..008Z7 0.48 
B28 ae 0.414** -0.00322%** 0.98% 
B29 20350 Os oul -0.00341 0.74 
B30 30.54 0.040 -0.00047 0.34 
B31 27.02 —0.261%* 0.00228 0.80 
B32 See coes Us 008 0.00070 0.04 
Bos hie) 0), O36 0.00065 0.34 
B34 48.20 0.086 -0.00104 O32 
B35 Siieee 0.067 =O. O00 7a Omang) 
B36 SUIS Ts) 0.55 -0 00323 0.87% 
B37 kone) 0.054 -0.00008 Oc 22 
B38 Sere a) 0.064 -0.00053 0.14 
B39 Spada C2202 ~O 50g 11 0.80 
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Table G-1. Barley Sites (cont.) 





Site b b, oR R? 
E01 Papo) ie 0.246 =O2 00185 0.88% 
E02 37-193 0.043 —0.,00072 Di. 72 
E03 a6.) --0.042 -0.00024 0,75 
E04 33854 OS + -0.00089 0.03 
E06 S200 0 0.096 = 00025 0.80 
E07 24.34 0.079 -0.00045 0.56 
E08 22.9 0.334** -0.00171% 0.98*x 
E09 34 03 0.142 -0.00058 0.98*x 
E10 23.46 0.430 -G 00330 o.o8 
E11 20.74 0.290** -0.00247%* 0.95% 
E13 36,20 0.148 —9 00079 0.38 
E14 £0.67 »-0.054 0.00028 0.45 
E15 25.59 0.130 =000132% 0.68 
E17 242935 0.036 0.00086 Go30 
E20 25252 02056 -0.00036 Oey ee 
E21 29.36 0.186 (0.00184 0.60 
E22 28.67 OVS Tree = =U F002 tae 0.88%*x 
E23 23.08 0.287** -0.00224** 0.95% 
E24 Ze naG 0.046 0.00080 0.61 
E25 35.46 02.262 = 00200 Posies 
E26 26.92 1. 140% ~0 00035 0.76% 
E27 PE. 32 0 2284% —0. 00221 0.83% 
E29 44,48 Oetetee -0.00089 O20 
E30 Boeese) > Oe 036 0.00028 Ow15 
E32 24.84 0.335** -0.00288x 0.83% 
E33 22.06 0.184** -0.00046 0.96** 
E34 Cae 0.047 -0.00067 Os. 
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Table G-1. Barley Sites (cont.) 


Site jen b, bo R? 
L0 1 eyes —~O.070 0.00144 Pats 
L03 See0G  —Ws08e 0.00222 5 a os 
LO4 25762 0.309 -0.00389 0.50 
LO5 Saeco —-~O2248 0.00433 0.88 
LO06 et Ish 0.598 = 20101 7 0.83 
L08 S30) ~0.t64 -0.00478 0.99 
L10 23630 0.472 -0.00589 0.97 
L11 23298 GSTs = Ovi 1.00% 
ie 28.28 On 2? = 00267 0.99 
V3 38.95 0.374 -0.00456 0.84 
L14 16.93 v.O%6 -0.01589 0.94 


L15 21.54 0.846 ~0.01122 0.99 
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w02 42.11 
Ww03 30.09 
wo4 Sik. 20 
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wW10 27.16 
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W13 45 523 
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Barley Sites (cont.) 
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R? 
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Table G-1. Barley Sites (cont.) 


Site bo b, b, R? 

TO1 20.67 0.400 -0.00405 0.86 
TO2 wavece  —0.080 0.00142 0.70 
T03 31.48 0.083 0.00005 0.77 
TO7 39.24 0.106 =OS0G0g 0.88 
TO8 28.58 0.144 -0.00074 1.00 
TO9 48.83 Oy255 0.00256 0.97 
12 52.66 0.009 0.00019 0.86 
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Table G-2. Rapeseed Sites 


Site bo b; b, Ke 
B4 1 6.67 0.093 -0.00073 0.68 
B42 veo 0.023 -0.00003 Os32 
B43 4.00 -0.034 0.00049 O23 
B44 9.48 eels -0.00163 Uso0 
B45 7206 0.030 -0.00034 0.18 
B46 Paare -O,0 TS 0.00036 Uae 
B47 Tea55 ears =O, 00271 0.64 
B48 S421 Ue236 =O 200490 Ue oe. 
B49 Seon 2 -O0.058 0.00101 0.81 
B50 Gard 0.092 -0 00073 0.74 
BS1 4.85 0.141 -0.00072 Cen 
B52 15.68 0.038 OV0OGT? Sad 7 
B53 F269 iat ane 0.00135 0.67 
B54 14.54 0.166 -0.00140 O56 
B55 11.14 GUE —O 200 hes Oet6 
B56 6535 0.214% =0 <QO nT Ge9s 
B57 reo O Ooh re -0.00138 0.95% 
B58 10.86 0.148 -0.00114 Vea9 
B59 ro 0.091% -0.00076* 0.92% 
B60 TU205 «60—=- 0.2046 0.0004 Gres9 
B61 cya ke 0.028% -0.00028 Ureks 
B62 Pahl 0.114 -0.0005 0.84 
B63 Una ay 0.144 =O 0.68 
B64 6.04 O2039 -0.00036 U226 
B65 Te TNS, 0.028 -0.00034 Dales 
B66 Geos 0.225% -0.00224 Cag 
B67 eave 0.242 =0,. 00207 0.81 
B68 6.38 0.191% =O.00736 0.90% 
B69 eke 0.048 =0.00035 0.42 
B70 PRR, 07035 =072,00029 MAAS 
B71 10,02 —0,0036 -0.00007 Ory 
B72 aco 0.098 -0.0007 0.88 
B73 ills 05.099 -0.00083 0.34 
B74 (Rey 0.166 -0.00114 0.65 
Bis Hen00s <0 055 =0'. 00073 Liss 
B76 Gos 0.171% -0.00144 0.83 
B77 19,40: -=0.092 -0.00085 Cot 


B78 8.64 0.311% -0.00265x 0.88 
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Table G-2. Rapeseed Sites (cont.) 


Site bo b, bs R? 
E37 1.40 0.140% -~G ..00135s 0.85 
E38 1h. 34 Got3s -0.00105 0.70 
E39 14.96 6.033 -0.0001 0.40 
E40 20.34 0.245% -0.00266* 0.83 
E41 Ohaus 0.038 =~ Couls 0.45 
E42 11.81 GO. 164 =. COL5S6 Shes Fes 
E44 14.73 0.039 —-O 00056 0230 
E45 MG ae Or Ca 0.00049 0.41 
E46 18.46 0.048 =0:, 00021 0.47 
E47 face? -—-O.019 0.00026 O07 
E48 gets} C207 =O..000 7-1 Oem fo 
E49 Poe eo OSG 0.00069 0.29 
E55 14.22 02055 -0.00062 0.28 
E57 14.28 0.398 =F 2.0002 Dele 
E58 1820.1 Ovd32 -0. 00007 rico 
E60 | oe at OSUBZ -0.0006 0.68 
E61 7.14 0. 184e2-.-0.001282 0.86** 
E62 eiaGo 0.055 —0'. 00.056 4 ets 
E63 6.62 0. 1044% —-0.00052 0.92% 
E64 15.74 0.033 =0,,00005 DIATE 
E65 10.29 0.160* = Oo sas O66 
E66 vee i 0.098 =) 00065. 0.56 
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Table G-2. Rapeseed Sites (cont.) 


Site bo 5 Be R? 

L48 9.47 Oe =) 00867 0.93 
L49 10.23 0.14 -0.00189 0.46 
L50 Peace -OV242 Des 2¢ 0.93 
L51 Zone.  —-O,300 0.00567 0.69 
L52 boe8a Ono -0.00244 0.68 
L53 mits) 0,034 —0 00033 0.84 


L54 Pada e208 -0.00089 0.86 
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Table G-2. Rapeseed Sites (cont.) 


Site Do b, b, R? 

w49 8.42 0.053 =0.00079 0.70 
Ww50 3.64 0.008 ~0,00025 0.60 
W52 11.00 -0.082+ 0.00078 0570 
W53 13.31 0.036 =O. 00035 OF39 
W54 obs he O.ts72e% 2=—0. COR 0.85% 





¢ Sasseqgan 
no} 28502 Basseg 
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Table G-2. Rapeseed Sites (cont.) 


Site Do b, b, R° 
THEY foews  —~Ol 055 -0.0004 ase 
720 12.84 0.047 -0.00056 0.42 
ee ped 8 he (0S2 -0.00061 0.81 
22 8.47 0.109** -0.0019*x 1.00%* 
T17 Gyagd ~—O.025 -0.00038 C509 
T18 re. OU 0.10 -0.00174 0.62 
219 wooo oOo isg 0.00123 o. 97 
T14 0% a? 0.046 -0.00033 94 
715 3.44 0.090 -0.00073 Uo 
T16 2.06 0.020 -0.00042 0.43 
T26 pray os 0.062 -0.00076 Onaa 
2 ara a .09 0.066 -0.00068 G.98 
T28 Bor 7 -02033 0.00018 0.23 
229 2.90 0.074 -0.00068 re 
T30 Poo (0. 0n4 -0.0003 0.70 
31 L2.a8 0.004 =0. 00075 0290 


T32. 14.19 -0.056 0.00065 0.31 





Table G-3. Wheat Sites 





Site bo by Be R? 

J01 16.50 0.066 -0.00032 0.94 
J02 howe 6 03 F077 -0.00063 0.97 
J03 12.34 0.038 0.00016 1.00% 
J04 2s2o5). -Oebn 0.00016 0.09 
J06 35236 0.082 af) 00732 0.59 
J07 23568 0.349 -0.00424 Use5 
J08 16.44 0.104 -0.00028 0.98 
J09 Py ig sie OR oe: -0.00044 0.20 
J10 63.700 e=09156 0.00088 0.57 


J11 NS ye 0.222% -0.00206 0.99% 
J12 21.86 0.128% -0.00140 0.98% 


J13 Phe eo FS 0.184 -0.00184 Gih2 
J14 16.26 0.494**x -0.00676x 0.98% 
J15 15.61 0. 104 -0.00016 ete 
J16 Pa Ni fe) 0.034 0.00098 1.00** 


J17 ioe 5) 0.361% -0.00492 0.96% 
J18 20.10 0.642** -0.00673** 1.00%* 


J19 21.44 0.131% =O 00 tay, 0.97% 
J20 Lie 2 O25 =JeU0005 Gere 
J22 ao0 9% 0. O90 0.00073 0.86 
J23 Zhao t, 0.069 0.00095 Dsis 
J24 18.05 0.223% -0.00029 1.00** 
J25 P2516 0.156% -0.00184 0.96% 
J26 b igeal oh, 0.130 -0.00067 ees 
J27 16.48 0.053 -0.00022 Jee 
J28 Pa Es 0.347% -O0055e2 0.96% 
J29 2OadU 0.164 =0.00216 0.46 
J30 S0e15 0.059 =) 00079 0.30 
J31 S2570 OSnAbS —aOug2 0261 
J32 28.64 0-46 -0.00098 0.64 
J33 20,000 —O,004 0.00149 0.42 
J34 202.20 0.004 —0 S00 156 Caac 
J35 Fie bes es tee -0.00343 0.83 
J36 22259 0.138 ~Om00365 0236 
J37 25.94 0.386 =O 00565 ey, 
J38 ee 0.003 0.00000 0.01 
J39 So.20 0.134 —0.00238 0.42 


J40 2960 0 Beam Bo A, -O05,.00159 0:.65 
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APPENDIX H 


Crop Response Calculations for Experimental Sites 
List of Abbreviations 
Y-max Calculated Maximum Site Yield (100 kg/ha) 
X-max Phosphate Rate for Y-max (kg P,0;/ha) 


Y-90%max 90% of Y-max (100 kg/ha) 
X-90%max Phosphate Rate for Y-90%max (kg P,0;/ha) 


Yield Increase (100 kg/ha) 
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Table H-1. Barley Sites 


Yield %Yield 
Site Y-max X-max X-90%max Y-90%max Increase Increase 


BO1 40.2 101 62 36.2 Sif 24.2 
B02 PA 76 26 26.4 337 13a 
BO3 tO.0 0 0 thes 0.0 0.0 
BO4 20.8 0 0 OEE: 0.0 0.0 
BO5 40.3 101 35 36.3 4.7 13.0 
BO6 ES) 0 0 10.8 O35 0.0 
B07 41.4 ns 24 eR, 4.9 322 
B08 SM re, 81 48 52.0 3220 6273 
BOS SS a2 40 569 8.1 24.3 
B10 36.6 74 32 3320 8.8 26.6 
B11 ahs D 87 37 42.7 eee 22er 
SM 35..8 76 43 32935 14.4 44,8 
BS 49.4 25 45 44.5 11.4 2547 
B14 vA! Pee 0 0 2962 0.0 0.0 
B15 043.3 re 36 48.9 toe0 30.6 
B16 42.6 74 24 36.3 aeU Rk 
B17 39.6 0 0 a7 e0 0.0 0.0 
B18 Sacks: 78 38 3536 iZe3 3443 
B19 ae 0 83 ae 3363 bas 41.8 
B20 43.2 101 da 36.9 24.2 6275 
B21 24a8 oo 0 +226 0.0 0.0 
B22 40.9 0 0 40.9 0.0 0.0 
B23 2550 0 0 23 95 0.0 0.0 
B24 36,7 oF 2] sre Seo 160% 
B25 3241 52 16 25.29 See tet 
B26 2 13 0 G23 0.0 0.0 
B27 a ge 0 0 3D 0S 0 au 0.0 
B28 SoU 73 38 RS ee Those 39.4 
B29 Jo ou) 64 30 oo) ou 30 
B30 Ae 8 48 0 oiled 0.0 0.0 
B31 27.6 0 0 eich = 0.0 0.0 
B32 aed & 0 0 2226 06.0 0.0 
aes 52. | 0 0 ciel 0.0 0.0 
B34 Baleas 46 0 49.1 0.0 0.0 
B35 it 0 0 Pal ong 0.0 0.0 
B36 ote 33 49 mB 19.8 38.3 
B37 EB on) 0 0 18.1 0.0 0.0 
B38 41.2 67 0 aon 0.0 0.0 
B39 49.4 104 34 44,5 eG 1a 
B40 40.8 87 6 36,7 1s 1.4 
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Table H-1. Barley Sites (cont.) 


; Yield %Yield 
Site Y-max X-max X-90%max Y-90%max Increase Increase 


E01 1 Re 74 28 31.8 aed Ye ae 
E02 acs 34 0 34.0 Go 0.0 
E03 Sono 0 0 S050 Oa) 0.0 
E04 38.0 69 0 32.2 O..0 0.0 
E06 40.9 101 4 36.8 eS: SoG 
E07 31: ae | ag 0 24.8 a..0 D0 
E08 41.8 101 of 28.6 14.0 Dike! 
E09 43.8 101 38 39.4 4.7 ies 
E10 aes aS 36 35.6 1B ef 8 
E11 au. 1 66 29 25.6 and 23.4 
E13 40.2 46 0 36.8 0.0 0.0 
E14 44,4 101 0 41.4 0.0 0.0 
E15 22.0 ae 6 26.6 0.6 2.4 
E17 SLRS: 0 0 Shi es 0.0 0.0 
E20 30.0 56 0 Phe’ 0.0 0.0 
E21 Spee 6 11 ot a | Ved aD 
E22 42.4 84 37 ore ee) Pee 
E23 2 ie id a3 30.4 629 Zant 
E24 Pe = 0 2ere0 0.0 0.0 
E25 ee 74 22 41.1 Darl Lan 
E26 £1 bids 73 8 2005 Ueg Ore 
E27 Ue Ne 58 ie a0 Sree Plies 
E29 49.6 74 0 4562 0.0 0.0 
E30 aes Oo 0 0 3140 0.0 0.0 
E32 a6). 2 65 2) a2 .0 Dad 22.0 
E33 a7 wD 101 68 a320 10c7 Silead 
E34 rhs eye | a4 0 O83 ee 0.0 0.0 





Table H-1. 

Site Y-max 
LO1 2576 
L03 34.2 
L04 33 wl) 
LO5 35.3 
L06 46.0 
L08 see 
L10 ae 63 
L11 3.5 
L12 i fe ate 
Bt3 48.3 
L14 a2 a9 


L15 


X-max 


0 
0 


Barley Sites (cont.) 


X-90%max 


0 
0 


Y-90%max 


Yield 
Increase 
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e © 
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*%Yield 
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Table H-1. Barley Sites (cont.) 


Yield %Yield 
Site Y-max X-max X-90%max Y-90%max Increase Increase 


WO 1 45.9 67 34 41.3 Sar 14.4 
WO02 43.1 1S 0 42.9 0.0 0.0 
WO3 oy 0 0 a 0.0 0.0 
wo4 62.5 67 20 3540 3.0 10.6 
WO5 She eee) 48 0 Biers 0.0 0.0 
WO06 = a 60 30 Zone oeo 3340 
W07 40.1 0 0 40.1 0.0 0.0 
w08 34.8 101 14 OAwS a 2.0 
wos S063 101 12 34.5 0.6 ae, 
W10 S05 69 0 Lhe eu 0.0 
Wi2 41.8 47 0 41.2 0.0 0.0 
W13 AS 94 0 46.0 0.0 0.0 
wi4 pe et 0 0 46.7 0.0 0.0 
W15 46.5 101 36 41.8 Gia? Lo.2 
W16 36.5 48 0 Saat 0.0 0.0 
W17 oO ye, 68 = 20a Ono no 
W18 41.4 66 2 BMS hate! 1052 
w19 41.0 0 0 41.0 0.0 0.0 
W20 B0</ 65 i S252 Fagen, 6.3 
W22 ey 64 0 Poh te 0.0 0.0 
W23 au She 10 46.8 ees 2.0 
W24 Sane 0 0 34.2 0.0 0.0 
W25 A201 0 0 4254 0.0 0.0 
W26 48.0 58 4 42.4 0.0 0.0 
W27 39.4 38 0 Bde ss 0.0 0.0 
W28 Gari, 26 0 cP ANe| 0.0 0.0 
W29 2oeU 82 0 2500 O20 0.0 
W31 ene 0 0 25.4 0.0 0.0 
W34 28.8 67 19 PAS ES, re iL Shwe 
W36 te, 60 2 Dilene 0.0 0.0 
W37 S2<U 78 0 Sa, 0.0 0.0 
W38 a 91 0 40.5 0.0 0.0 
wW41 S026 34a) 0 28.4 0.0 O50 
W42 2 3) 0 0 PAA 0.0 0.0 
wW43 49.4 101 25 44,5 Ul© ox 
w44 Le as 64 24 aa Ga 16.1 
W46 2aieet 67 24 21.4 i P| ie 
W47 PRE eM, 69 8 265) 0.8 beat 
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Table H-1. 

Site Y-max 
TO1 go 
TO02 36.5 
TO3 39.9 
TO7 PISS: 
TO8 3657 
TOS 56.9 
B10 5236 
Et2 55 3G 


X-max 


25 


Barley Sites (cont.) 


X-90%max 


Zo 
0 
45 
0 
3A 
i 
tS 
0 


Y-90%max 


53's 


NPN OW) ND) 


Yield 


Increase 


ON = WOW OO ~) 
e ° e e . 
ONAWMWOWO DO 
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%Yield 
Increase 


ine) 
OPN =—-OOCOO 0 


° e e 


—_d 
° ° . e 


e 
Os) 10 10 © '=1'© OF 





Table H-2. 


Site 


X-max 


72 
101 
0 
67 


On 


ONAN S| SH LSPWWNHWOEP DWE AHAWONAITOAINWMOD WN SAIN — OOW ON OO)N 


Y-max 


Rapeseed Sites 


X-90%max 


32 
60 

0 
35 


0 
19 
34 


y-90%max 


9. 


Yield 
Increase 


. e ° . e 


DOWOPODO—$-DOHKHAIVAPOOKFUTODOONWATOKP WH WANOANAWOO V10O = ND 
° . & . . ° ° 
WOPDWWMWODONODWAWAWMWOODOKHOAIWOH]-ANWAHAAAIADAH LH ODHNITUIGCDOOW 


%Yield 
Increase 
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—s 


—_ —— 
° e . e . 


WMWONDADWODNNNNDTDOVNA—|WwOO 


ip 


244 


SH sider 


% 


2 stoae 


ne 


_ 


~ 


“atte 


aDa-F 


teat @ 


oe Tirs 





¢ 


245 


Table H-2. Rapeseed Sites (cont.) 


Yield *Yield 


Site X-max Y-max xX-90%max Y-90%max Increase Increase 
E37 58 file nae a2 tice “ng 1&5 
E38 Bb (Ged 29 pr ee Se 2D hbase 
E39 101 PPok Zit hie ue Of 4.2 
E40 52 16.9 ao howe 4.7 3Q7 
E41 107 24.2 0 2acu 0.0 0.0 
E42 58 16.4 24 Toad yO 20st 
E44 39 15.8 0 15.0 O..0 0.0 
E45 0 16.5 0 16.5 0.0 0.0 
E46 101 Zea iy, 19.6 0.8 4.0 
E47 0 Ib 0 19.6 Ong 0.0 
E48 56 ioea 9 12.4 sel, 5.4 
E49 0 15.9 0 15.9 O70 0.0 
E53 40 745 0 6.9 0.0 0.0 
E54 39 Hy 3 0 Dine ORE, 0.0 
E55 49 15.9 0 14.4 0.0 0.0 
E57 101 16.8 jhe 15.1 0.6 Sal 
E58 101 20.9 12 18.8 Or 0.6 
E60 76 16.5 on 14.9 1.0 1005 
E61 8 1 leat | 45 ocr aye. 44.9 
E62 ate 13.4 3 Wi ob O52 5 toe 
E63 101 es 59 i hers 4.6 40.6 
E64 101 tao 25 L655 0.8 4,7 
E65 67 ISeg 31 14.3 300 26.6 
E66 84 LOR 43 Ove So 34.1 





Table H-2. 
Site X-max 
L48 39 
L49 4] 
L50 0 
L51 0 
L52 29 
L53 0 
L54 50 


Y-max 


12.0 
Lage 
iW aN: 
22.2 
16.0 
Les 
eo 


1 
1 


Rapeseed Sites (cont.) 


X-90%max 


DONDWDAOW — 


Y-90%max 


Yield 
Increase 
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%Yield 
Increase 


1 
1 





Table H-2. 


Site 


wW49 
W50 
W52 
W53 
W54 


X-max 


38 
18 

0 
=H 
59 


Y-max 


X-90%max 


Rapeseed Sites (cont.) 


Y-90%max 


Yield 
Increase 


%Yield 
Increase 


247 


Table H-2. Rapeseed Sites (cont.) 


Yield %Yield 
Site X-max Y-max X-90%max Y-90%max Increase Increase 
T19 30 537 0 14.4 0.0 0.0 
T20 24 14.2 0 13.9 0.0 0.0 
T21 29 228 0 12.0 0.0 OO 
T22 0 1092 0 rece 0.0 0.0 
47 0 6277 0 6.7 O20 0.0 
T18 20 18.8 0 18.3 0.0 CoO 
TT3 0 aa? 0 8.7 0.0 0.0 
T14 78 236 48 2.4 1.6 66.7 
T15 69 6.6 a0 5.9 2.5 41.5 
T16 Ze 9.5 0 9.2 0.0 0.0 
T26 0 1231 0 14.1 0.0 020 
T27 56 4.6 0 4.55 0.0 0.8 
T28 0 S79 0 3.9 0:6 0.0 
T29 90 al 0 8.8 0.0 0.0 
T30 z doe 0 12.2 ‘oPae) 0.0 
731 0 hese 0 12-52 0.0 0-0 
T32 0 14.4 0 14,4 0.0 0.0 


as 
, & €% 
\ = 
aint ae 
i= 
. 
= 
‘j 
Ed wv 


= ~~ ae 
wae 
~~ oy 
4 
ell 
vw 
1 . 
a 
J g 
| to 


. 
Te ’ 
> 
* 
> - 7 7 
~ * a 





249 


Table H-3. Wheat Sites 


Yield %Yield 
Site Y-max X-max X-90%max Y-90%max Increase Increase 


J01 20.6 ris 30 1526 Peo 9.6 
J02 Ziad a5 34 20.4 Fal 2) Lah 
J03 ZU 134 106 isi ate ats 3 
J04 Pe ea | 0 0 rae at | 0.0 0.0 
J06 37.4 hs 0 36.1 0.0 0.0 
J07 S251 46 17 2950 4.9 re 
J08 Pie ed 84 50 Pee 4.5 Zien! 
J09 abe ae 68 0 28.0 O20 0.0 
J10 34.3 0 0 34.3 Goo 0.0 
J11 20.6 60 20 a6 4.6 24 ag 
J12 PAS a 6 Za Caz 2.9 
J13 B059 56 tz Zee 2.1 has 
J14 20.1 4] 20 24.0 14 SUG 
J15 PN 67 4] Za G ao ie 
J16 31.0 67 44 2 Siac! nige 
J17 Pe A 4] 17 Pa i. 4.7 1352 
J18 SU pile 54 28 S340 Van3 59.4 
J19 26.0 63 13 23.4 A Go 
J20 22.6 67 21 20.4 2.4 neg 
J22 742 AM | 0 0 20a 0.0 0.0 
J23 fF tS 0 0 PAY Obs 0.0 a0 
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Effect Coding of Site Classification Systems 
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Table I-1. Barley Sites 


Agro-climatic Area 


Di D2 D3 Dé 
1 1 0 0 0 
2A 0 1 0 0 
2H 0 0 1 0 
3H 0 0 0 1 
3Ha my = = | = 
Soil Zone 
Di D2 D3 D4 
Gray 1 0 0 0 
Black 0 | 0 0 
Dark Gray 0 0 1 0 
Dark Brown 0 0 0 1 
Thin Black - | —1 —4 a 
Soil Order 
Die D2 oan Da 
Chernozemic 1 0 0 
Luvisolic 0 1 0 
Solonetzic 0 0 1 
Gleysolic = <a =i 
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Table I-2. Rapeseed Sites 


Agro-climatic Area 


Die 02 0D) D4 
1 1 0 0 0 
2A 0 1 0 0 
2H 0 0 1 0 
3H 0 0 0 1 
3Ha = son) se | aon 
Soil Zone 
Pi sb2e Dseey Ds DS 
Gray 1 0 0 0 0 
Black 0 1 0 0 0 
Dark Gray 0 0 1 0 0 
Dark Brown 0 0 0 1 0 
Brown 0 0 0 0 | 
Thin Black =| = at =i oa 
Soil Order 
Tig Dee es 
Gleysolic 1 0 0 
Luvisolic 0 1 0 
Solonetzic 0 0 1 
Chernozemic eit =| =f 
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Table I-3. Wheat Sites 


Agro-climatic Area 


D1 
1 1 
2A -1 
Soil Zone 
DI D2 
Black 1 0 
Dark Brown 0 1 
Thin Black - - 
Soil Order 
Di 
Chernozemic 


Solonetzic 
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Sab ; - APPENDIX J 


Calculation of Total Discriminatory Power 
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Total Discriminatory Power’ 


TDP =1- (N/ (N- K) (1 +&) + 1) 


Z 
il 


Total Sample Size 


mn 
W 


Number of Groups 


A = Eigenvalue 


"Reference: Tatsuoka, 1970 
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